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EXECUTIVE  SUMMARY 


The  HEV  toolbox  was  initially  conceived  as  a  “standard”  set  of  modules  for  use  in  evaluating 
the  performance  of  concept  HEV  architecture.  This  effort  resulted  in  the  development  of  a 
general-purpose  HEV  modeling  and  simulation  toolbox,  named  PATHS  (Performance 
Assessment  Toolbox  for  Hybrid  Systems).  PATHS  is  based  on  principles  of  physics  models, 
developed  in  Matlab®  and  Simulink®,  is  modularized  into  several  libraries  and  uses  the 
forward-looking  approach.  It  can  be  used  to  obtain  detailed  performance  characteristics  of 
the  modeled  HEV  or  its  sub-components,  as  well  as  perform  parametric  trade-off  studies. 
Particularly  useful  is  the  ability  to  size  components  and  evaluate  their  effects  on  vehicle 
performance  in  computer  simulation  prior  to  prototyping  or  fabricating  hardware.  The 
following  are  possible  applications  of  the  toolbox:  evaluation  of  new  drivetrain  options  / 
architectures;  test  and  debug  hybrid  vehicle  control  strategies;  sub-componentssizing  and 
specifications;  evaluation  of  hardware  modifications  on  HEV  performance;  impact  of 
software  modifications  on  HEV  performance  and  other  sub-components;  build  a  computer 
model  of  a  HEV  before  hardware  prototyping. 

The  following  is  the  current  list  of  libraries  within  PATHS:  engine  library,  electric  machine 
library,  energy  storage  library,  vehicle  library,  controllers  library,  couplings  library,  and  a 
miscellaneous  library.  The  multiple  components  within  a  library  have  the  same  input-output 
structure  so  that  they  can  be  interchanged  in  an  HEV  model.  The  user  can  also  add  to  the 
library  as  long  as  the  input-output  structure  of  the  particular  library  is  maintained.  PATHS  is 
folly  customizable  as  long  as  the  user  adheres  to  the  input-output  structure  of  each  library 
and  is  aware  of  the  existing  naming  convention  of  the  variables. 

Validation  of  PATHS  was  performed  by  first  developing  a  computer  model  of  a  22-foot 
series  hybrid  shuttle  bus,  followed  by  comparing  its  output  to  data  collected  from  the  same 
bus  on  a  chassis  dynamometer.  The  hybrid  shuttle  bus  was  placed  on  a  heavy-duty  chassis 
dynamometer  and  driven  by  a  human  operator,  while  data  was  collected  from  the  drive 
motors,  Auxiliary  Power  Unit  (APU),  batteries,  and  vehicle  speed  sensors.  This  data  was 
compared  with  model  predictions,  executing  the  same  driving  profile  as  that  of  the  hybrid 
bus.  Both  steady  state  and  transient  comparisons  were  made  between  the  bus  and  the  model. 
Fuel  consumption  predicted  by  the  APU  model  was  within  1%  over  6  EPA  cycles  of 
measured  data.  The  battery  state  of  charge  predicted  by  the  lead-acid  battery  model  was 
within  5%  over  6  EPA  cycles  of  measured  data.  The  motor  model  predictions  (transient  / 
steady  state)  are  within  10%  of  measured  data  for  low  speeds  and  lower  throttle  positions. 
Larger  deviations  between  the  measured  and  model  predicted  data  at  larger  speeds  and  higher 
throttle  positions.  The  dynamic  vehicle  model  validated  to  within  2%  of  measured  data. 

This  report  is  divided  into  the  following  sections:  User’s  Guide,  Technical  Notes,  and 
PATHS  Validation. 
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1.0  OVERVIEW  OF  DARPA  HEV  SIMULATION  TOOLBOX 


1.1  Background 


The  DARPA  HEV  Simulation  Toolbox  was  initially  conceived  as  a  “standard”  set  of 
modules  for  use  in  evaluating  the  performance  of  concept  Hybrid  Electric  Vehicle  (HEV) 
architectures  in  simulated  driving  conditions.  Two  versions  of  the  toolbox  are  in  existence: 
one  for  commercial  applications,  and  one  for  military  applications.  Both  toolboxes  contain 
the  same  modules.  However,  the  simulation  toolbox  for  military  applications  contains  an 
additional  library  of  military  specific  modules.  The  HEV  simulation  toolbox  was  collectively 
developed  by  the  following  organizations: 


Southwest  Research  Institute  (SwRI) 

Univ.  of  TX  Inst,  for  Adv.  Technology  (IAT) 
Georgia  Institute  of  Technology  (GT) 

Univ.  of  Hawaii  (UH) 

Univ.  of  TX  Center  of  Electro-Mechanics  (CEM) 


Dr.  A.  Nedungadi 
Dr.  Scott  Fish 
Dr.  Thomas  Habetler 
Dr.  Vassilis  Syrmos 
Mr.  John  Kajs 


(210)  522-3965 
(512)  471-9060 
(404)  894-8929 
(808)  956-3432 
(512)471-4496 


1.2  Introduction 


This  user’s  guide  provides  a  concise  description  of  the  commercial  version  of  the  DARPA  HEV 
simulation  toolbox.  Use  of  this  toolbox  is  restricted  to  licensees  of  the  software.  An 
application  for  a  license  to  the  commercial  version  of  the  simulation  toolbox  is  available  by 
contacting  SwRI  (Dr.  A.  Nedungadi  (210)  522-3965).  All  applications  for  a  license 
agreement  to  use  the  military  version  of  the  simulation  toolbox  must  be  addressed  to  IAT 
(Dr.  Scott  Fish  (512)  471-9060). 


The  simulation  toolbox  is  organized  into  a  suite  of  libraries,  each  containing  modular  computer 
models  of  typical  HEV  components.  All  the  models  were  developed  using  MATLAB®  and  its 
graphical  click-and-pick  modeling  tool,  Simulink®.  The  objective  of  the  developed  toolbox  is  to 
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provide  the  user  the  essential  building  blocks  to  simulate  and  analyze  the  performance  of  a 
custom  built  or  pre-existing  electric  hybrid  vehicle. 


□ 

Engines 
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Electric 

Machines 


□ 


Energy 

Storage 


Q 


□ 


□ 


O 


Vehicles  Controllers  Coupling  Miscel 

aneous 


Figure  1.  Hybrid  Electric  Vehicle  Simulink  Toolbox  (Commercial) 


The  following  list  describes  all  the  libraries  and  models  available  within  the  toolbox: 

ENGINE  LIBRARY 

I  |  Table  lookup  diesel  engine  without  emissions 
I  |  Table  lookup  gas  turbine  engine  without  emissions 
I  |  Spark  ignition  dynamic  engine  with  emissions 

ELECTRIC  MACHINE  LIBRARY 

I  |  AC  induction  steady-state  motor 
I  |  AC  induction  dynamic  motor 
I  |  Permanent  magnet  steady-state  generator 
I  [  Permanent  magnet  dynamic  generator 
I  |  Table  lookup  steady-state  generator 

ENERGY  STORAGE  LIBRARY 

|  [  Lead  acid  battery 
I  |  flywheel  battery 
|  |  Capacitor 


4 


VEHICLE  LIBRARY 

I  |  Tracked  vehicle 
I  I  Wheeled  vehicle  (version  1  and  2) 

I  |  Wheel  (version  1  and  2) 

I  |  Driver 

CONTROLLERS  LIBRARY 

I  |  Energy  management  (2) 

•  Energy  storage  system  controller 

•  APU  controller 
I  |  System  controller 

COUPLINGS  LIBRARY 
|~1  Gear  box 

I  I  Engine-generator  coupling 
I  |  Motor-wheel  coupling  (version  1  and  2) 

MISCELLANEOUS  LIBRARY 
I  I DC-DC  converter 
I  |  AC- AC  converter 
I  |  AC-DC  converter 
I  I  DC-AC  converter 
|  |  Capacitor 
[~1  Roller  bearing 
I  I  Magnetic  bearing 
|~1  Windage  model 
□  Driving  profile 

A  rigid  naming  convention  was  used  for  all  the  variables  used  within  the  developed  models. 
The  first  three  alphabetic  characters  represent  the  abbreviation  of  the  model  name.  The  fourth 
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single  alphabetic  character  represents  the  possible  device  type,  if  one  exists.  The  fifth  single 
numeric  character  between  0  and  9  represents  the  number  of  instances  of  this  item.  The  sixth 
character  is  a  required  underscore,  followed  by  the  variable  name.  Table  1  summarizes  all  the 
model  names,  the  respective  three  character  abbreviations  that  describe  the  particular  model 
name,  and  all  possible  device  types  for  each  model. 


|  Table  1 .  Summary  Of  Model  Names,  Abbreviations  And  Device  Types  [ 

MODEL 

ABBREVIATION 

POSSIBLE  DEVICE  TYPES 

Engine 

eng 

s  (spark  ignition) 
d  (diesel) 
t  (turbine) 
s  (Stirling) 

Motor 

mot 

I  (induction  -  steady-state) 
k  (induction  -  dynamic) 

Generator 

gen 

p  (permanent  magnet  -  steady-state) 
q  (permanent  magnet  -  dynamic) 

Energy  Storage  System 

ess 

b  (lead  acid  battery) 
f  (flywheel) 
c  (capacitor) 

Axle 

axl 

none 

Coupling 

cpI 

none 

Wheel 

whl 

none 

Vehicle 

veh 

none 

Transmission 

trm 

m  (manual) 
a  (automatic) 

Controller 

ctl 

none 

For  example,  the  variable  motIl_VAr  represents  the  first  instance  of  the  variable  name  VAr  of 
the  induction  motor  model.  The  variable  genp2_J  represents  the  second  instance  of  the  variable 
name  J  of  the  permanent  magnet  generator  model.  MATLAB®  is  case  sensitive,  therefore 
motll  VAr  is  not  the  same  as  motil_VAr. 


The  following  format  is  used  to  describe  each  model.  A  table  summarizes  the  required  inputs 
and  outputs  of  the  model,  followed  by  a  list  of  constants  that  are  used  within  the  model.  Any 
auxiliary  calculations  that  are  not  part  of  the  output  are  listed  in  the  section  of  auxiliary 
calculations.  This  is  followed  by  the  first  three  lines  of  the  mask  (consult  your  Simulink 
manual)  of  the  model.  This  is  provided  to  enable  the  user  to  recreate  the  mask  if  it  is 
accidentally  deleted. 
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A  list  of  assumptions  is  listed  for  each  model,  enabling  the  user  to  determine  if  the  complexity 
of  the  model  would  suffice  for  the  given  application.  Finally,  typical  values  of  the  constants 
used  within  each  of  the  models  are  presented  to  give  the  user  an  indication  of  the  appropriate 
order  of  magnitude  of  the  constants  used  during  simulation  of  the  different  models. 

This  user’s  guide  is  accompanied  by  a  technical  notes  document  that  contains  all  the  differential 
equations/equations  used  within  each  of  the  models.  These  equations  are  provided  to  enable  the 
user  to  assess  the  applicability  of  the  particular  model  to  the  application  in  question.  If  certain 
modifications  or  extensions  are  required,  the  user  will  have  to  retrieve  the  particular  model  from 
the  Simulink®  environment  and  make  the  required  changes. 

It  is  assumed  that  the  user  is  an  experienced  Simulink  user  since  many  of  the  developed  models 
use  similar  advanced  features.  It  is  also  assumed  that  the  user  has  working  knowledge  of 
system  modeling  and  interfacing  modules  with  multiple  inputs  and  outputs.  The  user  interacts 
with  the  developed  toolbox  by  simply  dragging  the  required  item  from  the  appropriate  toolbox 
library  and  connecting  it  with  other  modules  while  adhering  to  the  input  and  output  connections 
described  in  this  manual. 

The  developed  models  have  all  been  tested  through  numerous  computer  simulations.  The 
equations  used  to  model  the  different  components  were  derived  from  classical  text  books,  using 
sound  engineering  methods.  Validation  of  the  most  essential  components  has  been  completed 
(see  accompanying  validation  report).  It  is  possible  to  tweak  the  parameters  used  within  each 
model  to  make  its  performance  compare  to  the  physical  system.  This  calibration  task  is  left  to 
the  user  if  he/she  needs  to  calibrate  a  particular  model  to  a  known  physical  system.  The 
numerous  simulations  performed  during  the  verification  phase  indicate  that  the  simulation  of  a 
complex  system,  such  as  an  electric  hybrid  vehicle,  can  be  significantly  affected  by  either 
parameter  mismatch  within  a  model  or  intercomponent  mismatch.  Typically,  this  is  diagnosed 
by  the  simulation  taking  very  small  time  steps  or  simulation  instabilities. 
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2.0  ENGINE  LIBRARY 


The  engine  library  contains  three  types  of  engines:  (1)  table  lookup  steady-state  diesel  engine 
without  emissions,  (2)  table  lookup  steady-state  gas  turbine  engine  without  emissions,  and 
(3)  a  dynamic  model  of  a  spark  ignition  engine  with  emission  calculations.  In  an  HEV 
application,  the  engine  model  is  typically  connected  to  a  generator  using  a  mechanical 
coupling  from  the  coupling  library  (see  Section  6).  The  first  two  engine  models  incorporate  a 
lag  response  time  with  steady-state  performance  tables.  Their  use  maintains  the  shape  of  the 
torque  vs.  speed  and  throttle  setting  function  while  scaling  for  the  desired  shaft  speed  and 
torque  at  design  conditions.  The  third  engine  model  is  a  more  detailed  dynamic  model  of  a 
spark-ignition  engine. 


>  > 

[ 

^  i 

►  > 

s 

► 

S'  ► 

_ _ ] 

>  / 

►  > 

✓ 

► 

V1.0  > 

► 

diesel 

gas  turbine  | 

engine  engine  si  Engine 

(Version  1) 


Figure  2.  Engine  Model  Library 


2.1  Diesel  Engine 


diesel  engine 


INPUTS 

OUTPUTS 

Non-Dimensional 
Throttle  Control  (0-1) 

Vector  of  Engine  Torque  (N-m) 
and  Inertia  Constant  (kg*mA2) 

Engine  Shaft  Speed  (RPM) 

Fuel  Mass  (kg) 

Purpose:  The  purpose  of  the  diesel  engine  model  is  to  accurately  represent  the  engine  torque 
and  fuel  consumption  characteristics  to  its  input  throttle  setting  and  shaft  speed.  The  model 
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is  quasi-steady  with  a  user-defined  first  order  lag  in  response  to  the  throttle.  Other  than  the 
engine  lag,  the  performance  is  determined  by  table  lookup  with  steady-state  performance 
data.  The  fuel  consumption  data  is  in  lb/(hp  hr)  and  is  assumed  not  to  vary  with  engine 
power  rating.  The  torque  table  is  scaled  with  the  engine  power  and  speed  to  remain 
consistent. 


Parameters:  Note:  These  values  are  loaded  from  disk  through  the  mask  or  inputted  by  the 
user. 


engdl_pwr  engine  design  power  (hp) 

engdl_spd  engine  design  speed  (rpm) 

engdlJ  inertia  constant  (kg*mA2) 

fuelmassi  initial  fuel  mass  (kg) 


sample_N  max  number  of  samples 

sample_T  simulation  time  between  samples 


Auxiliary  Calculations:  None 


Diesel  Engine  Mask 


(1) ENGD 

(2)  Diesel  Engine  Block  Definitions\n(Click  on  "Help"  to  list  Inport  and  Outport 
Parameters.)|[pwr(hp),spd(rpm),J(kgm2)]|Initial  Fuel  Mass  (kg)|[sample_N,sample_T] 

(3)  load('engd');pwr=@l(l)*746;spd=@l(2)*pi/30;J=@l(3);fuel_i=@2; 
sample_N=@3  ( 1 ) ;  sample_T-@3  (2); 

(4)  Drawing  Command. 

(5)  Diesel  Engine  Block  Port  Parameters.\n\nInport  1 :  Nondimensional  command  signal 
(throttle  setting)  from  diesel  engine  controller  blockAnlnport  2:  Speed  from  the  generator 
coupler  block.\n\nOutport  1 :  Vector  of  [ENGTT,  ENGT  J]  output  to  the  generator 
coupler  block. 


Model  Equations:  The  diesel  model  here  is  composed  of  three  parts.  The  first  part  applies 
a  1/ts  transfer  function  to  the  throttle,  modeling  the  response  time  lag  in  the  actual  machine 
to  inputs  from  the  throttle.  The  second  two  parts  determine  the  output  torque  and  specific 
fuel  consumption  from  steady-state  lookup  tables  that  are  scaled  from  a  known  diesel  engine 
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performance  set.  It  is  assumed  that  the  steady-state  performance  tables  are  adequate  when 
combined  with  the  engine  lag  to  describe  the  actual  performance  of  a  dynamic  model.  No 
calculations  are  included  for  emissions  in  this  model. 


2.2  Gas  Turbine  Engine 


INPUTS 

OUTPUTS 

Non-Dimensional 

Throttle  Control  (0-1) 

Vector  of  Engine  Torque  (N-m) 
and  Inertia  Constant  (kg*mA2) 

Engine  Shaft  Speed  (RPM) 

Fuel  Mass  (kg) 

Purpose:  The  purpose  of  the  gas  turbine  model  is  to  accurately  represent  the  engine  torque 
and  fuel  consumption  characteristics  to  its  input  throttle  setting  and  shaft  speed.  The  model 
is  quasi-steady  with  a  user-defined  first  order  lag  in  response  to  the  throttle.  Other  than  the 
engine  lag,  the  performance  is  determined  by  table  lookup  with  steady-state  performance 
data.  The  fuel  consumption  data  is  in  lb/(hp  hr)  and  is  assumed  not  to  vary  with  engine 
power  rating.  The  torque  table  is  scaled  with  the  engine  power  and  speed  to  remain 
consistent.  The  output  fuel  mass  represents  remaining  mass  of  fuel  in  tanks  connected  to  this 
engine.  This  data  can  be  used  to  account  for  weight  loss  in  time  by  vehicle  dynamics  models. 

Parameters:  Note:  These  values  are  loaded  from  disk  through  the  mask  or  inputted  by  the 
user. 


engtl__pwr 

engtl_spd 

engtlJ 

fuel_mass_i 


engine  design  power  (hp) 
engine  design  speed  (rpm) 
inertia  constant  (kg*mA2) 
initial  fuel  mass  (kg) 


sample_N  max  number  of  samples 

sample  T  simulation  time  between  samples 


Auxiliary  Calculations:  None 
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Gas  Turbine  Engine  Mask 


(1)  ENGT 

(2)  Gas  Turbine  Engine  Block  Definitions\n(Click  on  "Help"  to  list  Inport  and  Outport 
Parameters.)!  [pwr(hp),spd(rpm),J(kgm2)]|Initial  Fuel  Mass(kg)|  [sample_N,sample_T] 

(3)  load('engtdata');pwr=@l  (1  )*746;spd=@l  (2)*pi/30;  J=@l  (3);fuel_i=@2; 
sample_N=@3(l);sample_T=@3(2); 

(4)  Drawing  Commands 

(5)  Gas  Turbine  Engine  Block  Port  Parameters.\n\nInport  1:  Throttle  Command 
(0->l).\nlnport  2:  Shaft  Speed  (RPM).\n\nOutport  1 :  Vector  of  [Shaft  Torque 
(N*m), Inertia  Constant  (kg*mA2)]. 

Model  Equations:  The  gas  turbine  model  here  is  composed  of  three  parts.  The  first  part 
applies  a  1/xs  transfer  function  to  the  throttle,  modeling  the  response  time  lag  in  the  actual 
machine  to  inputs  from  the  throttle.  The  second  two  parts  determine  the  output  torque  and 
specific  fuel  consumption  from  steady-state  lookup  tables  that  are  scaled  from  a  known  diesel 
engine  performance  set.  It  is  assumed  that  the  steady-state  performance  tables  are  adequate, 
when  combined  with  the  engine  lag,  to  describe  the  actual  performance  of  a  dynamic  model. 
No  calculations  are  included  for  emissions  in  this  model. 

2.3  Spark  Ignition  Engine 


Universal  gas  constant  [J/kmol/K] 

Standard  pressure  [  Pa] 

Standard  temperature  [K] 

Molecular  weight  of  the  air  [kg/kmol] 

Specific  gas  constant  for  air  [J/kg/K] 

Density  of  standard  air  [  kg/m3] 

Molecular  weight  of  exhaust  gases  [  kg/kmol] 
Specific  gas  constant  for  exhaust  gases  [  J/kg/K] 


INPUTS 

OUTPUTS 

Throttle  Fraction  [0-1] 

Output  Torque  [Nm] 

Desired  Engine  Speed  [rpm] 

Engine  Inertia  [kgm2] 

Engine  ON/OFF  Switch  [0/1] 

Total  Fuel  Consumed  [kg] 

%  General  constants: 

eng_universal_gas_const: 

eng_std_pres: 

eng_std_temp: 

eng_mol_wt_air: 

eng_gas_const_air: 

engstdairdensity: 

eng_mol_wt_exhaust: 

eng_gas_const_exhaust: 
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%  Primary  engine  scaling  parameters 

engdisplacement:  Engine  displacement  [m3] 

eng_rated_speed:  Maximum  engine  rated  speed  [rad/s] 


%  Secondary  engine  scaling  parameters 


engcompressionratio : 

engratedbmep: 

engratedtorque: 

eng_rated_power: 

eng_rated_vol_eff: 

engratedairvolflow: 

eng  rated  air  mass  flow: 


Engine  compression  ratio  [  ] 

Brake  mean  effective  pressure  at  rated  conditions  [Pa] 
Rated  engine  torque  [N-m] 

Rated  engine  power  [watts] 

Volumetric  efficiency  at  rated  conditions  [fraction] 
Air  volume  flow  rate  at  rated  conditions  [  m3/s] 

Air  mass  flow  rate  at  rated  conditions  [  kg/s] 


%  SI  Engine  2  (Main  model) 


eng_ambient_pressure : 
eng_ambient_temperature : 
eng_ambient_air_density: 
enginertia: 
engofftorquefactor: 


Ambient  pressure  [Pa] 

Ambient  temperature  [K] 

Density  of  ambient  air  [kg/m3] 

Engine  rotary  inertia  [kg-m2] 

Factor  for  engine  torque  when  engine  is  off  [N-m  /  rad/s] 


%  Air  fuel  system  (1) 

eng  manifold  volume  ratio:  Ratio  of  intake  manifold  to  engine  displacement  [  ] 
eng  manifold  volume:  Volume  of  intake  manifold  [m3] 


%  Air  system  (1.1) 


engthrottlebypassratio : 

eng_rated_intake_vacuum: 
eng_rated_pres_ratio_infl : 

eng_max_air_flow: 

eng_speed_density_weight: 


Ratio  of  bypass  flow  area  around  throttle  to  total  throttle  area 
[fraction  ] 

Intake  manifold  vacuum  at  rated  conditions  [Pa] 

Pressure  ratio  influence  factor  (in  flow  equation)  at  rated 
conditions  [fraction  ] 

Maximum  air  flow  through  throttle  (choked)  [kg/s] 

Fractional  weight  for  speed  density  air  flow  compared  to  mass 
flow  sensor  [fraction] 


%  Fuel  system  (1.2) 

engafratiospeed: 
eng  af  ratio  Pi : 
eng_af_ratio_table: 


Vector  of  engine  speed  for  air  fuel  ratio  [rad/s  ] 

Vector  of  intake  manifold  pressure  for  air  fuel  ratio  table  [Pa]; 
Air  fuel  ratio  table  [number] 
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%  Mixing  and  dynamics  (1.3) 

eng_liqmd_fuel_fraction :  Fraction  of  fuel  that  stays  in  liquid  form  [fraction] 

eng_fuel_evaporation_time_const:  Time  constant  for  evaporation  of  liquid  fuel  in  intake 

manifold  wall  [sec] 

%  Reciprocator  air  flow  (2.1) 

eng_vol_eff_speed:  Vector  of  engine  speed  for  volumetric  efficiency  table  [rad/s] 

eng_vol_eff_Pi:  Vector  of  intake  manifold  pressure  for  volumetric  efficiency 

table  [Pa] 

engvoleffteble:  Volumetric  efficiency  table  [number] 


%  Torque  production  (2.2) 


engnominalisfc: 
engfmepfactorO : 
engfinepfactorl : 
eng_fmep_factor_2 : 


(Nominal)  indicated  specific  fuel  consumption  [kg/J] 
Correlation  factor  0  for  friction  mean  effective  pressure  [  ] 
Correlation  factor  1  for  friction  mean  effective  pressure  [  ] 
Correlation  factor  2  for  friction  mean  effective  pressure  [  ] 


%  Engine  exhaust  (2.3) 

engstoichafratio: 

eng_exh_temp_af_coef : 

engexhtempspeed: 

engexhtempPi: 

eng_exh_temp_table: 

engemissionsspeed: 

engemissionsPi: 

eng_emissions_HC_table: 
eng_emissions_CO_table: 
eng_emissions_NOx_table : 
engemissionsPMtable: 


Stoichiometric  air  fuel  ratio  [  ] 

Correlation  factor  for  effect  of  air  fuel  ratio  on  exhaust  temp.[  ] 
Vector  of  engine  speed  for  exhaust  temperature  table  [rad/s] 
Vector  of  intake  manifold  pressure  for  exhaust  temp,  table  [Pa] 
Exhaust  temperature  table  [K] 

Vector  of  engine  speed  for  engine  emission  tables  [rad/s] 

Vector  of  intake  manifold  pressure  for  engine  out  emission 
tables  [Pa] 

Engine  out  hydrocarbon  emissions  table  [kg/kg  fuel] 

Engine  out  carbon  monoxide  emissions  table  [kg/kg/fuel] 
Engine  out  oxides  of  nitrogen  emissions  table  [kg/kg  fuel] 
Engine  out  particulate  matter  emissions  table  [kg/kg  fuel] 


%  Muffler  and  aftertreatment  (3) 


eng_rated_exhaust_pressure : 
eng_exh_pres_scaling: 
eng_rated_exhaust_temperature: 
engratedexhaustdensity : 
eng_muffler_pressure_factor: 
engcatalysteffNOx: 
engcatalysteffHC : 
engcataly  st_eff_CO : 
engcatalysteffPM: 


Exhaust  pressure  at  rated  conditions  [Pa] 

Scaling  factor  for  matching  exhaust  pressure  [  ] 
Exhaust  temperature  at  rated  conditions  [K] 

Density  of  exhaust  at  rated  conditions  [kg/m3] 
Correlation  factor  for  engine  muffler  pressure  [  ] 
Catalyst  conversion  efficiency  for  oxides  of  nitrogen  [  ] 
Catalyst  conversion  efficiency  for  hydrocarbons  [  ] 
Catalyst  conversion  efficiency  for  carbon  monoxide  [  ] 
Catalyst  conversion  efficiency  for  particulate  matter  [  ] 
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The  following  variables  (last  value)  are  saved  in  work  space  for  post  run  analysis: 


eng_s_speed: 

eng_s_throttle: 

eng_torque: 

eng_total_emissions : 

eng  s  volumetric  efficiency: 

engs jDunipingtorque : 

eng_s_gross_torque : 

engscombustiontorque: 

eng_therm_eff_gross: 

engsjfrictiontorque: 

eng_s_exhaust_temp : 

eng_s_exhaust_gage_pressure: 

eng_s_exhaust_pressure : 

eng_engine_emissions: 

engvehicleemissions: 


Engine  speed  [rad/sec] 

Throttle  [fraction  ] 

Engine  torque  [Nm] 

Vector  of  total  engine  emissions  [g  ] 

Volumetric  efficiency  [number  ] 

Pumping  loop  torque  [Nm] 

Gross  torque  [Nm] 

Combustion  torque  [Nm] 

Gross  thermal  efficiency  [fraction] 

Friction  torque  [Nm] 

Exhaust  temperature  [K] 

Exhaust  gage  pressure  [Pa] 

Exhaust  absolute  pressure  [Pa] 

Vector  of  engine  out  emissions,HC,CO,NOx,PM  [g/hr] 
Vector  of  vehicle  out  emissions,HC,CO,NOx,PM  [g/hr] 


Mask:  None 


Table  2  presents  typical  values  used  for  the  constants  in  the  engine  model 


j  Table  2.  Typical  Constant  Values  in  Engine  Model 

Constant 

Value 

eng_universal_gas_const: 

8314.3  [J/kmol/K] 

eng_std_pres: 

101.3  [Pa] 

eng_std_temp: 

298.17  [K] 

eng_mol_wt_air: 

28.96  [  kg/kmol] 

eng_gas_const_air: 

287.096  [J/kg/K] 

eng_std_air_density: 

1.1834  [  kg/m3] 

eng_mol_wt_exhaust: 

28.96  [  kg/kmol] 

eng_gas_const_exhaust: 

287.096  [  J/kg/K] 

eng_displacement: 

0.0018  [m3] 

eng_rated_speed: 

628.3185  [  rad/s] 

eng_compression_ratio: 

8 

eng_rated_bmep: 

1 E6  [Paj 

eng_rated_torque: 

1.131E4  [Nm] 

eng_rated_power: 

7.106E6  [watts] 

eng_rated_vol_eff: 

80  [  %] 

eng_rated_air_vol_flow: 

0.072  [  m3/s] 

eng_rated_air_mass_flow: 

0.0852  [  kg/s] 

eng_ambient_pressure : 

101300  [Pa] 

eng_ambient_temperature: 

298.17  [K] 

eng_ambient_air_density: 

1.1834  [kg/m3] 

engjnertia: 

0.1454  [kg-mA2] 

eng_off_torque_factor: 

-0.1  [N-m/ rad/s] 

eng_manifold_volume_ratio: 

0.0018 

eng_manifold_volume: 

0.0018  [m3] 
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%  Air  system  (1.1) 

eng_throttle_bypass_ratio:  0.005 

engratedintakevacuum:  3000  [Pa] 

eng_rated_pres_ratio_infl:  0.234 
eng_max_air_flow:  0.3642  [kg/s] 

eng_speed_density_weight:  1 

engafratiospeed:  [0.1667, 1] 

eng  af  ratio  Pi:  [0,  88000, 102000]; 

eng_af_ratio_table:  [14.6, 14.6,13] 

[14.6, 14.6, 12.5] 

eng_liquid_fuel_fraction :  0.5 

engfuelevaporationtimeconst:  0.25  [sec] 

eng_vol_eff_speed:  [0.1667  0.3333  0.5000  0.6667  0.8333  1.0000] 

engvoleffPi:  [20000  100000]  [Pa] 

eng_vol_eff_table:  [0.5000  0.7000 

0.5000  0.7500 

0.5000  0.8000 

0.5000  0.8900 

0.5000  0.8700 

0.5000  0.7700] 

eng_nominal_isfc:  6.67E-8  [kg/J] 

engfmepfactorO:  0.75 

eng_fmep_factor_l :  0. 1 5 

eng_fmep_factor_2:  0.05 

engstoichafratio:  14.6 

eng_exh_temp_af_coef :  0 

eng_exh_temp_speed:  [0.1667  1.0000  ]  [rad/s] 

eng  exh  temp  Pi:  [20000  100000]  [Pa] 

engexhtemptable:  [673  823] 

[773  1173]  [deg  K] 

eng_emissions_speed:  [673  823  773  1173]  [rad/s] 

eng  emissions  Pi:  [40000  50000  60000  70000  80000  ]  [Pa] 

engemissionsHCtable:  [0.0162  0.0146  0.0135  0.0124  0.0113] 

[0.0151  0.0139  0.0134  0.0121  0.0111] 

[0.0140  0.0136  0.0130  0.0120  0.0110] 

[0.0130  0.0133  0.0128  0.0118  0.0107] 

[0.0124  0.0131  0.0125  0.0112  0.0100] 

[kg/kg  fuel] 

engemissionsCOtable:  [0.0802  0.0831  0.0851  0.0837  0.0855] 

[0.0840  0.0860  0.0890  0.0887  0.0898] 

[0.0861  0.0894  0.0935  0.0941  0.0933] 

[0.0892  0.0921  0.0969  0.0987  0.0987] 

[0.0909  0.0959  0.0995  0.1013  0.1016] 

[kg/kg  fuel] 
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engemissionsNOxtable: 


engemissionsPMtable : 


[0.0279 

0.0465 

0.0584 

0.0560 

0.0517] 

[0.0384 

0.0553 

0.0633 

0.0583 

0.0539] 

[0.0482 

0.0635 

0.0673 

0.0598 

0.0558] 

[0.0569 

0.0698 

0.0695 

0.0608 

0.0571] 

[0.0632 

0.0755 

0.0715 

0.0631 

0.0593] 

[kg/kg  fuel] 

[0  0 

0  0 

0] 

[0  0 

0  0 

0] 

[0  0 

0  0 

0] 

[0  0 

0  0 

0] 

[0  0 

0  0 

0] 

[kg/kg  fuel] 


eng_rated_exhaust_pressure:  40000  [Pa] 

eng_exh_pres_scaling:  0.5 

eng_rated_exhaust_temperature:  1173  [K] 

eng_rated_exhaust_density:  0.4196  [kg/m3] 

eng_muffler_pressure_factor:  4.85E5 

eng_catalyst_eff_NOx:  0.5 

eng_catalyst_eff_HC:  0.5 

eng_catalyst_eff_CO:  0.5 

eng_catalyst_eff_PM:  0.5 


Assumptions: 


The  following  assumptions  have  been  used  for  the  quasi  steady-state  model: 

■  May  be  used  with  any  type  of  engine: 

Spark  ignited  (SI)  liquid  fuel  (gasoline) 

Spark  ignited  (SI)  gaseous  fuel  (natural  gas) 

Compression  ignited  (Cl),  typically  called  diesel  engines 
Gas  turbine 
Any  other 

■  Uses  steady-state  maps  for  torque  curves,  fuel  consumption  and  emissions. 

■  Any  torque  in  the  range  of  minimum  to  maximum  is  assumed  to  be  available  instantly. 

■  No  transients  of  ANY  KIND  are  included. 

■  Fastest  model  execution.  Adequate  in  many  vehicle  models,  including  series  hybrid 
applications  where  the  engine  is  operating  at  steady  or  nearly  steady  conditions. 
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The  following  assumptions  have  been  used  in  the  Spark  Ignited  (SI)  engine  transient  model: 


■  Liquid  fueled  (e.g.  gasoline).  Should  not  be  used  for  gas  fueled  (e.g.  natural  gas)  engines. 

■  Naturally  aspirated  (no  turbocharging  or  supercharging).  Most  SI  engines  are  naturally 
aspirated. 

■  Mean  value  torque  model.  Estimates  average  torque  over  engine  cycle. 

■  Does  not  compute  torque  variation  within  a  cycle. 

■  Assumes  continuous  flows.  Does  not  consider  manifold  pressure  waves  or  tuning  Tuning 
may  be  included  through  effect  on  volumetric  efficiency. 

■  Three  internal  mode  states  associated  with  intake  manifold  emptying/filling  manifold  wall 
wetting  by  fuel  catalyst  bed  temperature. 

■  Model  does  not  consider  any  time  delays  associated  with: 

fueling  decisions 
fuel  transport  lag 
torque  production 

■  Engine  is  assumed  to  be  fully  warmed  up  at  all  times.  Cold-start  emissions  are  not 
calculated. 

■  Model  does  not  simulate  cranking  or  associated  emissions. 
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3.0  ELECTRIC  MACHINE  LIBRARY 


The  electric  machine  library  contains  models  of  an  AC  induction  motor  as  well  as  a  permanent 
magnet  generator.  The  AC  induction  motor  model  is  connected  to  the  wheel  of  an  HEV  either 
directly  or  via  a  motor  wheel  coupling.  The  generator  is  connected  to  an  engine  via  an  engine- 
generator  coupling. 


□ 


Permanent  Magnet 
Machines 


□ 


Induction 

Machines 


Note:  The  modules  with  a  "D"  in  the 
upper  corners  are  for  dynamic  simula¬ 
tions  with  a  maximum  time  step  of  1  ms. 


Figure  3.  Machine  Model  Library 


3.1  Steady-State  Permanent  Magnet  Generator  System  Module  with  Fixed 
Efficiency  Rectifier 


INPUTS 

OUTPUTS 

DC  Bus  Voltage  [V] 

Rectifier  Current  [A] 

DC  Command  Current  [A] 

[Output  Torque  [Nm],  Machine  Inertia  [kgm2]] 

Engine  Speed  [rpm] 

Constants  as  defined  by  the  mask: 

VAr:  Rated  Power  of  Generator  [W] 

Vr:  Rated  Line-Neutral  Voltage  of  Generator  [V] 

nr:  Rated  Speed  of  Generator  [rad/s] 

poles:  Number  of  Generator  Poles  [  ] 

lb:  Rated  Line  Current  of  Generator  [A] 

Zb:  Base  Impedance  of  Generator  for  Per  Unit  [Q] 

wb:  Base  Speed  of  Generator  for  Per  Unit  [rad/s] 
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Rspu: 

Per  Unit  Stator  Resistance  [  ] 

Rmpu: 

Per  Unit  Iron  Loss  Resistance  [  ] 

Lqspu: 

Per  Unit  Stator  q-axis  Inductance  [  ] 

Ldspu: 

Per  Unit  Stator  d-axis  Inductance  [  ] 

Llspu: 

Per  Unit  Stator  Leakage  Inductance  [  ] 

Lmdpu: 

Per  Unit  Stator  d-axis  Magnetizing  Inductance  [  ] 

Idfpu: 

Per  Unit  PM  “Field”  Current  [  ] 

recteff: 

Rectifier  Efficiency  [  ] 

Inertia_out: 

Machine  Inertia  [kg-m2] 

sample_N: 

Number  of  Samples  to  Save  to  the  MATLAB®  Workspace 

sample_T: 

Length  of  Time  Between  Each  Sample 

Auxiliary  Calculations:  None 


Mask: 

1 .  SS  PM  Generator  System 

2.  SS  PM  Generator/Rectifier  System|Gen 
Rtgs[P(hp),V(V),n(rpm)]:|Rtgs[poles,NLV(V),Core  L(hp)]:|[Rs  (ohms), Lqs,Lds, Lis 
(H)]:|Rectifier  Efficiency:|Machine  Inertia(kg*m*m):|  [sample_N,  sample_T]: 

3.  mvr=gencinit(@l,@2,@3);  VAr  =  mvr(l,l);Vr  =  mvr(l,2);  nr  =  @1(1,3); 

poles  =  @2(1,1);  Rmpu  =  mvr(l,3);  lb  =  mvr(l,4);  Zb  =  mvr(l,5);  wb  =  mvr(l,6);  Rspu  = 
mvr(l,7);  Lqspu=mvr(l,8);  Ldspu  =  mvr(l,9);  Llspu  =  mvr(l,10);  Lmdpu  =  mvr(l,l  1); 
Idfpu=mvr(l,12);  recteff  =  @4;  Inertia_out  =  @5;  sample_N  =  @6(1,1); 
sample_T=@6(  1 ,2) ; 


This  module  simulates  a  permanent-magnet  generator  connected  to  a  dc  bus  through  a  three- 
phase  rectifier.  The  module  mask  contains  the  example  values  for  a  50  hp,  230  V  generator, 
which  is  user  changeable.  The  generator  is  connected  to  the  dc  bus  through  a  three-phase, 
fully  controlled  PWM  voltage  source  (boost-type)  rectifier  as  shown  in  Figure  4.  Please  refer 
to  the  following  assumptions  for  other  important  simulation  details. 


Assumptions: 

■  This  module  will  only  be  used  in  simulations  with  a  minimum  time  step  of  100  ms  or 
larger. 
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Figure  4.  -  Generator  System  Configuration.  Ideal  Switches  Can  Be  IGBTS,  BJTS  or 

GTOS  With  Anti-Parallel  Diodes 

■  The  generator  sign  convention  is  that  positive  current  is  leaving  the  terminals  of  the 
machine.  This  current  produces  positive  torque  for  positive  speeds.  For  the  module 
itself,  a  positive  actual  dc  rectifier  current  (output  #1)  indicates  the  generating  mode. 

■  Model  does  not  include  rectifier  harmonic  effects. 

■  The  rectifier  is  assumed  to  have  a  fixed,  user  definable  efficiency.  The  losses  in  the 
generator  are  modeled  as  stator  copper  and  core  loss  (R.  and  R^,,  respectively  in  Figure  6). 

■  Within  the  module,  a  set  of  power  values  is  computed  at  each  time  step  in  order  to 
facilitate  energy  balance  equations.  The  power  values  are  as  follows:  the  dc  bus  power, 
the  mechanical  power  of  the  generator,  the  electrical  power  at  the  terminals  of  the 
machine,  the  rectifier  loss  power,  and  the  total  motor  losses  (core  &  copper  losses). 
These  values  are  not  saved  to  the  workspace  automatically. 

■  The  generator  is  controlled  at  all  times  with  conventional  “brushless  dc”  (also  known  as 
“field-oriented”  or  “vector”)  control.  That  is,  the  current  in  the  generator  is  controlled  to 
be  in  phase  with  the  internal  generator  voltage  (back  emf). 

■  No  provision  is  made  for  field  weakening  of  the  PM  machine.  If  sufficient  bus  voltage  is 
unavailable  to  maintain  the  control  (this  typically  occurs  when  the  speed  is  above  the 
rated  speed  with  rated  bus  voltage),  then  the  generator  is  turned  off  (output  power  to  the 
bus  is  zero)  due  to  an  “overspeed  condition”  for  the  given  bus  voltage. 
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■  The  generator  current  is  limited  to  1.5  times  the  value,  which  results  in  rated  torque  at 
rated  speed.  This  also  limits  the  direct  axis  flux  in  the  machine. 

■  The  user  is  cautioned  that  the  machine  parameters  and  rated  values  are  not  necessarily 
independent  of  one  another.  It  is  certainly  possible  to  input  values  in  the  mask  which 
results  in  an  uncontrollable  machine  or  in  control  over  a  very  small  operating  range. 

Figure  5  shows  the  inputs  and  outputs  of  the  steady-state  generator  model.  Figure  6  depicts 
the  steady-state  synchronous  d-q  reference  frame  model  used  in  the  simulation.  The 
electrical  terminal  voltage  is  calculated  directly  from  this  equivalent  circuit  using  the 
command  values  of  current  and  speed.  The  electrical  and  mechanical  equations  used  in  the 
simulation,  and  given  by  the  equivalent  circuit,  are, 


iqs  iqm  n  r^ds^dm  ) 

m 

•e  _  je  _  ®rLqsiqm 
* ds  l dm 

m 

Vqs  =  ~RJqs  +  ^r^md^df  ~  ^r^dJdm 

V*  =  ~RJds  +  COrLqfqm 

^  ~  dfiqm 


where  the  superscript  e  denotes  the  synchronous  reference  frame,  and  Idf  is  the  fictitious 
current  that  corresponds  to  the  strength  of  the  permanent  magnet.  This  current  is  internally 
calculated  from  the  no  load  speed  and  voltage.  The  iron  loss  resistor,  R„, ,  is  calculated  from 
the  rated  core  losses  value  provided  by  the  user. 
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Figure  5  -  Inputs  and  Outputs  of  PM  Generator  System 


Figure  6  -  Steady-State  dq-Equivalent  Circuits 

The  generator  is  controlled  using  conventional  synchronous  reference  frame  field  orientation 
(i.e.,  “brushless  dc”  control)  such  that  i^  =  0. 

When  the  mask  is  activated,  the  window  in  Figure  7  is  displayed.  The  first  input  line  asks  the 
user  to  list  the  rated  power,  rated  line,  voltage,  and  rated  speed  of  the  generator.  The  next 
input  line  requests  the  values  for  the  number  of  poles,  the  no-load  voltage  at  rated  speed,  and 
the  core  losses  at  rated  speed.  The  machine  parameters  of  the  per-phase  stator  resistance, 
per-phase  q-  and  d-axis  stator  self  inductances,  and  the  per-phase  leakage  inductance  are 
defined  in  the  third  line  of  the  mask.  The  fixed  efficiency  of  the  rectifier  is  specified  in  line 
4,  and  the  system  inertia  is  defined  in  the  fifth  input  vector.  Finally,  the  user  can  specify  the 
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number  of  samples  and  the  time  step  between  each  sample  for  use  in  “To  Workspace”  blocks 
within  the  module.  The  default  values  listed  in  Table  3  describe  a  50  horsepower  permanent- 
magnet  synchronous  generator  with  a  fixed  efficiency  rectifier. 


Figure  7  -  Window  Displayed  When  Mask  Is  Activated  For  The  Permanent  Magnet 
Generator  Module  With  Fixed  Rectifier  Efficienc 


Table  3.  Machine  Definitions  and  Values  I 

Machine  Definitions 

Values 

Rated  Power 

50  hp 

Rated  Voltage 

230  V 

Rated  Speed 

1800  rpm 

Number  of  Poles 

4 

No-Load  Voltage  @  Rated  Speed 

93  V 

Core  Losses  @  Rated  Voltage  &  Speed 

0.25  hp 

Per-Phase  Stator  Resistance,  Rs 

1.419E-2  Q 

Per-Phase  Stator  q-axis  Inductance,  L0 

1.882E-3H 

Per-Phase  Stator  d-axis  Inductance,  Ld 

1.505E-3  H 

3.763E-4  H 

Rectifier  Efficiency 

0.98 

Machine  Inertia 

0.1  kgm2 

Number  of  Samples  Saved  to  Workspace 

20000 

Time  Between  each  Saved  Sample 

0.1 
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3.2  Steady-State  Permanent  Magnet  Generator  System  Module  with  Calculated 
Rectifier  Losses 


INPUTS 

OUTPUTS 

DC  Bus  Voltage  [V] 

Rectifier  Current  [A] 

DC  Command  Current  [A] 

Output  Torque  [N-mj 

Engine  Speed  [rpm] 

Machine  Inertia  [kg  m2] 

Constants  as  defined  by  the  mask: 


VAr: 

Vr: 

nr: 

poles: 

lb: 

Zb: 

wb: 

Rspu: 

Rmpu: 

Lqspu: 

Ldspu: 

Llspu: 

Lmdpu: 

Idfpu: 

fc: 

EQD: 

Vnom: 

rQ: 

rD: 

VQ: 

VD: 

PWM: 

Inertia_out: 

sample_N: 

sampleT: 


Rated  Power  of  Generator  [W] 

Rated  Line-Neutral  Voltage  of  Generator  [V] 

Rated  Speed  of  Generator  [rad/s] 

Number  of  Generator  Poles  [  ] 

Rated  Line  Current  of  Generator  [A] 

Base  Impedance  of  Generator  for  Per  Unit  [Q] 

Base  Speed  of  Generator  for  Per  Unit  [rad/s] 

Per  Unit  Stator  Resistance  [  ] 

Per  Unit  Iron  Loss  Resistance  [  ] 

Per  Unit  Stator  q-axis  Inductance  [  ] 

Per  Unit  Stator  d-axis  Inductance  [  ] 

Per  Unit  Stator  Leakage  Inductance  [  ] 

Per  Unit  Stator  d-axis  Magnetizing  Inductance  [  ] 

Per  Unit  PM  “Field”  Current  [  ] 

PWM  Switching  Frequency  of  the  Rectifier  [Hz] 

Total  Switching  Energy  Loss  per  Amp  of  the  IGBT-Diode  Combination  [J/A] 
Voltage  at  which  EQD  is  specified  [V] 

On  State  IGBT  Resistance  [Q] 

On  State  Diode  Resistance  [Q] 

Forward  Voltage  Drop  of  the  IGBT  [V] 

Forward  Voltage  Drop  of  the  Diode  [V] 

Determines  if  Sinusoidal(l)  or  Space  Vector  (0)  PWM  is  used.[  ] 
Machine  Inertia  [kg.m2] 

Number  of  Samples  to  Save  to  the  MATLAB  Workspace  [  ] 

Length  of  Time  Between  Each  Sample  [  ] 


Auxiliary  Calculations:  None 


Mask: 

1 .  SS  PM  Generator  System 

2.  SS  PM  Generator  /  Rectifier  System|Gen  Ratings[P(hp),V(V),n(rpm),J(kgm2)]| 

Rtgs [poles, NLV(V), Core  L(hp)]|[Rs(ohms),Lqs(H),Lds(H),Lls(H)]|  Rect 
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[rQ,rD(ohms),VQ,VD(V)]|[PWM,f(Hz),EQD(J/A),Vnom(V)]|  [sample_N, 
sample_T]: 

3.  mvr=gencinit(@  1  ,@2,@3);VAr=mvr(  1 , 1 ); Vr=mvr(l ,2);nr=@  1  ( 1 ,3);poles-@2(  1,1); 
Rmpu=mvr(l,3);Ib=mvr(l,4);Zb=mvr(l,5);wb=mvr(l,6);Rspu=mvr(l,7); 
Lqspu=mvr(l  ,8);Ldspu=mvr(  1 ,9);Llspu=mvr(  1 , 1 0);Lmdpu=mvr(l ,  11); 
IdQ)u=mvr(l,12);rQ=@4(l,l)3-D=@4(l,2);VQ=@4(l,3);VD=@4(l,4);PWM=@5(l,l); 
fc=@5(  1 ,2);Inertia_out=@  1(1 ,4);EQD=@5(1 ,3);Vnom=@5(l  ,4);sample_N=@6(l  ,1); 
sample_T=@6(l,2); 

This  module  simulates  a  permanent-magnet  generator  connected  to  a  dc  bus  through  a  three- 
phase  rectifier.  The  module  mask  contains  example  values  for  a  50  hp,  230V  generator, 
which  are  user  changeable.  The  generator  is  connected  to  the  dc  bus  through  a  three-phase 
fully  controlled  PWM  voltage  source  (boost-type)  rectifier  as  shown  in  Figure  8.  The  rectifier 
uses  either  the  sinusoidal  PWM  method  [1]  or  the  space  vector  PWM  method  [2], 


Assumptions: 

■  This  module  will  only  be  used  in  simulations  with  a  minimum  time  step  of  100  ms  or 
larger. 

■  The  generator  sign  convention  is  such  that  positive  current  leaves  the  terminals  of  the 
machine.  This  current  produces  positive  torque  for  positive  speeds.  A  positive  dc 
rectifier  current  (output  #1)  indicates  the  generating  mode. 


Figure  8  -  Generator  System  Configuration: 

Ideal  Switches  Can  Be  IGBTs,  BJTs  or  GTOs  With  Antiparallel  Diodes 


25 


■  Within  the  module,  a  set  of  power  values  is  computed  at  each  time  step  in  order  to 
facilitate  energy  balance  equations.  The  power  values  are  as  follows:  the  dc  bus  power, 
the  mechanical  power  of  the  generator,  the  electrical  power  at  the  terminals  of  the 
machine,  the  rectifier  loss  power,  and  the  total  motor  losses  (core  &  copper  losses). 
These  values  are  not  saved  to  the  workspace  automatically. 

■  Loss  calculations  for  the  Alt-Rev  space  vector  PWM  method  [2]  are  based  on  the  third 
harmonic  PWM  method  since  both  control  methods  result  in  similar  steady-state 
operation. 

■  Phase  currents  and  phase  voltages  of  the  generator  are  assumed  to  be  purely  sinusoidal 
for  rectifier  loss  calculations.  This  assumption  is  justifiable  if  the  PWM  frequency  is 
reasonably  high. 

■  The  machine  model  does  not  include  rectifier  harmonic  effects. 

■  The  losses  in  the  generator  are  modeled  as  stator  copper  loss  (R.,  in  Figure  10)  and  core 
loss  (R™  in  Figure  10). 

■  The  generator  is  controlled  at  all  times  with  conventional  “brushless  dc”  (also  known  as 
“field-oriented”  or  “vector”)  control.  That  is,  the  current  in  the  generator  is  controlled  to 
be  in  phase  with  the  internal  generator  voltage  (back  emf). 

■  No  provision  is  made  for  field  weakening  of  the  PM  machine.  If  sufficient  bus  voltage  is 
unavailable  to  maintain  the  control  (this  typically  occurs  when  the  speed  is  above  the 
rated  speed  with  rated  bus  voltage),  the  generator  is  turned  off  (output  power  to  the  bus  is 
zero)  due  to  an  “overspeed  condition”  for  the  given  bus  voltage. 

■  The  generator  current  is  limited  to  1.5  times  the  value,  which  results  in  rated  torque  at 
rated  speed.  This  also  limits  the  direct  axis  flux  in  the  machine. 

■  The  user  is  cautioned  that  the  machine  parameters  and  rated  values  are  not  necessarily 
independent  of  one  another.  It  is  certainly  possible  to  input  values  in  the  mask,  which 
result  in  an  uncontrollable  machine,  or  in  control  over  a  very  small  operating  range. 
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Figure  10  shows  the  inputs  and  outputs  of  the  steady-state  permanent  magnet  generator 
model  with  calculated  rectifier  losses.  Figure  11  depicts  the  steady-state  synchronous  d-q 
reference  frame  machine  model  used  in  the  simulation.  The  electrical  terminal  voltage  is 
calculated  directly  from  this  equivalent  circuit  using  the  command  values  of  current  and 
speed.  The  electrical  and  mechanical  equations  used  in  the  simulation,  and  given  by  the 
equivalent  circuit,  are, 


V*  V>  D  ^r^dshm  ^r^md^df) 


1 

R, 

n-hA 

K 

V?i  =  ~Rjq*  +  ^r^md^df  ~  ^r^ds^dm 

-RA  +  coLJL 


ie  —  ie  —  r  qs  Qm 
lds  1  dm 


^  Rntd^dfigm 


where  the  superscript  e  denotes  the  synchronous  reference  frame  and  Idf  is  the  fictitious 
current  which  corresponds  to  the  strength  of  the  permanent  magnet.  This  current  is  internally 
calculated  from  the  no  load  speed  and  voltage.  The  iron  loss  resistor,  R^,,  is  calculated  from 
the  rated  core  losses  value  provided  by  the  user. 
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Figure  9.  Steady-State  Permanent  Magnet  Generator  with  Calculated  Rectifier  Losses 
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Figure  10  -  Inputs  and  Outputs  of  PM  Generator  System 
with  Calculated  Rectifier  Losses. 


The  generator  is  controlled  using  conventional  synchronous  reference  frame  field  orientation 
(i.e.,  “brushless  dc”  control)  such  that  =  0. 


The  following  are  the  machine  variables  used  for  the  device  loss  calculations  within  the 
rectifier: 


where 


VP 

b 

M 


IP  =y/2Ib  ■yj2  {iqspu2  +  Idspu2 ) 
Vp=JlVr  ^(Vqspu2  +Vdspu2) 


(j>  =  tan' 


Vdspu 

Vqspu 


■tan' 


Idspu 

Iqspu 


M  = 


El 

Vdc 


Peak  machine  output  phase  voltage 
Peak  output  current  of  the  machine 
Modulation  index 
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* 

Phase  lag  of  the  machine  line  current 

Vdc 

Voltage  of  the  DC  bus. 

Iqspu 

RMS  q-axis  machine  output  current 

Idspu 

RMS  d-axis  machine  output  current 

Vqspu 

RMS  q-axis  machine  output  phase  voltage 

Vdspu 

RMS  d-axis  machine  output  phase  voltage 

Table  4  shows  the  equations  governing  the  rectifier  device  currents. 


Table  4.  Equations  Governing  Rectifier  Device  Currents 

Device  Currents 

|  Sinusoidal  PWM 

Space  Vector  PWM 

h 

'l 

1  (  Mcos^ 

2n+  8  J 

/ 

P 

1  |  a/JMcos^ 
^2;r  +  12  , 

h 

i 

1  M  cos^ 

2  n  8  J 

■ 

\lQ\ 

\  V\ RMS 

i 

p 

/l  (  Mcos</> 

V8  3  n 

j  1 1  (  ^30cos^-cos3^ 

H8  45Sn 

K\ 

1 1  M  cos  <f> 

18  3n 

I  P  30COS^-COS3^ 

H8  45y/3n 

where 


IQ  Average  current  through  the  IGBT 

I o  Average  current  through  the  diode 

| lQ  |  RMS  current  through  the  IGBT 

\Id\rms  RMS  current  through  the  diode 


The  rectifier  losses  are  then  calculated  using  the  following  equations: 
Conduction  losses: 
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PQ(cond)  -Iq-VQ  +  l^el^  •  rQ 
PD(cond)  =  I d'VD  +  \ID\rms  ■  rD 

where  Pq^  Conduction  power  loss  in  the  IGBT 

P  Conduction  power  loss  in  the  diode 

D(cond) 


Switching  losses: 


p  = 
rQ(sw) 


f  '  En  ■!„(  Vdc  \ 


*[ 


Vnom 


D(sw) 


_f-^D-Jp(Vdc) 


n 


PtOT(sw)  ~  Pq(sw)  +  Pd(sw) 


Vnom  J 

f  -EQD-Ip 


7T 


Vdc  N 
Vnom j 


where 


Pq(sw)  Switching  power  loss  in  the  IGBT 

PD(sw)  Switching  power  loss  in  the  diode 

Ptot(sw)  Total  switching  power  loss 

Vdc  Voltage  of  the  DC  bus. 

EQ  Switching  energy  loss  per  Amp  of  the  IGBT 

Ed  Switching  energy  loss  per  Amp  of  the  diode 

EQD  Eg  +  Ed 


The  total  rectifier  loss  power(P/0^)  is  six  times  the  sum  of  the  conduction  and  switching 
losses  of  the  IGBT  and  the  diode. 


[  loss 


-  6(p, 


+  Pn 


Q(sw)  1  *  D(sw) 


P()(cond)  Pofcond)  j 


IGBT  Selection  Guide: 

Table  5  shows  the  IGBT  modules  that  contain  the  antiparallel  diodes. 
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1  Table  5.  IGBT  Modules  tl 

lat  Contain  the  Antiparallel  Diodes  I 

IGBT  Module# 

Frequency 

Vces 

Ic© 

Pdmax 

VQmax  @ 

VDmax© 

EQD(@V™) 

rQ  @ 

rD@ 

Range 

(V) 

Tj=25°C 

(W) 

T  j=25°C 

T  j=25°C 

@Tj=125°C 

Tj=125°C 

Tj=125°C 

(kHz) 

(V) 

(V) 

(mJ/A) 

(mQ) 

(mQ) 

IRGDDN300M06 

3-10 

600 

400 

1563 

2.0 

2.0 

0.2@Vmn=300V 

— 

_ 

IRGDDN400M06 

3-10 

600 

600 

1984 

2.0 

2.0 

0.2  @  Vnom=300V 

_ 

IRGDDN600M06 

3-10 

600 

800 

2604 

2.7 

2.0 

0.2@Vnm=300V 

_ 

— 

IRGDDN200M12 

3-10 

1200 

420 

1800 

2.7 

3.4 

0.6@Vim=600V 

2.00 

5.50 

IRGDDN300M12 

3-10 

1200 

560 

2400 

2.7 

3.4 

o.eov^oov 

2.50 

4.00 

IRGDDN400M12 

3-10 

1200 

400 

2770 

2.7 

3.4 

0.6@Vnom=600V 

_ 

IRGDDN300K06 

10-30 

600 

340 

1563 

2.7 

2.0 

0.15@Vnom=300V 

— 

_ 

IRGDDN400K06 

10-30 

600 

520 

1984 

2.7 

2.0 

o.isfgv^oov 

_ 

_ 

iRGDDN600K06 

10-30 

600 

680 

2604 

2.7 

2.0 

0.15@V_=300V 

- 

- 

where 


Vces  Maximum  Continuous  Collector  to  Emitter  Voltage 

Ic  Maximum  Continuous  Collector  Current 

Pdmax  Maximum  Power  dissipation 

VQmax  Maximum  Collector  to  Emitter  ON  Voltage 

VD^  Maximum  Forward  Voltage  Drop  of  the  Diode 

EQD  Total  Switching  Energy  Loss  per  Amp  of  the  IGBT-Diode  Combination 

V„nn  Collector  Voltage  at  which  EQD  is  specified 

rQ  On  State  IGBT  Resistance 

rD  On  State  Diode  Resistance 

Tj  Device  Junction  Temperature 


The  value  of  on-state  resistance  of  the  IGBT  and  diode  is  found  by  calculating  the  inverse  of 
the  slope  of  the  device  current  versus  device  voltage  curve  (during  conduction). 


When  the  mask  is  activated,  the  window  in  Figure  1 1  is  displayed.  The  first  input  line  asks 
the  user  to  list  the  rated  power,  rated  line  voltage,  and  rated  speed  of  the  generator  as  well  as 
the  system  inertia.  The  next  input  line  requests  the  values  for  the  number  of  poles,  the  no 
load  voltage  at  rated  speed,  and  the  core  losses  at  rated  speed.  The  machine  parameters  of 
per-phase  stator  resistance,  per-phase  q-  and  d-axis  stator  self  inductances,  and  the  per-phase 
leakage  inductance  are  defined  in  the  third  line  of  the  mask.  The  rectifier  is  described  by  the 
next  two  input  vectors  of  the  mask.  In  the  first  vector,  the  on-state  resistance  of  the  switching 
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devices(rQ)  and  the  diodes(rD)  and  their  respective  forward  voltage  drops  are  defined.  The 
second  input  vector  for  the  rectifier  specifies  the  PWM  method  (1  for  Sinusoidal/O  for  Space 
Vector),  the  switching  frequency,  the  total  switching  energy  loss  per  amp  of  the  switching 
device/diode  combination,  and  the  collector  voltage  at  which  this  switching  loss  occurs. 
Finally,  the  user  can  specify  the  number  of  samples  and  the  time  step  between  each  sample 
for  use  in  “To  Workspace”  blocks  within  the  module.  The  default  values,  which  describe  a 
50  horsepower  permanent-magnet  synchronous  generator  with  a  three-phase  IGBT  boost 
rectifier,  are  listed  in  Table  6  on  the  following  page. 


Figure  12.  Window  Displayed  When  Mask  Is  Activated  For  The  Permanent  Magnet 
Generator  Module  With  Calculated  Rectifier  Losses 
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Table  6.  Default  Values  I 

Machine  and  Rectifier  Definitions 

Values 

Rated  Power,  P 

Rated  Voltage,  V 

230  V 

Rated  Speed,  n 

Number  of  Poles,  poles 

No-Load  Voltage  @  Rated  Speed,  NLV 

f  93V 

Core  Losses  @  Rated  Voltage  &  Speed,  Core  L 

0.25  hp 

Per-Phase  Stator  Resistance,  Rs 

1.419E-2Q 

Per-Phase  Stator  q-axis  Inductance,  Las 

1.882E-3H 

Per-Phase  Stator  d-axis  Inductance,  Lds 

1.505E-3H 

Per-Phase  Stator  Leakage  Inductance,  Lb 

3.763E-4  H 

Machine  Inertia,  J 

PWM  Switching  Frequency  of  the  Rectifier,  fc 

10000  Hz 

Total  Switching  Energy  Loss  per  Amp  of  the  IGBT-Diode  Combination,  EQD 

160  DJ/A 

Voltage  at  which  EQD  is  specified,  Vnom 

350  V 

On  State  IGBT  Resistance,  rQ 

4.0E-3  □ 

On  State  Diode  Resistance,  rD 

4.0E-3  □ 

Forward  Voltage  Drop  of  the  IGBf ,  VQ 

1  V 

Forward  Voltage  Drop  of  the  Diode,  VD 

1  V 

type  of  PWM  used  [Sinusoidal(l)  or  Space  Vector(0)] 

1 

Number  of  Samples  to  Save  to  the  MATLAB  Workspace,  sample  N  “1 

20000 

Length  of  Time  Between  Each  Sample,  sample  T 

0.1 

References: 

[1]  Ken  Berringer,  Jeff  Marvin  and  Philippe  Perruchoud,  “Semiconductor  Power  Losses 
in  AC  Inverters,”  Conference  Record-IAS  Annual  Meeting  (IEEE  Industry  Applications 
Society),  pp.  882-886, 1995. 

[2]  Philippe  J.P.  Perruchoud  and  Peter  J.  Pinewski,  “Power  Losses  for  Space  Vector 
Modulation  Techniques,”  IEEE  Workshop  on  Power  Electronics  in  Transportation,  pp. 
167-173,  Oct.  1996. 


3.3  Steady-State  Induction  Generator  System  with  Fixed  Efficiency  Rectifier 


INPUTS 

OUTPUTS 

DC  Bus  Voltage  [V] 

Rectifier  Current  [A] 

DC  Command  Current  [A] 

[Output  Torque  [N  m],  Machine  Inertia  [kgm2]] 

Engine  Speed  [rpm] 
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Constants  as  defined  by  the  mask: 


VAr: 

Rated  Power  of  Motor  [W] 

Vr: 

Rated  Line-Neutral  Voltage  of  Motor  [V] 

wb: 

Base  Speed  of  Motor  for  Per  Unit  [rad/s] 

wmr: 

Rated  Speed  of  Motor  [rad/s] 

poles: 

Number  of  Motor  Poles  [  ] 

corel: 

Core  Losses  of  Motor  @  Rated  Voltage  &  Speed  [W] 

SLL: 

Stray  Load  Losses  of  Motor  @  Rated  Voltage  &  Speed  [W] 

lb: 

Base  Line  Current  of  Motor  [A] 

Zb: 

Base  Impedance  of  Motor  for  Per  Unit  [Q] 

Rspu: 

Per  Unit  Stator  Resistance  [  ] 

Rrpu: 

Per  Unit  Rotor  Resistance  [  ] 

Lspu: 

Per  Unit  Stator  Inductance  [  ] 

Lrpu: 

Per  Unit  Rotor  Inductance  [  ] 

Lmpu: 

Per  Unit  Magnetizing  Inductance  [  ] 

fluxr: 

Per  Unit  Rotor  Flux  @  Rated  Conditions  [  ] 

recteff: 

Inverter  Efficiency  [  ] 

Inertiaout: 

Machine  Inertia  [kgm2] 

sample_N: 

Number  of  Samples  to  Save  to  the  MATLAB®  Workspace 

sampleT: 

Length  of  Time  Between  Each  Sample 

Auxiliary  Calculations:  None 


Mask: 

1 .  SS  Indue.  Generator  Sys. 

2.  SS  Induction  Generator/Rectifier  System|Gen 

Rtgs[P(hp),  V  (V),f(Hz)] :  |Rtgs[n(rpm),CL(hp),  SLL(hp)]:|[Rs,Rr  (ohms),  Lls,Llr,Lm 
(H)]:|Rectifier  Efficiency:|Machine  Inertia(kg*m*m):|[sample_N,  sample_T]: 

3.  mvr=genlinit(@l ,@2,@3);VAr=mvr(l ,l);Vr=mvr(l ,2);wb=mvr(l ,3);wmr=mvr(l,4); 
poles=mvr(l  ,5);corel=mvr(l  ,6);SLL=mvr(l ,  1 5);Ib=mvr(l  ,7);Zb:=mvr(l  ,8);Rspu=mvr(l  ,9 
);  Rrpu=mvr(  1 , 1 0);Lmpu=mvr(  1 , 1 1  );Lspu=mvr(  1 , 1 2);Lrpu=mvr(  1,13  );fluxr=mvr(  1,14); 
recteff=@4;  Inertia_out=@5;sample_N=@6(l,l);sample_T=@6(l,2);') 


This  module  simulates  a  three-phase  squirrel  cage  induction  generator  and  a  three  phase 
boost-type  or  voltage-sourced  (IGBT)  rectifier  supply.  The  module  mask  contains  the 
example  values  for  a  50  hp,  230  V  generator,  which  are  user  changeable.  The  core  losses  and 
stray  load  losses  are  included  in  the  model  and  can  be  set  in  the  module  mask  as  well.  Please 
refer  to  the  assumptions  below  for  other  important  simulation  details. 
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Assumptions: 


■  This  module  will  only  be  used  in  simulations  with  a  minimum  time  step  of  100  ms  or 
larger. 

■  The  generator  is  controlled  at  all  times  with  conventional  “field-oriented”  or  “vector” 
control.  That  is,  the  slip  frequency  and  current  in  generator  are  regulated  to  provide 
dynamic  control  of  the  generator  torque. 

■  The  flux  in  the  generator  is  limited  to  its  rated  value  (~  rated  V/f). 

■  The  torque  is  limited  to  the  peak  or  pullout  torque  for  rated  flux.  The  current  is  therefore 
limited  to  the  current  at  peak  torque. 

■  When  the  bus  voltage  is  not  sufficient  to  provide  rated  flux,  field  weakening  is  enabled 
which  maintains  maximum  bus  voltage  on  the  generator.  In  this  mode,  the  flux  is 
inversely  proportional  to  the  generator  speed,  and  the  peak  torque  is  inversely 
proportional  to  the  square  of  the  generator  speed. 

■  At  low  speeds,  the  rotor  copper  losses  significantly  affect  the  performance  of  the 
induction  generator.  When  in  the  motoring  mode,  the  actual  DC  current  can  overshoot 
the  requested  value  substantially.  To  protect  the  HEV  system  and  eliminate  large 
overshoot  in  the  motoring  mode,  the  q-axis  stator  current  is  limited  to  1.5  pu  in  addition 
to  the  previously  discussed  limits. 

■  The  copper  losses  are  modeled  by  R^.  and  R/  in 

■  Figure  13.  The  core  losses  are  calculated  to  be  proportional  to  the  generator  flux.  The 
stray  load  losses  are  proportional  to  the  generator  torque. 

■  The  generator  sign  convention  is  that  positive  current  is  leaving  the  terminals  of  the 
machine.  This  current  produces  positive  torque  for  positive  speeds.  For  the  module 
itself,  a  positive  actual  dc  rectifier  current  (output  #1)  indicates  the  generating  mode. 

■  Within  the  module,  a  set  of  power  values  is  computed  at  each  time  step  in  order  to 
facilitate  energy  balance  equations.  The  power  values  are  as  follows:  the  dc  bus  power, 
the  mechanical  power  of  the  generator,  the  electrical  power  at  the  terminals  of  the 
machine,  the  rectifier  loss  power,  and  the  total  motor  losses  (core,  stray  load,  and  copper 
losses).  These  values  are  not  saved  to  the  workspace  automatically. 
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■  The  model  does  not  include  rectifier  harmonic  effects. 

■  The  rectifier  is  assumed  to  have  a  fixed,  user-definable  efficiency. 


Figure  13  shows  the  inputs  and  outputs  of  the  steady-state  induction  generator  model.  Figure 
14  depicts  the  per-phase  equivalent  circuit  model  of  the  induction  machine,  and  defines  the 
generator  parameters.  Figure  15  shows  the  induction  generator  model  after  the  transformation 
to  the  synchronous  d-q  reference  frame.  The  electrical  terminal  voltages  and  currents  are 
calculated  directly  from  this  equivalent  circuit  and  the  mechanical  equations  using  the 
command  values  of  torque  and  speed.  The  synchronous  reference  frame  electrical  and 
mechanical  equations  used  in  the  simulation  are, 

Ks  =-Rfqs+®Ks 

vds  =-RJds-<aAeqs 

0  =  Rrieqr+(a>-a)r)A,edr 
0  =  Rri‘r-(<o-a>r)Aegr 

K*  =  + 

A'ds  ~  +  Lls)ids  +  Lmidr 

%r=(Lm+LlrK-LJeqs 

^ dr  =(An  +  Rlr  )hr  ~  Rm^ds 

r  =  Lm{iedsieqr-ieqfdr) 

where  the  e  superscript  denotes  the  synchronous  reference  frame. 


DC  Current 
Command  [A] 


Figure  13.  Inputs  And  Outputs  Of  The  Steady-State  Induction  Generator  Model  With 

A  Fixed  Efficiency  Rectifier 
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The  generator  is  controlled  using  conventional  rotor  flux  field  orientation  (“brushless  dc” 
control)  such  that  in  the  synchronous  reference  frame,  Xqre=  0. 

When  the  mask  is  activated,  the  window  in  Figure  15  is  displayed.  The  first  input  line  asks 
the  user  to  list  the  rated  power,  rated  line  voltage  and  rated  frequency  of  the  induction 
generator.  The  next  input  line  requests  the  values  for  the  rated  speed,  and  the  rated  core  and 
stray  load  losses  that  occur  under  rated  operating  conditions.  The  per-phase  machine 
parameters  of  stator  and  rotor  resistance,  stator  and  rotor  leakage  inductance,  and  the 
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magnetizing  inductance  are  defined  in  the  third  line  of  the  mask.  The  fixed  efficiency  of  the 
rectifier  is  specified  in  line  4,  and  the  system  inertia  is  defined  in  the  fifth  input  vector. 
Finally,  the  user  can  specify  the  number  of  samples  and  the  time  step  between  each  sample 
for  use  in  “To  Workspace”  blocks  within  the  module.  The  default  values  listed  in  Table  7 
below  describe  a  50  horsepower  induction  generator  with  a  fixed  efficiency  rectifier. 


Figure  16.  Window  Displayed  When  Mask  Is  Activated  For  The  Induction  Generator 


Module  With  Fixed  Rectifier  Efficienc 


Table  7.  Default  Values 

Machine  Definitions 

Values 

Rated  Power 

50  hp 

Rated  Voltage 

230  V 

Rated  Frequency 

60  Hz 

Rated  Speed 

1740  rpm 

Core  Losses  @  Rated  Voltage  &  Speed 

0.25  hp 

Stray  Load  Losses  @  Rated  Voltage  &  Speed 

0.25  hp 

Per-Phase  Stator  Resistance,  Rs 

2.8328E-2  Q 

Per-Phase  Rotor  Resistance,  Rr 

3.547E-2  Q 

Per-Phase  Stator  Leakage  Inductance,  L,s 

3.763E-4  H 

Per-Phase  Rotor  Leakage  Inductance,  Llr 

3.763E-4  H 

Per-Phase  Magnetizing  Inductance,  Lm 

7.527E-3  H 

Rectifier  Efficiency 

0.98 

Motor  Inertia 

0.012  [kgm2] 

Number  of  Samples  Saved  to  Workspace 

20000 

Time  Between  each  Saved  Sample 

0.1 
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3.4  Steady-State  Induction  Machine  -  Traction  Motor  System  with  Fixed 
Efficiency  Inverter 


INPUTS 

OUTPUTS 

DC  Bus  Voltage  [V] 

Inverter  Current  [A] 

Rotor  Speed  [rpm] 

[Output  Torque  [N-m],  Machine  Inertia  [kgm2] 

Per  Unit  Input  Torque  Command!  1 

Constants  as  defined  by  the  mask: 


VAr: 

Rated  Power  of  Motor  [W] 

Vr: 

Rated  Line-Neutral  Voltage  of  Motor  [V] 

wb: 

Base  Speed  of  Motor  for  Per  Unit  [rad/s] 

wmr: 

Rated  Speed  of  Motor  [rad/s] 

poles: 

Number  of  Motor  Poles  [  ] 

corel: 

Core  Losses  of  Motor  @  Rated  Voltage  &  Speed  [W] 

SLL: 

Stray  Load  Losses  of  Motor  @  Rated  Voltage  &  Speed  [W] 

lb: 

Base  Line  Current  of  Motor  [A] 

Zb: 

Base  Impedance  of  Motor  for  Per  Unit  [Q] 

Rspu: 

Per  Unit  Stator  Resistance  [  ] 

Rrpu: 

Per  Unit  Rotor  Resistance  [  ] 

Lspu: 

Per  Unit  Stator  Inductance  [  ] 

Lrpu: 

Per  Unit  Rotor  Inductance  [  ] 

Lmpu: 

Per  Unit  Magnetizing  Inductance  [  ] 

fluxr: 

Per  Unit  Rotor  Flux  @  Rated  Conditions  [  ] 

invefif: 

Inverter  Efficiency  [  ] 

Inertiaout: 

Machine  Inertia  [kg  m2] 

sample_N: 

Number  of  Samples  to  Save  to  the  MATLAB®  Workspace 

sample_T: 

Length  of  Time  Between  Each  Sample 

Auxiliary  Calculations:  None 


Mask: 

1 .  SS  Induction  Motor  Sys 

2.  SS  Induction  Motor/Inverter  System|Rated  [P(hp),V(V),f(Hz)]:|Rated 
[n(rpm),CL(hp),SLL(hp)]:|[Rs,Rr  (ohms),  Lls,Llr,Lm  (H)]:|Inverter  Efficiency:  |Machine 
Inertia(kg*m*m):|[sample_N,sample_T] : 

3.  mvr=motIinit(@l,@2,@3);  VAr=mvr(l,l);  Vr=mvr(l,2);  wb=mvr(l,3);  wmr=mvr(l,4); 
poles=mvr(l,5);  corel=mvr(l,6);  SLL=mvr(l,15);  Ib=mvr(l,7);  Zb=mvr(l,8); 
Rspu=mvr(l,9);  Rrpu=mvr(l,10);  Lspu=mvr(l,12);  Lrpu=mvr(l,13);  Lmpu=mvr(l,ll); 
fluxr=mvr(l,14);  inveff=@4;  Inertia_out=@5;  sample_N=@6(l,l);  sample_T=@6(l,2); 
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This  module  simulates  a  three-phase,  squirrel-cage  induction  motor  and  a  three-phase,  boost- 
type  or  voltage-sourced  (IGBT)  inverter  supply.  The  module  mask  contains  the  example 
values  for  a  50  hp,  230  V  motor,  which  are  user  changeable.  The  core  losses  and  stray  load 
losses  are  included  in  the  model  and  can  be  set  in  the  module  mask  as  well.  Please  refer  to 
the  assumptions  below  for  other  important  simulation  details. 

Assumptions: 

■  This  module  will  only  be  used  in  simulations  with  a  minimum  time  step  of  100  ms  or 
larger. 

■  The  motor  is  controlled  at  all  times  with  conventional  “field-oriented”  or  “vector” 
control.  That  is,  the  slip  frequency  and  current  in  motor  are  regulated  to  provide  dynamic 
control  of  the  motor  torque. 

■  The  flux  in  the  motor  is  limited  to  its  rated  value  (~  rated  V/f). 

■  The  torque  is  limited  to  the  peak  or  pullout  torque  for  rated  flux.  The  current  is  therefore 
limited  to  the  current  at  peak  torque. 

■  When  the  bus  voltage  is  not  sufficient  to  provide  rated  flux,  field  weakening  is  enabled, 
which  maintains  maximum  bus  voltage  on  the  motor.  In  this  mode,  the  flux  is  inversely 
proportional  to  the  motor  speed,  and  the  peak  torque  is  inversely  proportional  to  the 
square  of  the  motor  speed. 

■  The  copper  losses  are  modeled  by  and  R/  in  Figure  17.  The  core  losses  are  calculated 
to  be  proportional  to  the  motor  flux.  The  stray  load  losses  are  proportional  to  the  motor 
torque. 

■  The  sign  convention  for  the  motor  module  is  that  positive  current  enters  the  terminals  of 
the  machine.  Positive  current  produces  positive  torque,  i.e.,  positive  dc  inverter  current 
(output  #1)  indicates  the  module  is  motoring. 

■  Within  the  module,  a  set  of  power  values  is  computed  at  each  time  step  in  order  to 
facilitate  energy  balance  equations.  The  power  values  are  as  follows:  the  dc  bus  power, 
the  mechanical  power  of  the  motor,  the  electrical  power  at  the  terminals  of  the  machine, 
the  inverter  loss  power,  and  the  total  motor  losses  (core,  stray  load,  and  copper  losses). 
These  values  are  not  saved  to  the  workspace  automatically. 
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The  model  does  not  include  inverter  harmonic  effects. 

The  inverter  is  assumed  to  have  a  fixed,  user-definable  efficiency. 


Figure  17  shows  the  inputs  and  outputs  of  the  steady-state  induction  motor  model.  Figure  18 
depicts  the  per-phase  equivalent  circuit  model  of  the  induction  machine  and  defines  the 
motor  parameters.  Figure  19  shows  the  induction  motor  model  after  the  transformation  to  the 
synchronous  d-q  reference  frame.  The  electrical  terminal  voltages  and  currents  are  calculated 
directly  from  this  equivalent  circuit  and  the  mechanical  equations  using  the  command  values 
of  torque  and  speed.  The  synchronous  reference  frame  electrical  and  mechanical  equations 
used  in  the  simulation  are, 

VV  =  R/gs  +  <0%ds 

vL  =  RJeds-<»xqs 
0  =  ^4+(®"®r)4 
0  =  Rfdr-(co-a>r)Aeqr 

4  =(4  +  4)4 +44 
^=(4  +  4)4  +  44 
^4  =(4  +  4)4 +  44 
A4  =(4+4)4 +  44 
T  =  4(44  -t'dfgr) 

where  the  e  superscript  denotes  the  synchronous  reference  frame. 


Figure  17.  Inputs  And  Outputs  Of  Induction  Motor  System 
With  Fixed  Efficiency  Inverter 
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Figure  19.  Steady-State  dq  Equivalent  Circuit  Of  An  Induction  Motor 


The  motor  is  controlled  using  conventional  rotor  flux  field  orientation  (“brushless  dc” 
control)  such  that  in  the  synchronous  reference  frame,  A.qre=  0. 

When  the  mask  is  activated,  the  window  in  Figure  19  is  displayed.  The  first  input  line  asks 
the  user  to  list  the  rated  power,  rated  line  voltage  and  rated  frequency  of  the  induction  motor. 
The  next  input  line  requests  the  values  for  the  rated  speed,  and  the  rated  core  and  stray  load 
losses  that  occur  under  rated  operating  conditions.  The  per-phase  machine  parameters  of 
stator  and  rotor  resistance,  stator  and  rotor  leakage  inductance,  and  the  magnetizing 
inductance  are  defined  in  the  third  line  of  the  mask.  The  fixed  efficiency  of  the  inverter  is 
specified  in  line  4  and  the  system  inertia  is  defined  in  the  fifth  input  vector.  Finally,  the  user 
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can  specify  the  number  of  samples  and  the  time  step  between  each  sample  for  use  in  “To 
Workspace”  blocks  within  the  module.  The  default  values,  which  describe  a  50  horsepower 
induction  machine/traction  motor  with  a  fixed  efficiency  inverter,  are  listed  in  Table  8  on  the 
following  page. 


Figure  20.  Window  Displayed  When  Mask  Is  Activated  For  The  Induction  Machine  - 
Traction  Motor  Module  With  Fixed  Rectifier  Efficiency 


Table  8.  Default  Values  1 

Machine  Definitions 

Values 

Rated  Power 

50  hp 

Rated  Voltage 

230  V 

Rated  Frequency 

60  Hz 

Rated  Speed 

1740  rpm 

Core  Losses  @  Rated  Voltage  &  Speed 

0.25  hp 

0.25  hp 

Per-Phase  Stator  Resistance,  R, 

2.8328E-2  Q 

Per-Phase  Rotor  Resistance,  Rr 

3.547E-2  Cl 

Per-Phase  Stator  Leakage  Inductance,  Lk 

3.763E-4  H 

Per-Phase  Rotor  Leakage  Inductance,  Llr 

3.763E-4  H 

Per-Phase  Magnetizing  Inductance,  Lm 

7.527E-3  H 

Inverter  Efficiency 

0.98 

Motor  Inertia 

■mni 

Number  of  Samples  Saved  to  Workspace 

20000 

Time  Between  each  Saved  Sample 

0.1 
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3.5  Steady-State  Induction  Machine  -  Traction  Motor  System  with  Calculated 
Inverter  Losses 


INPUTS 

OUTPUTS 

DC  Bus  Voltage  [V] 

Inverter  Current  [A] 

Rotor  Speed  [rpm] 

[Output  Torque  [N  m],  Machine  Inertia  [kg-m2]] 

Per  Unit  Input  Torque  Command[  ] 

Constants  as  defined  by  the  mask: 


VAr: 

Rated  Power  of  Motor  [W] 

Vr: 

Rated  Line-Neutral  Voltage  of  Motor  [V] 

wb: 

Base  Speed  of  Motor  for  Per  Unit  [rad/s] 

wmr: 

Rated  Speed  of  Motor  [rad/s] 

poles: 

Number  of  Motor  Poles  [  ] 

corel: 

Core  Losses  of  Motor  @  Rated  Voltage  &  Speed  [W] 

SLL: 

Stray  Load  Losses  of  Motor  @  Rated  Voltage  &  Speed  [W] 

lb: 

Base  Line  Current  of  Motor  [A] 

Zb: 

Base  Impedance  of  Motor  for  Per  Unit  [Q] 

Rspu: 

Per  Unit  Stator  Resistance  [  ] 

Rrpu: 

Per  Unit  Rotor  Resistance  [  ] 

Lspu: 

Per  Unit  Stator  Inductance  [  ] 

Lrpu: 

Per  Unit  Rotor  Inductance  [  ] 

Lmpu: 

Per  Unit  Magnetizing  Inductance  [  ] 

fluxr: 

Per  Unit  Rotor  Flux  @  Rated  Conditions  [  ] 

Inertiaout: 

Machine  Inertia  [kg-m2] 

fc: 

PWM  Switching  Frequency  of  the  Inverter  [Hz] 

EQD: 

Total  Switching  Energy  Loss  per  Amp  of  the  IGBT-Diode  Combination 

[J/A] 

Vnom: 

Voltage  at  which  EQD  is  specified  [V] 

rQ: 

On  State  IGBT  Resistance  [Q] 

rD: 

On  State  Diode  Resistance  [Q] 

VQ: 

Forward  Voltage  Drop  of  the  IGBT  [V] 

VD: 

Forward  Voltage  Drop  of  the  Diode  [V] 

PWM: 

Determines  if  Sinusoidal(l)  or  Space  Vector  (0)  PWM  is  used.[  ] 

sample_N: 

Number  of  Samples  to  Save  to  the  MATLAB®  Workspace 

sample_T: 

Length  of  Time  Between  Each  Sample 

Auxiliary  Calculations:  None 
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Mask: 


1 .  SS  Induction  Motor  Sys 

2.  SS  Induction  Motor/Inverter  System|Rated  [P(hp),  V(V),  f(Hz),  Machine 
Inertia(kg*m*m)]:|Rated  [n(rpm),  CL(hp),  SLL(hp)]:|[Rs,Rr  (ohms),  Lis,  Llr,  Lm 
(H)]:|Inverter  [rQ,  rD  (ohms),  VQ,  VD(V)]:|Inverter  [PWM,  fc(Hz),  EQD(J/A),  Vnom 
(V)] :  |  [sample_N,sample_T] : 

3.  mdefh=motIinit(@  1  ,@2,@3);  VAr=mdefn(l,l);  Vr=mdefn(l,2);  wb=mdefii(l,3); 
wmr=mdefn(  1 ,4);  poles=mdefn(l,5);  corel=mdefn(l,6);  SLL=mdefn(l,15); 
Ib=mdefn(l,7);  Zb=mdefn(l,8);  Rspu=mdefh(l,9);  Rrpu=mdefn(l,10); 

Lspu=mdefn(l ,  1 2);  Lrpu=mdefn(l,13);  Lmpu=mdefh(l,l  1);  fluxr=mdefh(l,14); 
rQ=@4(l,l);  rD=@4(l,2);  VQ=@4(1,3);VD=@4(1,4);  PWM=@5(1,1);  fc=@5(l,2); 
EQD=@5(1,3);  Vnom=@5(l,4);  Inertia_out=@l(l,4);  sample_N=@6(l,l); 
sample_T=@6(  1,2); 

This  module  simulates  a  three-phase,  squirrel-cage  induction  motor  and  a  three-phase,  boost- 
type  or  voltage-sourced  (IGBT)  inverter  supply.  The  three  phase  inverter  employs  either  the 
sinusoidal  PWM  method  [1]  or  the  space  vector  PWM  method  [2].  The  module  mask 
contains  the  example  values  for  a  50  hp,  230  V  motor,  which  are  user  changeable.  The  core 
losses  and  stray  load  losses  are  included  in  the  model  and  can  be  set  in  the  module  mask  as 
well.  Please  refer  to  the  assumptions  below  for  other  important  simulation  details. 


Assumptions: 

■  This  module  will  only  be  used  in  simulations  with  a  minimum  time  step  of  100  ms  or 
larger. 

■  The  motor  is  controlled  at  all  times  with  conventional  “field-oriented”  or  “vector” 
control.  That  is,  the  slip  frequency  and  current  in  motor  are  regulated  to  provide  dynamic 
control  of  the  motor  torque. 

■  The  flux  in  the  motor  is  limited  to  its  rated  value  (~  rated  V/f). 

■  The  torque  is  limited  to  the  peak  or  pullout  torque  for  rated  flux.  Therefore,  the  current  is 
therefore  limited  to  the  current  at  peak  torque. 

■  When  the  bus  voltage  is  not  sufficient  to  provide  rated  flux,  field  weakening  is  enabled 
which  maintains  maximum  bus  voltage  on  the  motor.  In  this  mode,  the  flux  is  inversely 
proportional  to  the  motor  speed,  and  the  peak  torque  is  inversely  proportional  to  the 
square  of  the  motor  speed. 
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■  The  copper  losses  are  modeled  by  R.  and  R/  in  Figure  21.  The  core  losses  are  calculated 
to  be  proportional  to  the  motor  flux.  The  stray  load  losses  are  proportional  to  the  motor 
torque. 

■  The  sign  convention  for  the  motor  module  is  that  positive  current  enters  the  terminals  of 
the  machine.  Positive  current  produces  positive  torque,  i.e.,  positive  dc  inverter  current 
(output  #1)  indicates  the  module  is  motoring. 

■  Within  the  module,  a  set  of  power  values  is  computed  at  each  time  step  to  facilitate 
energy  balance  equations.  The  power  values  are  as  follows:  the  dc  bus  power,  the 
mechanical  power  of  the  motor,  the  electrical  power  at  the  terminals  of  the  machine,  the 
inverter  loss  power,  and  the  total  motor  losses  (core,  stray  load,  and  copper  losses). 
These  values  are  not  saved  to  the  workspace  automatically. 

■  Loss  calculations  for  the  Alt-Rev  space  vector  PWM  method  [2]  are  based  on  the  third 
harmonic  PWM  method  since  both  control  methods  result  in  similar  steady-state 
operation. 

■  Phase  currents  and  phase  voltages  of  the  machine  are  assumed  to  be  purely  sinusoidal  for 
inverter  loss  calculations.  This  assumption  is  justifiable  if  the  PWM  frequency  is 
reasonably  high. 

■  The  model  does  not  include  inverter  harmonic  effects. 

Figure  21  shows  the  inputs  and  outputs  of  the  steady-state  induction  motor  model.  Figure  22 
depicts  the  per-phase  equivalent  circuit  model  of  the  induction  machine  and  defines  die 
motor  parameters.  Figure  23  shows  the  induction  motor  model  after  the  transformation  to  the 
synchronous  d-q  reference  frame.  The  electrical  terminal  voltages  and  currents  are  calculated 
directly  from  this  equivalent  circuit  and  the  mechanical  equations  using  the  command  values 
of  torque  and  speed.  The  synchronous  reference  frame  electrical  and  mechanical  equations 
used  in  the  simulation  are, 


46 


V%  =  R/,s  + 

V4  =44  -*>^s 

0  =  Rrieqr+((D-G)r)Zdr 
0  =  Rfdr-(co-a>r)Xir 

4=(4  +  4)C+4C 
4,  =(4  +  4)4 +  44 
4=  (4  +4)4+44 
4,  =(4  +  4)4 +  44 
T  =  4(44  -44) 


where  the  e  superscript  denotes  the  synchronous  reference  frame. 


Per  Unit  Torque 
Command 


Figure  21.  Inputs  And  Outputs  Of  Induction  Traction  Motor  System 
With  Calculated  Inverter  Losses 
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The  motor  is  controlled  using  conventional  rotor  flux  field  orientation  (“brushless  dc” 
control)  such  that  in  the  synchronous  reference  frame,  A,qre=  0. 


where  Vp  Peak  machine  output  phase  voltage 

Ip  Peak  output  current  of  the  machine 

M  Modulation  index 

<f>  Phase  lag  of  the  machine  line  current 

Vdc  Voltage  of  the  DC  bus. 

Iqspu  RMS  q-axis  machine  output  current 

Idspu  RMS  d-axis  machine  output  current 

Vqspu  RMS  q-axis  machine  output  phase  voltage 

Vdspu  RMS  d-axis  machine  output  phase  voltage 


Table  9  shows  the  equations  governing  the  inverter  device  currents. 
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Table  9.  Equation  Governing  the  Inverter  Device  Currents 

Device  Currents 

Sinusoidal  PWM 

Space  Vector  PWM 

h 

I 

p 

I 

p 

(  1  (  V3Mcos^ 

[2^  12  J 

1 

td 

I 

p 

(  1  M  cos^ 

8  J 

I 

p 

1 

1  Iq 

\  V  RMS 

j  1 1  Mcos<f> 

Ws  +  3  n 

I  l  +  i\f3Ocos0~cos30 

p\8  45-v/T n 

Kol  RMS 

j  1 1  M  cos^ 

*>\8  3  n 

I  ft  M  ^0  cos  cos  3^ 

p\8  45^/3x 

where 


Average  current  through  the  IGBT 

Average  current  through  the  diode 
RMS  current  through  the  IGBT 

RMS  current  through  the  diode 


The  inverter  losses  are  then  calculated  using  the  following  equations: 


Conduction  losses: 

PQ(cond)  =  Iq  •  VQ  + 1  Iq  '  rQ 
^D(cond)  ~  |4L-ri> 

where  Pg^ondj  Conduction  power  loss  in  the  IGBT 

•Py  Conduction  power  loss  in  the  diode 

Switching  losses: 
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where 


f.Ea-I,(  Vdc  ) 
Q(sw)  n  Vnom ) 

P  _f-ED-I,(  Vdc' | 
D(sw)  n  V  Vnom  J 


^TOT(sw)  ~  Pq(sw)  +  Pd(sw) 


^f-EQD-ir  Vdc 
7t  l  Vnom 


Pq(sw)  Switching  power  loss  in  the  IGBT 

PD(sw)  Switching  power  loss  in  the  diode 

PT0T(sw)  Total  switching  power  loss 

Vdc  Voltage  of  the  DC  bus 

EQ  Switching  energy  loss  per  Amp  of  the  IGBT 

Ed  Switching  energy  loss  per  Amp  of  the  diode 

EQD  Eq  +  Ed 


The  total  inverter  loss  power  ( Pfoss )  is  si*  times  the  sum  of  the  conduction  and  switching 
losses  of  the  IGBT  and  the  diode. 

Pi oss  ^  ^Pq(sw)  "t”  Pd(sw)  P(Q(Lund)  ^Dfcond)  ) 


IGBT  Selection  Guide: 

Table  10  shows  the  IGBT  modules  that  contain  antiparallel  diodes. 


Table  10.  IGBT  Modules  that  Contain  Antiparallel  Diodes 


IGBT  Module# 

Frequency 

Range 

(kHz) 

VCES 

(V) 

lc@ 

Tj=25°C 

(A) 

Pdmax 

(W) 

vqmax@ 

Tj  =25°C 

(V) 

VDMAX@ 

Tj=25°C 

(V) 

EQD  (@Vnom ) 
@Tj=125°C 

(mJ/A) 

rQ(mQ) 

@Tj 

=125°C 

rD(mQ) 

@Tj=1250 

C 

IRGDDN300M06 

3-10 

600 

400 

1563 

2.0 

2.0 

0.2  @  Vnom=300V 

“ 

mn 

IRGDDN400M06 

3-10 

600 

600 

1984 

2.0 

2.0 

0.2  @  Vnoni=300V 

- 

IRGDDN600M06 

3-10 

600 

800 

2604 

2.7 

2.0 

0.2@VnOni=300V 

IRGDDN200M12 

3-10 

1200 

420 

1800 

2.7 

3.4 

0.6@Vnom=600V 

2.00 

5.50 

IRGDDN300M12 

3-10 

1200 

560 

2400 

2.7 

3.4 

0.6@VnOm=600V 

2.50 

4.00 

IRGDDN400M12 

3-10 

1200 

400 

2770 

2.7 

3.4 

0.6@Vnom=600V 

- 

IRGDDN300K06 

10-30 

600 

340 

1563 

2.7 

2.0 

015@Vnom=300V 

- 

B9H 

IRGDDN400K06 

10-30 

600 

520 

1984 

2.7 

2.0 

0.15@Vnom=300V 

- 

IRGDDN600K06 

10-30 

600 

680 

2604 

2.7 

2.0 

0.15  @  VnOm=300V 

- 

■ 
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where 

Ic 

Pdmax 

VQmax 

VDmax 

EQD 

V. 

rQ 

rD 

T, 


Maximum  Continuous  Collector  to  Emitter  Voltage 

Maximum  Continuous  Collector  Current 

Maximum  Power  dissipation 

Maximum  Collector  to  Emitter  ON  Voltage 

Maximum  Forward  Voltage  Drop  of  the  Diode 

Total  Switching  Energy  Loss  per  Amp  of  the  IGBT-Diode  Combination 

Collector  Voltage  at  which  EQD  is  specified 

On  State  IGBT  Resistance 

On  State  Diode  Resistance 

Device  Junction  Temperature 


The  value  of  on-state  resistance  of  the  IGBT  and  diode  is  found  by  calculating  the  inverse  of 
the  slope  of  the  device  current  versus  device  voltage  curve(during  conduction). 


When  the  mask  is  activated,  the  window  in  Figure  24  is  displayed.  The  first  input  line  asks 
the  user  to  list  the  rated  power,  rated  line  voltage,  and  rated  frequency  of  the  induction  motor 
as  well  as  the  system  inertia.  The  next  input  requests  the  values  for  the  rated  speed,  and  the 
rated  core  and  stray  load  losses  that  occur  under  rated  operating  conditions.  The  per-phase 
machine  parameters  of  stator  and  rotor  resistance,  stator  and  rotor  leakage  inductance,  and  the 
magnetizing  inductance  are  defined  in  the  third  line  of  the  mask.  The  inverter  is  described  by 
the  next  two  input  vectors  of  the  mask.  In  the  first  vector,  the  on-state  resistance  of  the 
switching  devices  (rQ)  and  the  diodes  (rD)  and  their  respective  forward  voltage  drops  are 
defined.  The  second  input  vector  for  the  inverter  specifies  the  PWM  method  (1  for 
Sinusoidal/O  for  Space  Vector),  the  switching  frequency,  the  total  switching  energy  loss  per 
amp  of  the  switching  device/diode  combination,  and  the  collector  voltage  at  which  this 
switching  loss  occurs.  Finally,  the  user  can  specify  the  number  of  samples  and  the  time  step 
between  each  sample  for  use  in  “To  Workspace”  blocks  within  the  module.  The  default 
values  listed  in  Table  11  on  the  following  page  describe  a  50  horsepower  induction 
machine/traction  motor  with  a  three-phase  IGBT  boost  inverter. 
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Figure  24.  Window  Displayed  When  Mask  Is  Activated  For  The  Induction  Machine  - 
Traction  Motor  Module  With  Calculated  Inverter  Losses 


Table  1 1 .  Default  Values 


Machine  and  Inverter  Definitions 

Values 

Rated  Power 

50  hp 

Rated  Power 

50  hp 

Rated  Voltage 

230  V 

Rated  Frequency 

60  Hz 

Rated  Speed 

1740  rpm 

Stray  Load  Losses  @  Rated  Voltage  &  Speed 

0.25  hp 

Core  Losses  @  Rated  Voltage  &  Speed,  Core  L 

0.25  hp 

Per-Phase  Stator  Resistance,  Rs 

2.8328E-2  Q 

Per-Phase  Rotor  Resistance,  Rr 

3.547E-2  Q 

Per-Phase  Stator  Leakage  Inductance,  L,s 

3.763E-4  H 

Per-Phase  Rotor  Leakage  Inductance,  L,r 

3.763E-4  H 

Per-Phase  Magnetizing  Inductance,  Lm 

7.527E-3  H 

Motor  Inertia 

0.012  [kg  m2] 

PWM  Switching  Frequency  of  the  Inverter,  fc 

10000  Hz 

Total  Switching  Energy  Loss  per  Amp  of  the  IGBT-Diode  Combination,  EQD 

160  DJ/A 

Voltage  at  which  EQD  is  specified,  Vnom 

350  V 

On  State  IGBT  Resistance,  rQ 

4.0E-3  □ 

On  State  Diode  Resistance,  rD 

4.0E-3  □  ! 

Forward  Voltage  Drop  of  the  IGBT,  VQ 

1  V 

Forward  Voltage  Drop  of  the  Diode,  VD 

1  V 

Type  of  PWM  used  [Sinusoidal(l)  or  Space  Vector(O)] 

1 

Number  of  Samples  to  Save  to  the  MATLAB  Workspace,  sample_N 

20000 

Length  of  Time  Between  Each  Sample,  sample_T 

0.1 
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3.6 


Table  lookup  steady-state  generator  (No  electrical  parameters) 


ale 
ig 

PM  Gen 


INPUT 

OUTPUT 

Desired  AC  Current  (Amps) 

Shaft  Torque  (N-m) 

Shaft  Speed  (rpm) 

DC  Current  (Amps) 

DC  Bus  Voltage  (V) 

Actual  AC  Current  (Amps) 

Purpose: 


Parameters: 


llimwmot 

llimw_gen 

llim_mot 

llim_gen 


vector  of  shaft  speed  (rad/s)  for  lookup  of  limit  current  (motor  mode) 
vector  of  shaft  speed  (rad/s)  for  lookup  of  limit  current  (motor  mode) 
vector  of  limit  currents  in  motor  mode 
vector  of  limit  currents  in  generator  mode 


mtrpl_Voc 

mtrplVd 

mtrplNp 

mtrpl_n 

mtrpll 

mtrplJ 


line  to  line  voltage  (V)  for  machine  at  design  speed  but  open  circuit 
line  to  line  voltage  (V)  for  machine  at  design  conditions 
number  of  poles 

shaft  speed  (rad/s)  for  machine  at  design  conditions 
line  to  line  current  (A)  for  machine  at  design  conditions 
motor  inertia  constant  (kg*mA2) 


mtrp  1  _eff  motor  efficiency  (0- 1 ) 

mtrp  1  conv  eff  converter  efficiency  (0- 1 ) 
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Auxiliary  Calculations: 

Ke  =  Voc/  cod  <Dd  is  design  speed  in  rad/s 


/*  V  *J?> 
Kt  =  - — — — — - 


co. 


lm 


Ls  = 


I 


V2  -V2 

Vd  ync 


(I*  cod*  Np/2)2 


genpl_modulel  Mask 


(1)  genp_simp 

(2)  Simple  Scaling  PM  Generator  Definitions  (Click  on  "Help"  to  list  Inport  and  Outport 
parameters.)|Limit  Current  Lookup  Vectors|Ke  &  Kt|  [Motor  Efficiency, Converter 
Efficiency]!  [Npoles,Ls]|Inertia  Constant  (kg*mA2) 

(3)  Ilimw_mot=@l(  1  ,:);Ilimw_gen=@l  (2,  :);Ilim_mot=@  1  (3,  :);Ilim _gen=@l  (4, :); 
Ke=@2(  1  );Kt=@2(2);eff_motor=@3  ( 1  );c  1  =sqrt(3 )/ @3(2);c2=sqrt(3)*  @3  (2) ; 
npLs=@4(l  )*  @4(2)/2;  J=@5 ; 

(4)  Simple\nScaling\nPM  Gen 

(5)  Simple  Scaling  PM  Generator  Block  Port  Parameters\n\nInport  1:  DC  Bus  Voltage 
(V)Anlnport  2:  Shaft  Speed  (RPM).\nInport  3:  AC  Current  Command  (A).\n\n 
Outport  1:  DC  Current  (A).\nOutport  2:  Vector  of  [Shaft  Torque  (N*m),  Inertia 
Constant  (kg*mA2)]. 


Modeling  Equations 

Vac  =  Ke  *  tod  where  is  the  machine  efficiency 


T  =  Iac*Kt*r)m  for  motor  mode  (lac  >0) 

Iac*K,  /  r|m  for  generator  mode  (lac<0) 


V*L*J3(  1  ^ 


j  _  ac  ac 

he  -  y 

VdC 


\n e) 


V*T*J 3 


ac  ac 


Vdc 


Vc 


for  motor  mode  (lac  >0) 
for  generator  mode  (lac<0) 


and  r|c  is  the  converter  efficiency 
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4.0  ENERGY  STORAGE  LIBRARY 


The  energy-storage  library  contains  modules  that  can  be  used  to  store  or  deliver  energy  to  an 
HEV .  The  energy-storage  device  is  typically  located  between  a  power  management  controller 
and  the  traction  motor. 

4.1  Lead  Acid  Battery 


INPUTS 

OUTPUTS 

System  Current  Sum  [A] 

DC  Bus  Voltage  [V] 

State-Of-Charge  [  ] 

Constants  as  defined  by  the  mask: 


totnum: 

Eoc: 

cell_r: 

cell_K: 

termRb: 

Ct: 

C5: 

chrg_eff: 

dischrgeff: 

amb_T: 

initECP: 

init_SOC: 

initQ: 


Total  Number  of  Series  Cells  in  System  [  ] 
Open-Circuit  Voltage  of  a  Single  Cell  [V] 

Internal  Resistance  of  a  Cell  [Q] 

Polarization  Resistance  of  a  Cell[Q] 

Terminal  Resistance  of  the  Battery  System  [Q] 
Terminal  Capacitance  of  the  Battery  System  [F] 
Discharge  Capacity  @  5  hour  rate  and  30°C  [A-hr] 
Charge  Efficiency  of  the  Battery  [  ] 

Discharge  Efficiency  of  the  Battery  [  ] 

Ambient  Temperature  in  °C  (-20°C  to  60°C) 

Initial  Polarization  Voltage  [V] 

Initial  State-of-Charge  of  Battery  System  (0-1) 
Initial  Energy  within  the  Battery 
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Auxiliary  Calculations:  None 


Mask: 

1 .  Lead- Acid  Battery  -  ESS 

2.  Lead- Acid  Battery-Series  Cell  Configuration!  [Num  of  Cells,  Eoc(V)  of  cell]:  |Cell 
Defiis[r,Kpr  (ohms)]:  |Term  Defhs[Rb(ohms),Ct(F)]:  |Chrg  Defns[Cap(A-hrs)]:  |[chg 
eff,dischg  eff,temp(C)]:  |[Init  K  V(V),Init  SOC(0-1)]: 

3.  tot_num=@  1(1,1);  Eoc=@l(l,2);  cell_r=@2(l,l);  term_Rb=@3(l,l);  cell_K=@2(l,2); 
Ct=@3(l,2);  C5=@4(l,l);  amb_T=@5(l,3);  chrg_eff=@5(l,l);  dischrg_eff=@5(l ,2); 
init_ECP=@6(l ,  1 );  init_SOC=@6(l,2);  init_Q=C5 * 3 600* ( 1  -@6(  1 ,2)); 

Assumptions: 

■  The  battery  model  is  valid  for  the  temperature  range  -20°C  to  60°C. 

■  The  internal  resistance  of  the  battery  is  considered  constant  throughout  the  battery 
operation.  Typically,  this  varies  with  battery  state  of  charge  and  temperature. 

■  When  modeling  a  collection  of  batteries,  each  battery  is  assumed  to  have  the  same  open 
circuit  voltage  and  state  of  charge.  In  actuality,  this  is  not  true,  since  each  battery  can  be 
at  a  different  state  of  charge  and  therefore  have  a  different  open  circuit  voltage. 

■  Positive  current  entering  the  battery  charges,  negative  current  discharges 

■  The  electrode  polarization  is  not  big. 

■  Electromotive  potential  E  is  the  function  of  SOC  and  temperature. 

■  Polarization  resistance  depends  on  both  the  SOC  and  the  direction  of  the  current. 

■  Capacity  is  the  function  of  temperature  and  current. 

■  Acid  electrolyte  mole  concentration  (if  flooded  battery)  is  5M  when  the  battery  is  fully 
charged. 


The  inputs  and  outputs  of  the  lead-acid  battery,  energy-storage  system  are  shown  in  Figure  26. 


Battery  SOC 
[0-1] 

Figure  26.  Inputs  And  Outputs  Of  Lead-Acid  Battery  Energy-Storage  System 
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The  mask  contains  default  values  for  300V  -  25  12V  cell,  lead-acid  battery  system.  These 
values  are  listed  in  Table  12.  An  example  of  the  mask  window  for  this  module  is  displayed 
in  Figure  27.  The  first  input  line  asks  for  the  number  of  cells  within  the  energy  storage 
system  and  the  open  circuit  voltage  of  a  single  cell.  The  next  input  vector  defines  the  internal 
cell  resistance  and  the  polarization  resistance  of  a  cell.  The  terminal  resistance  and 
capacitance  of  the  energy-storage  system  are  listed  in  next  input  vector.  The  total  capacity  of 
the  energy-storage  system  is  defined  in  the  fourth  vector.  The  fifth  input  vector  contains  the 
values  for  the  ambient  temperature  of  the  system  as  well  as  the  system’s  charge  and  discharge 
efficiencies.  The  initial  state  of  the  energy-storage  system  is  defined  by  the  initial 
polarization  voltage  and  the  initial  state  of  charge. 


Figure  27.  Window  Displayed  When  Mask  Is  Activated  For  The  Lead  Acid 
Battery  -  Energy  Storage  System 
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Table  12.  Default  Values  j 

System  Definitions 

Values 

Total  Number  of  Batteries 

25  cells 

Single  Battery  Open-Circuit  Voltage 

12  V 

Internal  Resistance(r)  of  a  Battery 

2.0E-3  Q 

Polarization  Resistance(K)  of  a  Battery 

2.6E-3  Q 

Terminal  Resistance(Rb)  of  the  Battery  Pack 

0.3  0 

Terminal  Capacitance  of  the  Battery  Pack 

1.0  F 

Discharge  Capacity  @  5  hour  rate  and  30°C  for  a  Battery 

100  [A-hr] 

Charge  Efficiency  of  the  Battery  System 

0.8 

Discharge  Efficiency  of  the  Battery  System 

1 

Ambient  Temperature  in  °C 

30°C 

Initial  Polarization  Voltage(ri) 

0 

Initial  State-of-Charge  of  the  Battery  System(0-1 ) 

6.85  j 

4.2  Flywheel  Battery 


INPUTS 

OUTPUTS 

Input  Current  (A) 

DC  Bus  Voltage  (V) 

Flywheel  State  of  Charge  (0-1) 

Purpose: 

This  block  models  a  flywheel  battery  and  an  associated  internal  controller.  The  model 
incorporates  a  simple  flywheel  model  with  a  steady-state  PM  motor/generator  model  and 
inverter/rectifier.  The  inputs  and  outputs  are  the  same  as  the  chemical  battery  making  this 
device  appear  as  a  controllable  DC  voltage  source.  In  order  to  make  voltage  an  output,  this 
model  includes  a  capacitor  and  resistor  in  parallel  across  the  terminal.  The  choice  of 
capacitance  to  maintain  stability  in  the  bus  voltage  is  dependent  on  the  peak  loads,  and  time 
step  being  used  in  the  simulation.  Note  that  this  capacitance  does  not  relate  directly  to  a 
capacitor  rating  required  on  the  bus.  Though  some  form  of  capacitance  will  be  required,  its 
exact  value  will  require  a  higher  fidelity  simulation  running  at  much  smaller  time  steps  than 
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the  1/10  sec  typically  used  for  this  model.  The  controller  included  in  this  model  is  a  sample 
of  a  simple  technique  for  maintaining  a  specified  voltage  on  the  bus.  Other  control  methods 
may  also  be  desired  and  may  be  implemented  by  the  user  by  unmasking  the  block. 

Parameters: 


Auxiliary  Calculations:  None 


Flywheel  Battery  Mask:  no  mask ) 


4.2.1  Flywheel  Battery  (Roller  Bearings) 


INPUTS 

OUTPUTS 

Voltage  (V) 

Output  Current  (A) 

Current  Command  (A) 

Speed  (RPM) 

Discharge  Efficiency  (0-1) 

Charge  Efficiency  (0-1) 

Purpose: 

There  are  two  major  components  and  one  auxiliary  component  of  this  system: 

1 .  Motor/Generator  and  Converter 

2.  Flywheel 

3.  Efficiency  Calculation 

This  subsystem  uses  a  permanent  magnet  generator.  The  'current  command’  is  the  current 
command  provided  to  the  generator  controller.  The  ‘output  current’  is  the  actual  current 
output  to  the  DC  Bus.  Inputs  for  the  ‘Flywheel  Subsystem’  are  DC  Voltage,  generator 
torque,  generator  inertia,  and  bearing  loads.  The  ‘Flywheel  Subsystem’  has  outputs  of 
auxiliary  current  draw,  stored  mechanical  power,  and  rotor  speed. 
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Parameters: 


Auxiliary  Calculations:  None 

Flywheel  Energy  Storage  Subsystem  1  Mask:  ( no  mask ) 


Purpose: 

There  are  two  major  components  and  one  auxiliary  component  of  this  system: 

1 .  Motor/Generator  and  Converter 

2.  Flywheel 

3.  Efficiency  Calculation 

This  subsystem  uses  a  permanent  magnet  generator.  The  ‘current  command’  is  the  current 
command  provided  to  the  generator  controller.  The  ‘output  current’  is  the  actual  current 
output  to  the  DC  Bus.  Inputs  for  the  ‘Flywheel  Subsystem’  are  DC  Voltage,  generator 
torque,  generator  inertia,  and  bearing  loads.  The  ‘Flywheel  Subsystem’  has  outputs  of 
auxiliary  current  draw,  stored  mechanical  power,  and  rotor  speed. 

Parameters: 

Auxiliary  Calculations:  None 

Flywheel  Energy  Storage  Subsystem  II  Mask:  (no  mask ) 
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Purpose: 

This  block  models  a  flywheel  battery  with  a  simple  feed-forward  control  and  low  gain 
feedback.  The  feed-forward  signal  exactly  balances  the  input  current.  There  is  a  low  gain 
feedback  to  keep  the  bus  voltage  near  its  reference  to  adjust  for  losses. 

Parameters: 

Auxiliary  Calculations:  None 
Flywheel  Battery  System  Mask:  (  no  mask ) 


4.3  Capacitor 


Constants: 

cap  capacitance:  Capacitance  [F] 

cap  resistance:  Self-discharge  resistance  [Q] 

cap_initial_voltage:  Initial  Charge  Voltage  [V] 
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Auxiliary  Calculations:  None 


Mask: 

(1)  Capacitor  Module 

(2)  Capacitor  Module|Capacitance  [F]|Self-Discharge  Resistance  [ohms]|Initial 

Voltage  [V] 

(3)  Capacitance=@l  ;Resistance=@2;Initial_Voltage=@3; 

Assumptions: 

■  For  most  capacitors  the  primary  loss  mechanism  is  the  series  resistance  of  the  capacitor. 
This  has  been  ignored  in  this  model  to  avoid  algebraic  loops  which  could  arise.  There  are 
other  secondary  effects  such  as  hysteresis  which  could  be  modeled,  but  in  general  these 
effects  are  neglected. 

The  following  presents  typical  values  of  the  constants  used  in  the  capacitor  model 


Constant  Value 

cap_capacitance:  20E-3  [F] 

cap  resistance:  2E6  [Q] 

cap_initial_voltage:  230  [V] 
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5.0  VEHICLE  LIBRARY 


Figure  28.  Vehicle  Model  Library 


This  library  contains  models  of  both  tracked  and  wheeled  vehicles.  There  are  two  versions 
for  the  wheeled  vehicle  and  wheel  modules.  The  first  version  contains  a  differentiator  block, 
and  the  second  version  avoids  this  differentaitor  block  through  reformulation  of  the  dynamic 
equations.  Both  versions  have  been  tested  to  produce  the  same  results,  although  version  2 
executes  faster  than  version  1 . 

5.1  Tracked  Vehicle 


INPUTS 

OUTPUTS 

Fuel  Mass(kg) 

Sprocket  (Axle)  Speed  (rpm) 

Vector  of  [Shaft  Torque(N-m), 
Eff  Motor  +  Coupling  Inertia 
(kg*mA2) 

Vehicle  Position  (m) 

Terrain  Gradient  (rad) 

Vehicle  Velocity  (rpm) 

Purpose: 

The  purpose  of  this  block  is  to  simulate  the  dynamics  of  a  tracked  vehicle  by  summing  forces 
from  various  inputs,  including  aerodynamic  drag,  rolling  resistance,  and  drag  due  to  slope. 
This  sum  is  divided  by  the  inertia  of  the  vehicle  to  determine  its  acceleration,  which  is  then 
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integrated  once  for  velocity  and  twice  for  position.  Note  that  this  model  does  not  include 
active  suspension  loads. 


Parameters: 

track_r 

radius  (m)  of  track  drive  sprocket 

trackb 

track  damping  constant  (N/(rad/s))  for  friction  torque 

track_J 

effective  inertia  of  track  (in  vehicle  reference  frame,  Kg*  m2) 

vehmass 

Gross  vehicle  mass  less  fuel  (kg) 

vehCd 

coefficient  of  aerodynamic  drag  of  the  vehicle 

veh_A 

Frontal  area  of  vehicle  (mA2) 

veh_pos_i 

Initial  position  of  vehicle  (m) 

veh_spd_i 

Initial  velocity  of  vehicle  (m/s) 

vehRRm: 

Rolling  resistance  of  vehicle  (velocity  independent)  [N/kg] 

veh_RR_v: 

Rolling  resistance  vs.  velocity  [N/(m/s)] 

sample_N 

max  number  of  samples 

sampleT 

simulation  time  between  samples 

Auxiliary  Calculations:  None 


Vehicle  Mask 

(1)  Tracked-Vehicle  Dynamics 

(2)  Tracked  Vehicle  Dynamics  Block  Definitions\n(Click  on  "Help"  to  list  Inport  and 
Outport  parameters.)!  [Sprocket  Radius  (m),  Damping  Constant  (N*s),  Inertia 
(kg* mA2)] | Vehicle  Mass  (kg)|{RR_m  (N/kg),RR_v  (N*s/m),Cd,Affront]|[Initial 
Position  (m), Initial  Speed  (m/s)]|[sample_N,sample_T] 

(3)  r=@l(l);b=@l(2);J=@l(3);mass=@2;RR_m=@3(l);RR_v=@3(2);Cd=@3(3); 
A=@3  (4);pos_i=@4(  1  );spd_i=@4(2);sample_N=@5(  1  );sample_T=@5(2); 

(4)  Drawing  Commands 

(5)  Tracked- Vehicle  Dynamics  Block  Port  Parameters\n\nInport  1:  Fuel  Mass 
(kg).\nInport  2:  Vector  of  [Shaft  Torque  (N*m),  Effective  Motor+Coupling  Inertia 
(kg*mA2)].\nInport  3:  Terrain  Gradient  (radians) .\n\nOutport  1:  Sprocket  (Axle) 
Speed  (RPM).\nOutport  2:  Vehicle  Position  (m)\nOutport  3:  Vehicle  Velocity  (m/s) 
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Model  Equations: 


This  model  is  appropriate  for  modeling  heavy  vehicles  and  in  particular,  vehicles  with  tracks, 
where  there  is  data  on  the  static  rolling  resistance.  It  assumes  that  there  are  five  forces  acting 
on  the  vehicle: 

The  aerodynamic  force:  Faero  =  Cd*Area*velA2 

The  static  rolling  resistance:  Frrm  =  Rm*vehicle  mass 

The  dynamic  rolling  resistance:  Frrv  =  Rv*vel 

The  gravity  force  (dependent  on  grade  in  direction  of  travel): 

Fgrad  =  vehicle  mass  *  g  *  sin(grad) 

The  effective  thrust  force  of  the  drive  train:  an  input  (see  note  below) 

These  forces  are  summed  and  divided  by  the  effective  mass  of  the  vehicle  to  determine  its 
acceleration.  This  acceleration  can  then  be  integrated  to  determine  the  velocity  of  the 
vehicle. 

The  effective  mass  of  the  vehicle  is  determined  by  the  following  equation,  which  accounts 
for  fuel  usage  and  rotational  inertia  of  the  drive  train. 

meffective  =  mdry  +  mfuel  mefifdrivetrain 

Note  that  the  effective  thrust  available  from  the  drive  train  has  had  all  velocity-dependent 
torques  deducted  in  previous  blocks,  but  has  not  had  any  deduction  for  inertial  effects.  These 
inertial  effects  are  collected  in  this  block  and  accounted  for  in  meffdrivetrain. 


65 


5.2  Wheeled  Vehicle  (Version  1) 


VEHICLE 
Ver  1.0 


Inputs 

Outputs 

Vehicle  Velocity  [m/s] 

Total  Load  on  Vehicle  [N] 

Road  Grade  [radian] 

Maximum  traction  Force  [N] 

Fuel  Mass  [kg] 

Total  Vehicle  Mass  [kg] 

Vehicle  Acceleration  [m/s2] 

Constants: 


vehairdensity: 

veh_aero_drag_coef: 

veh_front_area: 

veh_normal_wind: 

veh_parallel_wind : 

veh_coef_yaw_sens: 

veh_drag_coef: 

veh_fs: 

vehfO: 

Timestepsmax: 


Density  of  the  air  that  the  vehicle  is  exposed  to  [kg/m3] 

Drag  coefficient  of  vehicle  [  ] 

Frontal  area  of  vehicle  [m2] 

Wind  velocity  normal  to  vehicle  motion  [m/s] 

Head  wind  velocity  [m/s] 

Coefficient  of  yaw  sensitivity  of  the  vehicle  during  motion  [  ] 
Aerodynamic  drag  coefficient  [  ] 

Vehicle  speed  dependent  rolling  coefficient  [  ] 

Vehicle  speed  independent  rolling  coefficient  [  ] 

Maximum  number  of  the  array  size 

of  variables  saved  to  workspace  during  simulation 


Auxiliary  Calculations: 


yout_veh: 

vehRRtractloss: 
veh_RR_regen_loss : 
vehaeroregenloss: 

veh_aero_tract_loss : 

veh_inertial_tract_energy: 

veh_inertial_regen_energy: 


An  array  consisting  of  two  vectors:  actual  vehicle 
velocity  [km/hr]  and  distance  traveled  [km] 

Rolling  resistance  energy  loss  during  acceleration  [kWh] 
Rolling  resistance  energy  loss  during  deceleration  [kWh] 
Aerodynamic  resistance  energy  loss  during  acceleration 
[kWh] 

Aerodynamic  resistance  energy  loss  during  deceleration 
[kwh] 

Total  energy  required  to  accelerate  vehicle  [kwh] 

Total  recovery  energy  available  during  deceleration 
[kwh] 


Mask: 

(1)  Vehicle  Module 

(2)  Vehicle  Module  (SwRI)|[Air  Density  [kg/m3], Vehicle  Aero  Drag  Coeff  [ 
]]| Vehicle  frontal  Area  [m2]| [Normal  Wind  Speed  [m/s],  Parallel  Wind  Speed 
[m/s]]|[  Speed 

Effect  Friction  Coef  [],road  friction  coefficient  [j.  Vehicle  coef  of  yaw  sensitivity 
[  ]  |  Vehicle  Mass  [kg]|Maximum  Time  Steps  [  ] 
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(3)  veh_air_density=@  1  [  1  ] ;  veh_aero_drag_coef^@  1  [2];veh_front_area=@2; 

veh_normal_wind=@3  [  1  ];  veh__parallel_wind=@3  [2] ;  veh_fs=@4[  1  ] ;  veh_fO=@ 
4[2];veh_coef_yaw_sens=@4[3];veh_mass=@5;timesteps_max=@6;  ~ 


Assumptions: 

■  The  total  load  acting  on  the  vehicle  is  the  sum  of  all  the  driving  resistance  forces.  These 
consist  of  the  aerodynamic  drag,  tire  rolling  resistance,  and  road  grade  dependent  loads.  In 
the  calculation  of  aerodynamic  forces,  the  effect  of  head  wind  and  normal  wind  components 
have  been  included.  In  the  computation  of  the  tractive  resistance,  the  effect  of  the  vehicle 
speed  on  the  traction  coefficient  has  been  considered. 

m  The  total  vehicle  mass  includes  the  weight  of  the  fuel,  which  decreases  as  the  heat  engine 
produces  power. 


The  following  presents  typical  values  used  for  the  constants  in  the  vehicle  model 

Constant 

Value 

vehairdensity: 

1.23  [kg/m3] 

vehaerodragcoef: 

0.4  [] 

veh_ffont_area: 

3.0  [m2] 

vehnormalwind: 

0.0  [m/s2] 

veh_parallel_wind: 

0.0  [m/s2] 

veh_coef_yaw_sens: 

0.002  [] 

vehdragcoef: 

0.4  [] 

veh  fs: 

0.005  [] 

vehfD: 

0.01  [] 

Timestepsmax: 

10000 [] 

5.3  Wheeled  Vehicle  (Version  2) 


INPUTS 

OUTPUTS 

Driving  force  into  vehicle  [N] 

Velocity  of  vehicle  [m/s] 

Mass  of  fuel  consumed  [kg] 

Brake  fraction  [  ] 

Wheel  mass  [kg] 

Road  grade  [rad] 

> 

> 

VEHICLE 

> 

Ver2.0 

> 

> 

> 
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Constants: 


same  as  Wheeled  Vehicle  (Version  1)  (see  Section  5.2) 


Auxiliary  calculations: 

veh  1  _accel_hist:  V ehicle  acceleration  history 

vehl_speed_hist  Vehicle  speed  history 

vehl_position_hist  Vehicle  position  history 

Mask: 

(1)  Vehicle  Module 

(2)  Vehicle  Module|Air  Density  (kg/mA3),  Aero  Drag  Coeff  (  )|Vehicle  Frontal  Area 
(mA2)|Normal  Wind  Speed  (m/s), Parallel  Wind  Speed  (m/s)|Spd-dep  Rolling  Coef,V-indep 
Roll  Coeff, Yaw  Sens  Coeff|Vehicle  Chassis  Mass  (kg),  Road/Tire  Friction  Coefficient  [ 
]|Record  Variable  Quantity  and  Time  (#,sec) 

(3) veh_air_density=@l(l);veh_aero_drag_coef=@l(2);veh_front_area=@2;veh_normal_wind 
=@3(l);veh_parallel_wind=@3(2);veh_fs=@4(l);veh_f0=@4(2);veh_coef_yaw_sens=@4(3); 
veh_mass=@5(l  );roadtire_frict_coef=@5(2);sample_N=@6(  1  );sample_T=@6(2); 

Assumptions:  Same  as  Wheeled  Vehicle  (Version  1)  (see  Section  5.2) 


For  typical  values  for  the  constants,  please  refer  to  Section  5.2. 

5.4  Wheel  (Version  1) 


> 

> 

Wheel 

> 

> 

Ver  1.0 

> 

> 

> 

inputs 

Outputs 

Driving  force  into  vehicle  [N] 

Velocity  of  vehicle  [m/s] 

Mass  of  fuel  consumed  [kg] 

Brake  fraction  [] 

Wheel  mass  [kg] 

Road  grade  [rad] 

Constants: 

whll_tr: 

whll_mass: 

motor_N_max: 

timestepsmax: 


Wheel  Tire  Radius  [m] 

Wheel  mass  [kg] 

Maximum  motor  speed  [rpm] 

Maximum  number  of  the  array  size 

of  variables  saved  to  workspace  during  simulation 
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Auxiliary  Calculations: 

servicebrkenergy:  Energy  loss  during  application  of  the  mechanical  brake  [kWh] 
whl  tract  energy:  Energy  delivered  to  the  wheels  during  acceleration  [kWh] 

whl_regen_energy:  Energy  recovered  by  the  vehicle  during  deceleration  of  wheel 

[kWh] 


Mask: 

(1)  Wheel  Module 

(2)  Wheel  Module  (SwRI)|Wheel  Tire  Radius  [m]|Wheel  Tire  Mass  [kg]|Max 
Motor  Speed  [rpm]|Max  Timesteps  [  ] 

(3)  whl_tr=@l  ;whl_mass=@2;motor_N_max=@3;timesteps_max=@4; 

Assumptions: 


The  vehicle  velocity  and  acceleration  are  based  on  the  assumption  that  the  entire  mass  of  the 
vehicle  is  concentrated  at  the  wheel  center. 


The  following  presents  typical  values  used  for  the  constants  in  the  wheel  model. 


Constant 

whll_tr: 

whll_mass: 

motorNmax: 

timesteps_max: 


Value 
0.3  [m] 

25  [kg] 
7000  [rpm] 
10000 


5.5  Wheel  (Version  2) 


> 

> 

Wheel 

> 

> 

Ver  2.0 

> 

> 

> 

INPUTS 

OUTPUTS 

Input  Torque  [Nm] 

Driving  Force  [N] 

Coupler  Inertia  [kg  m2] 

Rotational  Speed  of  Wheel  [rad/sec] 

Vehicle  Velocity  [m/s] 

Equivalent  mass  [kg] 

Constants: 

whll_tr:  Wheel  Tire  Radius  [m] 

whll  mass:  Wheel  mass  [kg] 
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Auxiliary  Calculations:  None 


Mask: 

1.  WHEEL  MODULE 

2.  WHEEL  MODULE  (SwRI)\n\n  T1,V,J1  (in)  <=>  F,N2,m  (out)|Wheel  System's  Tire  Radius 
[m]|  Total  Wheel  Mass  (Wheel, Tire,  and  Brake  Sys)  [kg] 

3.  wheel_radius  =  @1;  wheel  mass  =  @2; 


Assumptions: 


The  vehicle  velocity  and  acceleration  are  based  on  the  assumption  that  the  entire  mass  of  the 
vehicle  is  concentrated  at  the  wheel  center. 


5.6  Driver 


INPUTS 

OUTPUTS 

Desired  Speed  [m/s] 

Driver  Power  Demand  [  ] 

Actual  Speed  [m/s] 

Constants: 


dvr_proport_gain : 
dvr_integraUgain: 
dvr_deriv_gain: 


Proportional  gain  of  PID  controller  within  driver  model  [  ] 
Integral  gain  of  PID  controller  within  driver  model  [  ] 
Derivative  gain  of  PID  controller  within  driver  model  [  ] 


Auxiliary  Calculations: 

yout  dvr:  vector  containing: 

1 .  PID  controller  output  (driver  power  demand)  [  ] 

2.  Desired  speed  [m/s] 

3.  Actual  speed  [m/s] 
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Mask:  None 
Assumptions: 


The  driver  is  modeled  as  a  PID  controller  with  fixed  gains.  The  controller  output  is  dependent 
on  its  gains  and  the  error  between  the  desired  and  actual  vehicle  speed. 


The  following  shows  typical  values  for  the  constants  of  the  driver  model. 


Constants  Value 

dvr_proport_gain:  0.2  [  ] 

dvr_integral_gain:  0.001  [] 

dvr  deriv^gain:  0  [  ] 
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6.0  CONTROLLER  LIBRARY 


System 

Controller 


MHi 

> 

mSm 

> 

Energy  Management 
controller 

>  Controller 

>  (lAT’s  demo) 

(SwRI's  demo) 

> 

> 

1 

► 


APU  Controller 
(Simplified) 


Figure  29.  Controller  Library 


The  controller  library  contains  examples  of  energy-management  controllers,  which  are 
dependent  on  the  HEV  application  at  hand.  Therefore,  these  controllers  are  provided  as 
examples  of  possible  energy-management  controllers,  which  the  user  can  use  as  a  basis  to 
develop  custom  power  management  controllers.  This  library  also  includes  controllers  for  the 
energy  storage  device  and  the  APU. 

6.1  Energy  Management  Controller  (I) 


SwRI 

System 

Controller 


INPUTS 

OUTPUTS 

Generator  Output  Current  [Amps] 

Battery  Current  Command  [Amps]  i 

Drive  Motor  Current  Draw  [Amps] 

Engine  Throttle  Command  [0-1] 

System  Voltage  [Volts] 

Generator  Current  Command  [Amps] 

ESS  Current  Requirement  [Amps] 

Engine  ON/OFF  [1/0] 

*ESS  State  of  Charge  [  ] 

Engine  Speed  [rpm] 

|  *ESS:  Energy  Storage  System  | 

Constants:  None 


Auxiliary  Calculations:  None 


Mask:  None 


Assumptions:  Specific  to  application 
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6.1.1  Energy  Storage  System 


1  Controller  f 


Energy  Storage 
System  Controller 


INPUTS 

OUTPUTS 

*ESS  State  of  Charge 

ON/OFF  signal  to  “APU  [1/0] 

*ESS:  Energy  Storage  System 
“APU:  Alternate  Power  Unit 

Constants:  None 


Auxiliary  Calculations:  None 

Mask:  None 


Assumptions:  Specific  to  application 


6.1.2  APU  Controller 


APU 

Controller 

> 


■> 

> 

> 

> 

.> 


|  INPUTS 

OUTPUTS 

ESS  State  of  Charge  [  ] 

Engine  Throttle  Command  10-11 

Engine  Speed  [rpm] 

Generator  Current  Command  [Amps] 

Engine  ON/OFF  [1/0] 

Constants:  None 


Auxiliary  Calculations:  None 


Mask:  None 


Assumptions:  Specific  to  application 
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6.2  System  Controller 


p 

1 

| 

X 

X 

£ 

II 

1 

1 

i 

I 

1 

I 

X 

s 

Cor 

ystem 

itrollers 

INPUTS 

OUTPUTS 

Generator  Output  Current  (A) 

Engine  Throttle  Command 

Drive  Motor  Output  Current  (A) 

Generator  Current  Command 

System  Voltage  (V) 

Battery  Power  Command 

ESS  Percent  Power  Command 

Engine  On/Off 

Engine  Speed  (RPM) 

Parameters: 


Auxiliary  Calculations:  None 
System  Controllers  Mask:  None 
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7.0  COUPLING  LIBRARY 


The  coupling  library  contains  mechanical  couplers  representing  shaft-to-shaft  connections 
between  blocks.  These  couplings  allow  for  a  gear  ratio  and  an  efficiency  of  power  transfer. 
Couplers  are  typically  located  between  the  engine  and  the  generator,  as  well  as  between  the 
traction  motor  and  the  wheel. 

7.1  Engine  Generator  Coupling 


Constants: 

cplspeedratio:  Ratio  of  input  shaft  speed  to  output  shaft  speed  of  the  coupler 

cpl_efficiency:  Efficiency  of  coupler 

Auxiliary  Calculations:  None 


Mask:  None 
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7.2  Motor  Wheel  Coupling  (I) 


> 

> 

Spd 
j  Rdcr 

> 

Ver  1.0 

> 

> 

> 

inputs  ~~ 

Outputs 

Input  Torque  [Nm] 

Output  Torque  [Nm] 

Output  Side  Rotational  Speed  [rpm] 

Input  Side  Rotational  Speed  [rpm] 

Inertia  of  the  driving  component 

Axle  inertia  at  input  [kg  -m2] 

Axle  inertia  at  output  [kg  -m2] 

Constants: 

axll_gear_ratio: 

axIlJl: 

axll_J2: 

axllefficiency: 


Axle  Gear  Ratio  (Output  Gear  Radius/Input  Gear  Radius)  [  ] 
Axle  Inertia  at  input  [kg-m2] 

Axle  Inertia  at  output  [kg-m2] 

Axle  Efficiency  from  input  to  output  [  ] 


Auxiliary  Calculations:  None 


Mask: 

(1)  Axle  Module 

(2)  Axle  Module  (SwRI)|Axle  Gear  Ratio  [  ]|Axle  Input  Side  Inertia  [kg-m2]|Axle 

Output  Side  Inertia  [kg-m2]|Axle  Efficiency  [  ] 

(3)  axl_gear_ratio=@  1  ;axl_Jl=@2;axl_J2=@3;axl_efficiency=@4; 

Assumptions: 

■  Internal  efficiencies  between  the  individual  gears  within  the  axle  subsystem  have  not  been 
considered.  These  are  modeled  by  the  variable  axl_efficiency,  which  represents  the  overall 
inefficiency  from  the  input  to  output  shaft. 

■  The  axle  module  dynamics  (see  Section  5.1,  Technical  Notes)  is  based  on  the  analysis  of 
two  spur  gears,  one  on  the  input  side  and  the  other  on  the  output  side.  The  dynamic  model 
will  need  additional  enhancements  if  more  complex  gear  systems  such  as  planetary  gears 
are  to  be  considered  for  the  axle. 

■  Backlash  effects  within  the  gears  have  not  been  considered  in  the  development  of  the  axle 
model. 
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The  following  presents  typical  values  used  for  the  constants  in  the  axle  model 


Constant 

axll_gear_ratio: 

axllJl: 

axll_J2: 

axll_efficiency: 

axl_mass: 


Value 

2[] 

10  [kg-m2] 
10  [kg-m2] 
96.13  % 
100  kg 


7.3  Motor  Wheel  Coupling  (II) 


Speed 
Rdcr  > 

Ver  2.0  > 


INPUTS 

OUTPUTS 

Input  Torque  [Nm] 

Output  Torque  [Nm] 

Wheel  Rotational  Speed  [rpm] 

Input  Side  Rotational  Speed  [rpm] 

Inertia  of  the  driving  component  [kgm2] 

Output  Equivalent  Inertia  [kg-m2] 

Power  Consumed  [kg-m2] 

Constants: 

axll_gear_ratio: 

axllJl: 

axll_J2: 

axll_efficiency: 


Axle  Gear  Ratio  (Output  Gear  Radius/Input  Gear  Radius)  [  ] 
Axle  Inertia  at  input  [kg-m2] 

Axle  Inertia  at  output  [kg-m2] 

Axle  Efficiency  from  input  to  output  [  ] 


Auxiliary  Calculations:  None 


Mask: 


(1)  Axle  Module 

(2)  AXLE  MODULE  (SwRI)\n\nTl,N2,Jl  <=>  T2,Nl,J2|Axle  Gear  Ratio  [  ]|Axle 
Efficiency  [  ]| Axle  Input  Gear  Inertia  [kg*m2]|Axle  Output  Gear  Inertia  [kg*m2] 

(3)  gear_ratio  =  @1  ^efficiency  =  @2;ingear_inertia  =  @3;outgear_inertia  =  @4 


Assumptions:  See  Section  7.2. 
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8.0  MISCELLANEOUS  LIBRARY 


DC  to  DC 
Converter 


Roller 

Bearing 


AC  to  AC 
Converter 


Magnetic 

Bearing 


Roller  Bearing  Act.  Mag  Bearing 


AC  to  DC 
Converter 


Rectifier 


>  Windage 
Windage 


DC  to  AC 
Converter 


Figure  31.  Miscellaneous  Components  Libra 


This  library  contains  electrical  couplings  that  convert  electrical  power  between  separate  buses. 
These  couplings  include  variations  in  voltage  and  account  for  power  transfer  inefficiencies. 
Other  bearing  models  such  as  roller  and  magnetic  bearings  are  also  included  in  this  library. 

8.1  DC-DC  Converter 

> 

DC  to  DC 
>  Converter 


Purpose: 

The  DC-DC  Converter  block  models  a  solid  state  converter  that  inverts,  transforms,  and 
rectifies  electricity  between  two  DC  buses.  This  model  is  a  simple  power  transformer  with 
efficiency  that  accounts  for  the  direction  of  power  flow.  The  frequency  of  the  inverted  power 
is  assumed  to  be  high  to  reduce  ripple  effects  on  the  DC  bus  with  small  capacitors.  The 
power  required  to  be  moved  from  the  high  voltage  bus  to  the  low  voltage  bus  is  input  as  a 
function  of  time,  and  the  model  determines  the  associated  DC  currents  on  each  bus.  Note  the 


INPUTS 

OUTPUTS 

Hiqh  Voltaae  Power  Command  (\N) 

High  Voltage  (V) 

Hiah  Side  Current  (A) 

Low  Side  Current  (A) 

Low  Voltage  (V) 
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sign  convention  on  the  currents  and  power  variables.  The  controller  that  supplies  the  power 
signal  in  watts  must  insure  that  the  bus  voltages  are  within  desired  ranges.  Note  that  the 
currents  are  determined  to  conserve  power,  and  may  become  large  as  voltages  are  decreased 
by  this  and  other  load  paths. 

Parameters: 

dcdc  l  eff  Efficiency  of  DC/DC  Converter  (0- 1 ) 

Auxiliary  Calculations:  None 

DC-DC  Mask 

(1)  DC  to  DC  Converter 

(2)  DC  to  DC  Converter  Block  Definitions\n(Click  on  "Help"  to  list  Inport  and  Outport 
parameters.)|Efficiency  (0-1) 

(3)  eff=@l; 

(4)  DC  to  DC\nConverter 

(5)  DC  to  DC  Converter  Block  Port  Parameters\n\nInport  1 :  High  Voltage  Power 
Command  (W).\nInport  2:  High  Voltage  (V)Anlnport  3:  Low  Voltage  (V).\n\n 
Outport  1 :  High  Side  Current  (A).\nOutport  2:  Low  Side  Current  (A). 

Model  Equations: 

Im  =  Ph/V M  Ihj  >  0  is  current  into  the  converter 

1,0  ~  Ih,*^  when  P  >  0  (power  flows  from  hi  to  lo  voltage) 

Ihi*(l/  ^l)  when  P  <  0 

■  note  that  Ilo>0  represents  current  flow  out  of  the  converter 
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8.2  AC-AC  Converter 


AC  to  AC 
Converter 


INPUTS 

OUTPUTS 

High  Voltage  Power  Command  (W) 

High  Side  Current  (A_rms) 

High  Voltage  (VrmsJI) 

Low  Side  Current  (A_rms)  j 

Low  Voltage  (VrmsJI) 

Purpose: 

The  AC-AC  Converter  block  models  a  solid  state  or  conventional  iron  transformer  for 
changing  the  voltage  between  two  AC  buses.  This  model  is  a  simple  power  transformer  with 
efficiency,  which  accounts  for  the  direction  of  power  flow.  It  follows  the  same  convention 
used  in  the  DC-DC  converter.  The  power  required  to  be  moved  from  the  high  voltage  bus  to 
the  low  voltage  bus  is  input  as  a  function  of  time,  and  the  model  determines  the  associated 
DC  currents  on  each  bus.  Note  the  sign  convention  on  the  currents  and  power  variables.  The 
controller,  which  supplies  the  power  signal  in  watts,  must  insure  that  the  bus  voltages  are 
within  desired  ranges.  Note  that  the  currents  are  determined  to  conserve  power  and  may 
become  large  as  voltages  are  decreased  by  both  this  and  other  load  paths.  The  power 
transferred  through  this  model  assumes  that  both  buses  are  at  a  Power  Factor  of  1 .  Later 
modifications  could  allow  for  more  detailed  treatment  of  not  PF=1  conditions. 

Parameters: 

aca_eff  Efficiency  of  AC/AC  Converter  (0-1) 

Auxiliary  Calculations:  None 

AC-AC  Mask 

(1)  AC  to  AC  Converter 

(2)  AC  to  AC  Converter  Block  Definitions\n(Click  on  "Help"  to  list  Inport  and  Outport 
parameters.)|Efficiency  (0-1) 
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(3)  eff=@l; 

(4)  AC  to  AC\nConverter 

(5)  AC  to  AC  Converter  Block  Port  Parameters\n\nInport  1 :  High  Voltage  Power 
Command  (W)\nInport  2:  High  Voltage  (Vrms_ll).\nInport  3:  Low  Voltage 

(V rms_ll) .\n\nOutport  1 :  High  Side  Current  (A_rms).\nOutport  2:  Low  Side  Current 
(Arms). 

Model  Equations: 


Ihi  ~  Ph/(sqrt(3)*Vhi)  Ihj  >  0  is  current  into  the  converter 

Iio =  Ihi*1!  when  P  >  0  (power  flows  from  hi  to  lo  voltage) 

Ihi*(l/  r|)  when  P  <  0 

■  note  that  I,o>0  represents  current  flow  out  of  the  converter 


8.3  AC-DC  converter 


/GDCRc 


INPUTS 

OUTPUTS 

Input  DC  Bus  Voltage  [Volts] 

Output  DC  Bus  Voltage  [Volts]  | 

Magnetizing  Per  Unitized  Current  [  ] 

Output  Magnetizing  Per  Unitized  Current  [  ] 

Per  Unitized  Stator  Current  along  q  axis  [  ] 

Output  Per  Unitized  Stator  Current  along  q  axis  [] 

Per  Unitized  Stator  Voltage  along  q  axis  [  ] 

Output  Per  Unitized  Stator  Voltage  along  q  axis  [  1 

Per  Unitized  Stator  Current  along  d  axis  [  ] 

Output  Per  Unitized  Stator  Current  along  d  axis  \  1 

Per  Unitized  Stator  Voltage  along  d  axis  [  ] 

Output  Per  Unitized  Stator  Voltage  along  d  axis  [] 

Auxiliary  Calculations:  None 

Mask:  None 

Assumptions: 

■  The  buck  rectifier  type  is  valid  for  simulations  that  require  positive  DC  rectifier  command. 

■  Rectifier  harmonic  effects  are  not  modeled. 

■  The  rectifier  is  assumed  to  have  a  fixed,  user-definable  efficiency. 
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8.4  DC-AC  Converter 


SBBcdJne 


INPUTS 

OUTPUTS 

Input  Bus  DC  Voltage  [Volts] 

Output  Per  Unitized  Stator  Voltage  in  the  q  axis  [  ] 

Per  Unitized  Synchronous  magnetic  field  speed  [  ] 

Output  Per  Unitized  Stator  Current  in  the  q  axis  [  ] 

Per  Unitized  Stator  Voltage  in  the  q  axis  [  ] 

Output  Per  Unitized  Stator  Voltage  in  the  d  axis  [  ] 

Per  Unitized  Stator  Current  in  the  q  axis  [  ] 

Output  Per  Unitized  Stator  Current  in  the  d  axis 

Per  Unitized  Stator  Voltage  in  the  d  axis  [] 

Output  Per  Unitized  Rotor  Current  in  the  q  axis  [  ] 

Per  Unitized  Stator  Current  in  the  d  axis  [  ] 

Output  Per  Unitized  Rotor  Current  in  the  d  axis  [  ] 

Per  Unitized  rotor  Current  in  the  q  axis  [  ] 

Per  Unitized  rotor  Current  in  the  d  axis  [  ] 

Auxiliary  Calculations:  None 


Mask:  None 


Assumptions: 

■  The  model  is  a  fully  controlled  boost  inverter. 

■  The  model  does  not  include  inverter  harmonic  effects. 

8.5  Capacitor 


Filter  Cap. 


INPUTS 

OUTPUTS 

Current  (A) 

Voltage  (V) 

Purpose:  To  model  a  parallel  capacitor  and  resistor. 

Parameters: 


CM_Cap  Capacitance  (F) 

CM_R.es  Self-Discharge  Resistance  (Ohms) 

CM_V  Initial  Voltage  (V) 
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Auxiliary  Calculations:  None 


Capacitor  Module  Mask 


(1)  Capacitor  Module 

(2)  Capacitor  Module| Capacitance  [F]| Self-Discharge  Resistance  [Ohms]|Initial  Voltage  [V] 

(3)  Capacitance=@l  ;Resistance=@2;Initial_Voltage=@3; 

(4)  Drawing  Commands 

(5)  Input  is  Current  in  Amps\n  Output  is  Voltage  in  Volts 


Model  Equations: 


dV 

dt 


8.6  Roller  Bearing 


> 

> 

> 

Roller 

Bearing 

> 

Roller  Bearing 

INPUTS 

OUTPUTS 

Radial  Force  (N) 

Torque  (N-m) 

Thrust  Force  (N) 

Bearing  Speed  (RPM) 

Purpose: 

To  models  a  roller  bearing.  Drag  torque  is  calculated  based  on  thrust  load,  radial  load,  and 
viscous  drag.  An  example  follows  showing  extraction  of  parameters  from  typical 
manufacturer’s  data. 


Parameters: 


RBrad 
RBthr 
RB  vis 


Radial  Load  Coefficients  [Tr  (N-m),  Fro  (N),  Er  (Exponent)] 
Thrust  Load  Coefficients  [Tt  (N-m),  Fto  (N),  Et  (Exponent)] 
Viscous  Drag  Coefficients  [Tv  (N-m),  No  (RPM),  Ev  (Exponent)] 
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Auxiliary  Calculations:  None 


Roller  Bearing  Mask 


(1)  Roller  Bearing 

(2)  T  =  Tr*(Fr/Fro)AEr  \n  +  Tt*(Ft/Fto)AEt\n  +  Tv*(N/No)AEv|Radial 
CoefFicients[Tr(N-m),  Fro(Force),  Er(Exponent)]|Thrust  Coefficients[Tt(N-m), 
Fto(Force),  Et(Exponent)]|Viscous  Coefficients[Tv(N-m),  No(RPM),  Ev(Exponent)] 

(3)  RE=@  1  (3);RK=@  1  ( 1  )/@  1  (2) ARE; AE=@2(3); AK=@2(  1 )/ @2(2) AAE; 

VE=@3  (3 ) ;  VK=@3  ( 1  )/@3  (2) A VE; 

(4)  Roller\nBearing 

(5)  This  bearing  model  calculates  the  drag  due  to  normal  ball  bearings.  Normally  radial  and 
thrust  exponents  (Er  &  Et)  are  near  1 .5,  viscous  exponent  (Ev)  is  usually  2/3.  Model 
uses  the  absolute  value  of  the  input  forces  and  speed  for  calculations.\n\nInput  lis  Radial 
Force  (Fr)\nlnput  2  is  Thrust  Force  (Ft)\nlnput  3  is  Speed  (N)\n\nOutput  is  Bearing  Drag 
inN-m 


Model  Equations: 


Example  of  Extraction  of  Inputs  from  Manufacturers  Data 


Table  13.  Typical  Manufacturer’s  Data 

Radial  Force 

Radial  Torque 

Thrust  Force 

Thrust  Torque 

Rotor  Speed 

Lube  Torque 

(lbs) 

(N-m) 

(lbs) 

(N-m) 

(RPM) 

(N-m) 

40 

0.0013 

0.0042 

1000 

0.0962 

50 

0.0018 

0.0058 

2000 

0.1528 

70 

0.0030 

0.0092 

EH 

0.2425 

100 

0.0051 

0.0153 

■H 

0.3178 

130 

0.0075 

0.0222 

8000 

0.3850 

180 

0.0121 

0.0352 

12000 

0.5045 

240 

0.0185 

0.0527 

16000 

0.6111 

350 

1 

450 

1 

600 

0.0729 

0.1897 

820 

0.1172 

0.2923 

1100 

0.1832 

0.4380 
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Choose  Nominal  Force  of  100  lbs,  Maximum  of  1 100  lbs  and  Speed=16000  RPM 


Fro=1001bs  Tro=.0051N-m 

£  _  *n(^ma xAra)  _  ln(.1832/.Q05 1)  _ 
*n(^ma x/^ro)  ln(l  100/100) 

Ft0  =100  lbs  Tto=.0153N-m 

E  -  ln(rn,a JTJZ  ln(.438/.0153) 

'  ln(Fmax/Fto)  ln(l  100/100)  * 

No= 16000  Tv=.611 

_  Infe/rJ  ln(.2425/.61l) 
v  ln(^/^)  ln(4/16)  ' 


8.7  Magnetic  Bearing 


> 

> 

Magnetic 

Bearing 

> 

> 

Act.  Mag  Bearing 

INPUTS 

OUTPUTS 

Force  (N) 

Electric  Power  (\N) 

Speed  (RPM) 

Drag  Torque  (N-m) 

Purpose: 

The  magnetic  bearing  approximately  models  an  active  magnetic  bearing.  Bearings  are 
assumed  to  have  a  defined  relationship  between  air  gap  flux  and  force.  This  model  linearly 
interpolates  the  electric  losses  and  rotating  losses  for  the  given  input  force  at  the  base  speed. 
The  rotating  losses  are  assumed  to  be  due  to  eddy  current  and  hysteresis  losses  in  the 
magnetic  core  of  the  rotor.  The  eddy  current  losses  will  be  proportional  to  speed  squared, 
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while  the  hysteresis  losses  will  be  proportional  to  the  speed.  This  allows  scaling  for  speed 
variations  of  these  rotating  torques  to  account  for  speed  effects. 


Parameters: 


MB_F 

MB_elc_loss 

MB_rot_loss 

MB_cur_rot_loss 

MB_spd 


Vector  of  Input  Forces  (N) 

Vector  of  Electric  Losses  (W) 

Vector  of  Hysteresis  Rotational  Losses  (W)  at  Base  Speed 
Vector  of  Eddy  Current  Rotational  Losses  (W)  at  Base  Speed 
Base  Speed  (RPM) 


Auxiliary  Calculations:  None 


Magnetic  Bearing  Mask 


(1)  Mag  Bearing 

(2)  Active  Magnetic  Bearing  Model|Input  Force|Electric  Losses  (W)|Hysteresis  Rotational 
Losses  (W)  at  Base  Speed|Eddy  Current  Rotational  Losses  (W)  at  Base  Speed|Base 
Speed  (RPM) 

(3)  Input=@l;Electric=@2;Hysteresis=@3/@5*30/pi;Eddy=@4/@5/@5*30/pi; 

(4)  MagneticVnBearing 

(5)  All  losses  are  calculated  as  curve  fits  as  a  function  of  input  force.  Electric  power  is 
assumed  independent  of  speed.  Hysteresis  losses  are  assumed  to  increase  linearly  with 
speed.  Eddy  current  losses  are  assumed  to  increase  as  the  speed  squared.  Hysteresis  and 
eddy  current  losses  appear  as  mechanical  drag  torque.\n\nInput  1  is  Force\nInput  2  is 
Speed  in  RPM\n\nOutput  1  is  Electric  Power  Draw  in  Watts\nOutput  2  is  Drag  Torque  in 
N-m 


Model  Equations: 


Pout  ~  Fit  (F,  Pelech  Fin) 
Hf  =  Fit  (F,  Physt,  Fin) 
Ef  =  Fit  (F,  Peddy>  Fin) 


TO)  =  H, 


\N.j 


+  E, 


/  \2 

N_ 

UJ 
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where 


Pout  is  the  required  electrical  power 
x  is  the  drag  torque 
F  is  the  input  force 
N  is  the  input  rotational  speed 

Fin  is  the  array  specified  in  the  mask  under  'Input  Force' 

Pelect  is  the  array  specified  in  the  mask  under  'Electric  Losses...' 
Physt  is  the  array  specified  in  the  mask  under  'Hysteresis...' 

Peddy  is  the  array  specified  in  the  mask  under  'Eddy  Current. . .' 
Fit(V,Y,X)  is  a  linear  interpolation  of  the  'Y'  vs.  'X'  data  at  x  value  'V' 


8.8  Windage 


W  Windage 


Windage 


INPUTS 

OUTPUTS 

Speed  (RPM) 

Drag  Torgue  (N-m) 

Purpose:  To  calculate  windage  as  a  function  of  speed  for  a  constant  vacuum  pressure. 


Parameters: 

W_spd  Base  Speed  (RPM) 

W_drag_spd  Drag  at  Base  Speed  (W) 

W_De  Speed  Exponent  for  Drag  Power 

Auxiliary  Calculations:  None 


Windage  Mask: 

(1)  Windage  Module 

(2)  Windage  Equation\nP=Po(N/No)ADe|Base  Speed  (RPM)  [No]|Drag  at  Base  Speed 
(W)  [Po]  |  Speed  Exponent  for  Drag  Power  [De] 

(3)  Drag_Const=@2  *  3  0/pi/ (@  1  A@3);Drag_Exp=@3  - 1 ; 

(4)  Windage 

(5)  Models  rotor  windage  assuming  a  constant  vacuum  level.  The  exponent  (De)  will 
normally  vary  from  De=2  at  reasonable  vacuum  (<.001  torr)  to  De=3  at  poor  vacuum 
(>1  torr).\n\nInput  is  speed  in  RPM\nOutput  is  drag  torque  in  N-m 
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Model  Equations: 


ICO  = 


8.9  Driving  Profile 


INPUTS 

OUTPUTS 

Driving  cycle 

Road  Gradient 

Head  Wind  data 

The  driving  profile  model  provides  the  environmental/operating  inputs  to  the  simulation.  These 
are:  driving  cycle  in  which  vehicle  velocity  is  described  as  a  function  of  time,  road  gradient,  and 
head-wind  information  as  a  function  of  time. 

Auxiliary  Calculations:  None 

Mask:  None 

Assumptions:  Application  dependent 
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LIST  OF  ALL  VARIABLES 


Variable 


axil _gear  ratio 


shft Jog  name 


Table  14.  Motor  Wheel  Coupling 
_  Description 


Inertia  at  input  side  of  coupling 


Inertia  at  output  side  of  coupling 


_ Coupling  efficiency 


Coupling  gear  ratio 


_ Mass  of  coupling 


_ Name  of  coupling 


Weight  of  coupling 


Name  of  coupling  shift  logic 


Units  |  Type 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Variable 


c  ft  m 


c_hp_kW 


c_inHg_Pa 


c  in  m 


c_lb_kg 


c  Ibf  N 


c_psi_Pa 


c_rpm_radps 


Table  15.  Constants 


_ Description 


Conversion  factor  from  bar  to  Pascal 


Conversion  factor  from  ft  to  m 


Conversion  factor  from  hp  to  kW 


Conversion  factor  from  in3  to  m 


Conversion  factor  from  in  Hg  to  Pascal 


Conversion  factor  from  in  to  m 


Conversion  factor  from  kPa  to  Pa 


Conversion  factor  from  lb-mass  to  kg 


Conversion  factor  from  lb-force  to  N 


Conversion  factor  from  ft-lb  to  Nm 


Conversion  factor  from  lb-mass  ft2  to  kgm2 


Conversion  factor  from  lb  per  hp  hr  to  kg  per  kWhr 


Conversion  factor  from  liter  to  m 


Conversion  factor  from  psi  to  Pa 


Conversion  factor  from  rpm  to  rad  per  sec 


Units  |  Type 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Table  16.  Engine-Generator  Coupling 


Variable  |  Description  Units  |  Type 


cpl1_speed_ratio  |  Coupling  ratio  |  Scalar 


Variable 


cycle 


cycle  name 


road_grade 


_ Table  17.  Driving  Profile 


_ _ Description 


Table  of  time  (first  column  in  sec)  versus  desired 
speed  (second  column  in  m/s) 


Name  of  driving  cycle 


Gradient  of  road 


Units 


Vector  (Nx2) 


Scalar 


Scalar 
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Variable 


dvr1_weight 


dvr_deriv_gain 


dvr_integral_gain 


dvrproportjain 


Table  18.  Driver 


Description 


Name  of  driver 


Driver  weight 


Derivative  gain  of  driver  PID  controller 


Integral  gain  of  driver  PID  controller 


Proportional  gain  of  driver  PID  controller 


Units  |  Type 


Scalar 


Scalar 


Scalar 


Scalar 


Scalar 


Variable 


eng  1_cat  cnvtr  name 


eng  1  _cat_cnvtr_weight 


eng1_fuel_HV 


eng1_fuel_density 


eng1_fuel_weight 


eng1_name 


eng1_weight 


eng_af_ratio_Pi 


eng_af_ratio_speed 


eng_af_ratio_table 


eng_ambient_air_density 


eng_ambient_pressure 


eng_ambient  temperature 


eng_catalyst_eff_CO 


eng_catalyst_eff_HC 


eng_catalyst_eff_NO 


eng_catalyst_eff_PM 


eng_compression  ratio 


eng_displacement 


eng_emissions_CO_table 


eng_emissions_HC_table 


eng_emissions_NOx_table 


eng_emissions_PM_table 


eng_emissions_Pi 


eng_emissions_speed 


eng_exh_pres_scaling 


eng_exh_temp_Pi 


eng_exh_temp  af  coeff 


Table  19.  Spark  Ignition  Engine 


Description 


Engine  catalic  converter 


Engine  catalic  converter  weight 


Engine  fuel  heat  value 


Engine  fuel  density 


Engine  fuel  weight 


Engine  name 


Engine  weight 


Vector  of  intake  manifold  pressure 
for  air  fuel  ratio 


Vector  of  engine  speed  for  air  fuel 
ratio 


Engine  air  fuel  ratio  table 


Engine  ambient  air  density 


Engine  ambient  pressure 


Engine  ambient  temperature 


Catalyst  conversion  efficiency  for 
carbon  monoxide 


Catalyst  conversion  efficiency  for 
hydro  carbons 


Catalyst  conversion  efficiency  for 
nitrogen  oxides 


Catalyst  conversion  efficiency  for 
particulate  matter 


Engine  compression  ratio 


Engine  displacement 


Table  of  engine  emissions  for  carbon 
monoxide 


Table  of  engine  emissions  for  hydro 
carbons 


Table  of  engine  emissions  for 
nitrogen  oxides 


Table  of  engine  emissions  for 
particulate  matter 


Vector  of  intake  manifold  pressure 


Vector  of  engine  speed  for  different 
emissions 


Scaling  factor  for  matching  exhaust 
pressure 


Vector  of  intake  manifold  pressure 
for  exhaust  temperature 


Engine  exhaust  temperature  to  air 


Units 


Type 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


vector  (1x3) 


rad/sec  I  vector  (1x2) 


Matrix  (2x3) 


scalar 


scalar 


scalar 


scalar 


kg/kg 

fuel 


kg/kg 

fuel 


kg/kg 

fuel 


kg/kg 

fuel 


Pa 


rad/sec 


scalar 


scalar 


scalar 


Matrix  (5x5) 


Matrix  (5x5) 


Matrix  (5x5) 


Matrix  (5x5) 


Vector  (1x5) 


Vector  (1x5) 
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Variable 


eng_exh_temp_speed 


engexhjemp  table 


eng_fmep_factor_0 


eng_fmep_factor_1 


eng_fmep_factor_2 


eng_fuel_evaporation_time_co 

nstant 


eng_gas_const  air 


eng_gas_const_exhaust 


eng  inertia 


eng_liquid_fuel_fraction 


eng  manifold  volume 


eng_max_air  flow 


eng_mol_wt_air 


eng_mol_wt_exhaust 


eng_muffler_pressure_factor 


eng_nominal_isfc 


eng_off_torque_fact 


eng_rated_air  mass  flow 


eng  rated  air  vol  flow 


eng_rated_bmep 


engratedexhaust_density 


engrated  exhaust_pressure 


eng_rated_exhaust_temperatu 

re 


eng_rated_intake_vacuum 


engrated  power 


eng_rated_pres_ratioJnfl 


eng_rated_speed 


eng_rated_torque 


eng_rated  vol  eff 


eng_speed_density_weight 


engstdairdensity 


Table  19.  Spark  Ignition  Engine 


_ Description 


fuel  ratio  coefficient 


Vector  of  engine  speed  for  exhaust 
temperature 


Exhaust  temperature  table 


Correlation  factor  0  for  friction  mean 
effective  pressure 


Correlation  factor  1  for  friction  mean 
effective  pressure 


Correlation  factor  2  for  friction  mean 
effective  pressure 


Time  constant  for  fuel  evaporation  in 
intake  manifold  wall 


Specific  gas  constant  for  air 


Specific  gas  constant  for  exhaust 
gases 


Engine  rotary  inertia 


Fraction  of  fuel  that  stays  in  liquid 
form 


Engine  intake  manifold  volume 


Engine  max  air  flow 


Molecular  weight  of  air  through 
throttle 


Molecular  weight  of  exhaust  gases 


Engine  muffler  pressure  correlation 
factor 


Engine  nominal  indicated  specific 
fuel  consumption 


Factor  for  torque  when  engine  is  off 


Rated  air  mass  flow  for  engine 


Rated  air  volume  flow  for  engine 


Rated  Brake  Mean  Effective 
Pressure 


Rated  exhaust  air  density 


Rated  exhaust  gas  pressure 


Rated  exhaust  gas  temperature 


Intake  manifold  vacuum  at  rated 
conditions 


Rated  engine  power 


Pressure  ratio  influence  factor  (in 
flow  equation)  at  rated  conditions 
[fraction] 


Maximum  engine  speed 


Rated  engine  torque 


Rated  engine  volumetric  efficiency 


Fractional  weight  for  speed  density 
air  flow  compared  to  mass  flow 
sensor  [fraction] 


Density  of  standard  air 


Units 


rad/sec  I  vector  (1x2) 


Matrix  (2x2) 


scalar 


J/kg/K 


J/kg/K 


kg/s 


kg/kmol 


kg/kmol 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


Table  19.  Spark  ignition  Engine 


Variable 

Description 

eng_std_pres 

Standard  pressure 

eng_std_temp 

Standard  temperature 

eng_stoich_af_ratio 

Correlation  factor  for  effect  of  air  fuel 
ratio  on  exhaust  temperature 

eng_throttle_bypass_ratio 

Ratio  of  bypass  flow  area  around 
throttle  to  total  throttle  area 

eng_universal_gas_const 

Universal  gas  constant 

eng_vol_eff_Pi 

Table  of  intake  manifold  pressure 
vs.  vol  efficiency 

eng_vol_eff_speed 

Table  of  engine  speed  for  vol 
efficiency 

eng_vol_eff_table 

Table  of  engine  volumetric 
efficiencies 

Units 


Pa 


K 


Type 


scalar 


scalar 


scalar 


J/kmol/K 


Pa 


scalar 


Vector  (1x2) 


Vector  (1x6) 


Matrix  (6x2) 


Variable 


essb_chrg_eff 


essb_dischrg_eff 


mss! 


_ Table  20.  Lead  Acid  Battery 


Description 


Battery  name 


Battery  weight 


Internal  capacitance  of  a  cell 


Discharge  capacity  @  5  hour  rate  at  30  Celsius  for  a  cell 


Terminal  capacitance 


Open  circuit  voltage  of  a  cell 


Discharge  current  @  5  hour  rate  for  a  cell 


Polarization  resistance  of  a  cell 


Initial  charge  in  battery  cell 


Terminal  resistance  for  the  battery 


Battery  charge  efficiency 


Battery  discharge  efficiency 


Battery  initial  polarization  voltage 


Battery  initial  state  of  charge 


Internal  resistance  of  a  cell 


Total  number  of  batteries 


Units 


Type 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


|  Table  21.  Permanent  Magnet  Generator 

Variable 

Description 

genp1_ctrl_dgain 

Derivative  gain  of  PID  controller  for  generator 

genp1_ctrl_igain 

Integral  gain  of  PID  controller  for  generator 

genp1_ctrl_pgain 

Proportional  gain  of  PID  controller  for  generator 

genp1_ctrlr_name 

Generator  controller  name 

genp1_ctrlr_weight 

Generator  controller  weight 

genp1_name 

Generator  name 

genp1_weight 

Generator  weight 

Units  |  Type 


scalar 


scalar 


scalar 


scalar 


N  |  scalar 


scalar 


N  scalar 
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Variable 


motll  ctrlr  name 


motllctrlr  weight 


motll  name 


mot!1_weight 


motor  N  max 


Table  22.  AC  Induction  Motor 


_ Description 


Motor  controller  name 


Motor  controller  weight 


Motor  name 


Motor  weight 


Maximum  rated  speed  of  motor 


Units 


scalar 


scalar 


scalar 


N  scalar 


rpm  scalar 


_  Table  23.  Simulation  Constants 


Variable _ Description 


sample_N  Number  of  samples  to  save  to  workspace 
sample_T  Time  step  length  between  saved  samples 


Units  |  Type 


scalar 


scalar 


Variable 


_aero_drag  coef 


air  density 


_chassis_weight 


_coef_yaw_sens 


fO 


front  area 


fs 


mass 


name 


normal  wind 


_parallel_wind 


=passengers  weight 


_roadtire  frict  coeff 


Table  24.  Vehicles 


_ Description 


Aerodynamic  drag  coefficient 


Air  density 


Chassis  Weight 


Coefficient  of  sensitivity  to  yaw 


Vehicle  speed  independent  rolling  coefficient 


Frontal  area  of  vehicle 


Vehicle  speed  dependent  rolling  coefficient 


Vehicle  mass 


Vehicle  name 


Wind  velocity  normal  to  vehicle  motion 


Head  wind  velocity 


Passenger  weight 


Road  tire  coefficient  of  friction 


Units  |  Type 


scalar 


kg/m3  scalar 


N  |  scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


scalar 


Variable 


weight  total 


Table  25.  Miscellaneous 


Description 

Units 

Total  vehicle  weight 

N 

scalar 


Table  26.  Wheel 


Variable  | 

Description 

Units 

Type 

Wheel  mass 

kg 

scalar 

lyim  ■  frMffPWB 

Wheel  radius 

m 

scalar 

E23KE9HHM 

Wheel  radius 

m 

scalar 
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HYBRID  ELECTRIC  VEHICLE 
SIMULINK  TOOLBOX 

(Commercial  Vehicle  Version) 

Technical  Notes 

Version  1.1 
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1.0  INTRODUCTION 


These  technical  notes  provide  the  equations  used  within  the  individual  modules  of  the  DARPA 
Hybrid  Electric  Vehicle  (HEV)  simulation  toolbox,  hereinafter  referred  to  as  PATHS 
(Performance  Assessment  Toolbox  for  Hybrid  Systems).  These  technical  notes  accompany 
the  User’s  Guide  for  the  HEV  simulation  toolbox.  Use  of  PATHS  is  restricted  to 
government-funded  projects,  following  explicit  permission  by  DARPA/SwRI  only.  All 
inquiries  to  obtain  permission  to  use  PATHS  must  be  directed  to  Dr.  A.  Nedungadi  at 
Southwest  Research  Institute  (210-522-3965). 


PATHS  is  developed  under  the  MATLAB/SIMULINK®  environment,  and  is  provided  to 
users  in  source  code  form.  An  experienced  MATLAB/SIMULINK®  user  can  easily  trace  the 
inputs  to  the  outputs  and  determine  the  transfer  function  between  these  two.  The  following 
notes  provide  the  user  with  the  basic  principles  of  physics  that  have  been  employed  to  model 
the  individual  modules  of  PATHS.  Modules  that  use  table  lookup  methods  are  not  described 
in  this  document,  since  these  are  not  based  on  principles  of  physics.  Cursory  equations 
describing  the  internal  operation  within  certain  modules  are  provided  in  the  User’s  Guide.  In 
this  case,  the  user  is  referred  to  the  User’s  Guide  for  information  on  the  underlying  technical 
equations.  All  equations  of  controller  modules  of  the  controllers  library  (refer  to  User’s 
Guide)  have  been  omitted  since  these  are  based  on  classical  PID  control  laws.  As  a  result, 
this  document  contains  technical  equations  for  the  following  modules,  which  is  a  subset  of 
the  entire  PATHS  software  (please  refer  to  the  User’s  Guide  for  a  complete  listing  of  all 
modules  of  PATHS): 

ELECTRIC  MACHINE  LIBRARY 

I  I  AC  induction  steady  state  motor 
I  I  Permanent  magnet  steady  state  generator 
I  I  AC  Induction  dynamic  motor 
I  I  Permanent  magnet  dynamic  generator 
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ENERGY  STORAGE  LIBRARY 

|  |  Lead  acid  battery 
|  |  flywheel  battery 
I  |  Capacitor 

VEHICLE  LIBRARY 

|  |  Wheeled  vehicle  (version  1  and  2) 

|  |  Wheel  (version  1  and  2) 

COUPLINGS  LIBRARY 

n  Gear  box 

|  |  Engine-generator  coupling 
[  |  Motor-wheel  coupling  (version  1  and  2) 

MISCELLANEOUS  LIBRARY 
I  |  Roller  bearing 
|  |  Magnetic  bearing 
I  |  Windage  model 
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2.0  ELECTRIC  MACHINE  LIBRARY 


2.1  Steady  State  AC  Induction  Motor 


Figure  32  shows  the  inputs  and  outputs  of  the  steady-state  induction  motor  model  - 
motll  module.  Figure  33  depicts  the  per-phase  equivalent  circuit  model  of  the  induction 
machine  and  defines  the  motor  parameters.  Figure  34  shows  the  induction  motor  model  after 
the  transformation  to  the  synchronous  d-q  reference  frame.  The  electrical  terminal  voltages 
and  currents  are  calculated  directly  from  this  equivalent  circuit  and  the  mechanical  equations 
using  the  command  values  of  torque  and  speed.  The  synchronous  reference  frame  electrical 
and  mechanical  equations  used  in  the  simulation  are, 

Vqs  =  K'ls  + 

Vds=RJds-0Xqs 

0  =  Rfdr-(a>-cor)Xqr 
K,  =  (Am  +  ArHfc  +  LJdr 

^=(4,+4)£+4& 

^=(4  +  4)C  +  4C 
r  =  A»  (&&“&»,%) 

where  the  e  superscript  denotes  the  synchronous  reference  frame. 


DC  Inverter 
Voltage  [V] 


„  I 


_r 


& 


Motor  Rotor  Speed  -) 
[fpm]  _ ^ 


DC  Inverter 
Current  [A] 


motl1_module 


Per  Unit  Torque 
Command 


Motor  Torque 
[N*m] 

Motor  Inertia 
[kg*mA2] 


->0 


Power  Vector  [Input  Power,  Output  Power, 
lnv_Loss  Power,  MachJ_oss  Power] 
[W] 


Figure  32.  Inputs  and  Outputs  of  motll  module 
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The  motor  is  controlled  using  conventional  rotor  flux  field  orientation  (i.e.,  “brushless  dc” 
control)  such  that  in  the  synchronous  reference  frame,  \{e=  0. 
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2.2  Steady  State  Permanent  Magnet  Generator 


Figure  35  shows  the  inputs  and  outputs  of  the  steady-state  generator  model  -  genpl_module. 
Figure  36  depicts  the  steady-state  synchronous  d-q  reference  frame  model  used  in  the 
simulation.  The  electrical  terminal  voltage  is  calculated  directly  from  this  equivalent  circuit 
using  the  command  values  of  current  and  speed.  The  electrical  and  mechanical  equations 
used  in  the  simulation,  and  given  by  the  equivalent  circuit,  are, 


4  iqm  r  ^r^ds^dm  ^md^df ) 


4  = 


ii- 


R 


V  44  ^r^md^df  ^r^ds^dm 

VI  =  -RA  +  <orLqfqm 

^  ~  ^md^df^qm 


where  the  superscript  e  denotes  the  synchronous  reference  frame,  and  Idf  is  the  fictitious 
current,  which  corresponds  to  the  strength  of  the  permanent  magnet.  This  current  is 
internally  calculated  from  the  no  load  speed  and  voltage.  The  iron  loss  resistor,  R,,,  ,  is 
calculated  from  the  rated  core  losses  value  provided  by  the  user. 


Figure  35.  Inputs  and  Outputs  of  genpl  module 
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Figure  36.  Steady-State  dq-Equivalent  Circuits 


The  generator  is  controlled  using  conventional  synchronous  reference  frame  field  orientation 
(i.e.,  “brushless  dc”  control)  such  that  i^  =  0. 


2.3  AC  Induction  Dynamic  Motor 


Figure  37  below  shows  the  inputs  and  outputs  for  the  induction  motor/boost  inverter  system. 
The  standard  steady-state  induction  machine,  per-phase  equivalent  circuit  is  shown  in  the 
previous  section.  Please  refer  to  this  section  for  information  about  what  each  resistance  or 
inductance  referenced  in  the  mask  relates  to  electrically. 


Figure  37.  Inputs  and  Outputs  of  motkl  module 
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Figure  38  depicts  the  synchronous  reference  frame  d-q  dynamic  model  of  the  induction 
machine  and  defines  the  motor  parameters.  The  electrical  terminal  voltages  and  currents  are 
calculated  directly  from  this  equivalent  circuit  and  the  mechanical  equations  using  the 
command  values  of  torque  and  speed.  The  synchronous  reference  frame  electrical  and 
mechanical  equations  used  in  the  simulation  are, 

Vjs=Rjeds-aXqs+%s 
Q  =  Rfqr+(<o-o>r)Xdr+$qr 
Q  =  Rfdr-(.<0-0)r)Xqr+%dr 

4=(4+4)C+4v 
=  (4  +  4  XL  +  44 
^,=(4+4)4  +  44 
**=(4  +  4)4  +  44 
'  =  4  (Wr-Wr) 

where  the  e  superscript  in  the  equations  denotes  the  synchronous  reference  frame. 

The  motor  is  controlled  using  conventional  rotor  flux  field  orientation  (i.e.,  “brushless  dc” 
control)  such  that  in  the  synchronous  reference  frame,  A.qre  =  0-  The  stator  currents,  iqse  and 
i^,  are  controlled  in  the  synchronous  reference  frame  through  a  standard  PI  current 
regulator.  The  PI  gains  are  user-definable. 

.  4  t-'lr  (to-1 C0r)A.drC 


R,’ 


.  rs  ®V  Lb  4  (CB-CDr)4 


r; 


Figure  38.  Dynamic  dq  equivalent  circuit  of  an  induction  motor 


103 


2.4  Permanent  Magnet  Dynamic 


Figure  39  shows  the  inputs  and  outputs  of  the  dynamic  generator  model  -  genql_module. 
Figure  40  depicts  the  dynamic  synchronous  d-q  reference  frame  model  used  in  the 
simulation.  The  electrical  terminal  voltage  is  calculated  directly  from  this  equivalent  circuit 
using  the  command  values  of  current  and  speed.  The  electrical  and  mechanical  equations 
used  in  the  simulation,  and  given  by  the  equivalent  circuit,  are, 

4  =  4  +  KV1  -  o)rLmdIdf  +  Ljt) 

Km 

4  =  4-7-  KV4  -  AJS.) 

M 

9 

^ qs  ~  ^s^qs  ^r^md^df  ^r^ds^dm  ^qs^qm 

9 

V*  =  -R&  +  0)rLq/qm  -  Ljfm 

^  Lmd^dfiqm 

where  the  superscript  e  denotes  the  synchronous  reference  frame,  and  1^  is  the  fictitious 
current,  which  corresponds  to  the  strength  of  the  permanent  magnet.  This  current  is 
internally  calculated  from  the  no  load  speed  and  voltage.  The  iron  loss  resistor,  R^, ,  is 
calculated  from  the  rated  core  losses  value  provided  by  the  user. 


Figure  39.  Inputs  and  Outputs  of  genql  module 
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Figure  40.  Dynamic  Generator  dq-Equivalent  Circuits 


The  generator  is  controlled  using  conventional  synchronous  reference  frame  field  orientation 
(i.e.,  “brushless  dc”  control)  such  that  i^  =  0.  The  q-axis  current,  iqm,  is  controlled  in  the 
synchronous  reference  frame  through  a  standard  PI  current  regulator.  The  PI  gains  are  user- 
definable. 


3.0  ENERGY  STORAGE  LIBRARY 


3.1  Lead  Acid  Battery 

Figure  41  shows  the  circuit  diagram  used  to  develop  the  dynamic  model  of  the  lead-acid 
battery.  A  single  battery  cell  is  described  by  the  following  parameters: 


Internal  resistance  of  cell  i 

E,M: 

Open  circuit  voltage  of  cell  i 

ci  [F]: 

Capacitance  of  cell  i 

Ki[Q]: 

Polarization  Resistance  of  cell  i 

SOQ: 

State  of  charge  of  cell  I 

r. 


Figure  41.  Circuit  Diagram  for  a  Single  Cell  of  a  Lead  Acid 
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In  a  battery  with  multiple  cells,  the  circuit  will  consist  of  N  single  cells  in  series  with  a  terminal 
resistance  Rx  [Q]  and  capacitance  CT  [F].  Figure  42  shows  2  single  cells  in  series,  a  terminal 
resistance  Rx  [Q]  and  capacitance  Cx  [F].  Positive  current  I  is  shown  by  the  direction  of  the 
arrow  in  Figure  41.  Referring  to  Figure  41,  and  assuming  that  c,  =  c2,  E,  =  E2  =  E,  K,  =  K2, 
KirchofF s  voltage  law  can  be  formulated  as  follows: 

N(Vc)  +  N(E)  +  N(Irr)  +  IrRx  =  VT  [1] 

where: 

Vc:  The  voltage  across  the  capacitor  c,  [V] 

N :  The  number  of  cells  in  the  battery  [  ] 

Ir:  The  portion  of  the  current  I  that  flows  into  the  internal  resistance  r  [A] 

Vx:  The  terminal  voltage  [V] 


r2 

-AAAAAA- 


E, 


K2/< 


soc, 


r, 

-AAAMA- 


^Vc,h 


Kj/( 


1 

|C2 

LPi 

1 

1 

1 

-VWW- 

-vwv\J 

-\aaaaa- 

Rt 


soc, 


Figure  42.  Circuit  Diagram  for  a  Lead  Acid  Battery  with  Two  Cells 


The  above  lead  acid  battery  model  has  two  states:  Vc  and  Vx.  Therefore,  this  is  a  second  order 
battery  model.  Assuming,  SOC,  =  SOC2  =  SOC,  the  above  equation  can  be  combined  with 
KirchofF s  current  law  to  determine  these  two  states  as  follows: 
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[2] 


I  =  C. 


dVT 

dt 


+  1. 


j  dVc  Vc 
It  =  c - +  ■ 


dt 


[3] 


/soc 

Solving  for  Ir  in  equation  (1 1)  and  substituting  it  in  equations  (12)  and  (13)  yield  the  following 
for  the  two  states  Vc  and  VT 


dVT 

dt 


I- 


(Vt-N(E)-N(Vc) 
N(r)  +Rt 


<  l' 


vcT; 


dVc  _ 

"(Vt-N(E)-N(Vc)' 

Vc 

dt 

N(r)  +  Rt 

K/ 

/socj 

[4] 


[5] 


Assuming  positive  current  (charging)  is  in  the  direction  shown  in  Figure  41,  the  state  of  charge 
SOC  of  the  battery  can  be  calculated  as  a  function  of  the  battery  capacity.  The  battery  capacity 
can  be  expressed  as  a  function  of  temperature  and  discharge  current  Id  as  follows: 


C(I,T)  =  C, 


/  \  -0.325 


UJ 


[l  +  0.008(T  -  30)] 


where: 

C5:  Internal  battery  capacity  [A-h] 

I5:  Discharge  current  at  5  hour  rate  for  the  battery 

Id:  Discharge  current 

T:  Temperature  at  which  the  battery  operates  [Deg  C] 


[6] 


During  the  discharge, (Id  <  0),  the  SOC  can  be  calculated  as  follows: 
dSOC  100 


dt  C(Id,T) 


-Id 


m 
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where: 


t|d:  Discharge  efficiency  [  ] 


During  charge  (Ic  >  0),  the  SOC  can  be  calculated  as  follows: 

dSOC  _  (100-SQC)ticIc 
dt  Q 

where: 


[8] 


lie 

Ic 

Q 


Charge  efficiency  [  ] 

Charge  current  [A] 

Total  Quantity  missing  in  the  battery  over  the  entire  operating  time  t 


Q  is  calculated  as  follows: 


Q  = 


Jl* 


[9] 


3.2  Flywheel  Battery 

There  are  numerous  possible  flywheel  energy  storage  system  configurations.  One  of  these 
configurations  is  shown  in  Figure  43. 
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Units 


Voltage: 

speed: 

current: 


Volts 


RPM 


Notes 

inputs: 

states: 

outputs: 


Amps 

2(voltage,  current  command) 
l+(speed  and  generator  states) 

4(  output  current,  speed,  discharge  &  charge  eff) 


direct  feedthrough :  yes 


There  are  two  major  components  and  one  auxiliary  component  of  this  system: 

1 .  Motor/Generator  and  Converter 

2.  Flywheel 

3.  Efficiency  Calculation 

This  subsystem  uses  a  permanent  magnet  generator.  The  'Current  Command'  is  the  current 
command  provided  to  the  generator  controller.  The  'Output  Current'  is  the  actual  current 
output  to  the  DC  Bus.  Inputs  for  the  'Flywheel  Subsystem’  are  the  DC  voltage,  generator 
torque,  generator  inertia  and  bearing  loads.  The  'Flywheel  Subsystem'  has  outputs  of 
auxiliary  current  draw,  stored  mechanical  power,  and  rotor  speed. 
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DC  Voltage 


Bearing 

Loading 


— ►  Magnetic  — 
^  Bearing  _ 

Radial  Bearing  1 


Magnetic 

Bearing 


Rect-Cyl 


Generator 

Inertia 


MuX>|u{2)/u(1)|_^|j-| 

Amps 

Draw 


Acceleration  Torque 
Jpi/3&\ 


ywnee 


PStore 

W21 


Power 

Stored 


Speed 


Figure  44.  Flywheel  Subsystem  with  Magnetic  Bearings 


DC  voltage 


Bearing 

Loading 


Rect-Cyl 


Generator 

Torque 


Generator 

Inertia 


Roller 

>  Bearing  — 

[~H 

Roller  Bearing  1 


Windage 


e 


Vacuum 

Pump 


Mux  Ju(2)/u(1)  [_►[!" 


Amps 

Draw 


Acceleration  Torque 


3 


ly  wheel 


n  PStore 

pS/3fr; 

Jt->2 

Power 

Stored 

Speed 

- JT1 

Figure  45.  Flywheel  Subsystem  with  Roller  Bearings 
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The  vacuum  system  is  assumed  to  have  a  constant  electrical  power  draw  specified  by  the 
user.  The  inputs  to  the  bearings  are  the  forces  on  the  bearings.  These  forces  will  be  highly 
dependent  on  the  mounting,  with  particular  attention  to  whether  gimbal  mounting  removes 
gyroscopic  effects.  The  torque,  which  must  be  reacted  by  the  bearings  as  radial  forces  to 
overcome  gyroscopic  effects,  follows  the  following  form: 

T  —  Jft)  rotor  0 

where 

t:  reaction  torque  to  rotate  the  flywheel  (N-m) 

J :  moment  of  inertia  of  the  rotor  (Kg-m) 

wrotor:  rotor  speed  (rad/sec) 

0 :  rate  of  change  of  rotor  position  (rad/sec) 

Because  of  the  fairly  high  rotor  speeds,  a  fairly  low  rate  of  change  of  rotor  angle  can  generate 
a  large  gyroscopic  torque.  If  a  small  bearing  span  is  used  this  could  generate  large  bearing 
forces.  The  equations  to  translate  vehicle  accelerations  from  the  six  degrees  of  freedom  (X, 
Y,  Z,  yaw,  pitch  and  roll)  into  bearing  forces  will  be  mounting  dependent  and  is  left  to  the 
system  integrator. 

3.2.1  Efficiency  Calculation  for  Flywheel  Energy  Storage  System 

Inputs  for  the  efficiency  calculation  subsystem  'Eff.  Calc.'  are  the  power  going  into  stored 
energy  and  the  power  going  into  electrical  energy.  The  efficiency  calculation  subsystem  uses 
the  following  equations. 


^discharge 


_  ^-out, discharge 
^in,discharge 


^charge 


~out, charge 
“in, charge 


There  is  a  case  where  electrical  power  flows  into  the  flywheel  energy  storage  system  but 
stored  energy  also  decreases.  With  the  above  equations,  these  losses  are  included  in  the 
efficiency  calculation  by  counting  them  as  charging  losses. 


Figure  46.  Flywheel  Energy  Storage  Efficiency  Calculation 


3.2.2  Flywheel  Battery  System 

A  flywheel  battery  system  can  be  defined  as  a  flywheel  energy  storage  system  with  an 
associated  filter  capacitor  and  controller  to  regulate  a  DC  bus  voltage.  The  relationship 
between  the  DC  voltage  and  input  current  could  be  similar  to  a  lead-acid  battery,  but  a 
control  scheme  that  maintains  a  fairly  constant  voltage  can  be  used.  A  flywheel  battery 
system  with  a  simple  feed-forward  control  and  low  gain  feedback  is  shown.  The  feed-forward 
signal  exactly  balances  the  input  current.  There  is  a  low  gain  feedback  to  keep  the  bus 
voltage  near  it's  reference  (300  volts  in  this  case)  to  adjust  for  losses. 
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Input 

Current 


Filter  Cap. 


Figure  47.  Simple  Flywheel  Battery  System 


3.3  Capacitor 


Mathematical  Model 

V  dV 

I-  — =  C — 

R  dt 

Units 

current: 

Amps 

voltage: 

Volts 

resistance: 

Ohms 

capacitance: 

Farads 

Notes 

inputs: 

1  (current) 

states: 

1  (Voltage) 

outputs: 

1  (Voltage) 

direct  feedthrough: 

no 
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4.0  VEHICLE  LIBRARY 


4.1  Wheeled  Vehicle  (Version  1) 

This  wheeled  vehicle  model  calculates  all  the  load  forces  acting  on  the  vehicle.  The  total 
load  force  is  as  follows: 

•^Load  ^*aero  F  ^roll  F  ^grade 

where: 

Fioad  [N];  Total  resistance  force  acting  on  the  vehicle 

Faero  [N]:  Aerodynamic  load  on  vehicle 

F roii  [N]:  Rolling  resistance  from  the  tires  on  the  vehicle 

Fgrade  [N] :  Resistance  as  a  result  of  the  slope  of  the  road  (road  grade  force) 

The  aerodynamic  drag  force  Faero  is  expressed  as  a  function  of  the  vehicle  frontal  area  (A), 
aerodynamic  drag  coefficient  (Q),  the  air  density  (r),  the  coefficient  of  yaw  sensitivity  (a),  the 
normal  (V„)  and  parallel  (Vp)  wind  velocity  component,  and  the  vehicle  velocity  (V)  as  follows: 

Faen,  =  0.5  pCd  A(V -V  P  f  (1  +  a  tm'  (-w~— ))  1 

('  V p) 


where: 

Fae™[N] : 

Aerodynamic  drag  force 

p  [kg/m3]: 

Air  density 

Cd  [  ]: 

Aerodynamic  drag  coefficient 

A  [m2]: 

Vehicle  frontal  area 

V  [m/s]: 

Vehicle  velocity 

V„  [m/s]: 

Normal  wind  velocity  component 

Vp  [m/s]: 

Parallel  wind  velocity  component 

a[]: 

Coefficient  of  yaw  sensitivity 
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The  rolling  resistance  force  is  expressed  as  follows: 


Fn,ii  =  //Weight  cos(r) 


where: 


V  [mph]: 

m  [  ] : 
Weight  [N]: 

fo[]: 

f,[J: 


Vehicle  speed 
Coefficient  of  friction 
Total  vehicle  weight 
Coefficient  of  rolling  resistance  [  ] 
Speed  dependent  coefficient  of  friction 


The  road  grade  force  is  expressed  as  follows: 

Fgrade  =  Weight  sin(a) 


(2) 


2 


where: 

Weight  [N]:  Total  vehicle  weight 

a  [rad]:  Road  gradient 


4.2  Wheeled  Vehicle  (Version  2) 

In  this  model,  the  vehicle  is  modeled  as  a  single  mass,  whose  acceleration  is  computed  using 
Newton’s  second  law  of  motion: 

^  eq  ^  ^driving  ^load 

where: 

meq  [kg]:  Equivalent  mass  of  vehicle 


116 


x  [rad/sec2]:  Vehicle  acceleration 

^driving  [N]'  Total  driving  force  on  vehicle 

Fload  [N] :  Total  resistance  load  on  vehicle 

The  mass  of  vehicle  takes  fuel  consumed  mass  (as  fuel  is  consumed,  total  fuel  mass  decreases) 
into  consideration.  The  vehicle  acceleration  is  integrated  once  to  determine  vehicle  speed,  and 
again  to  determine  vehicle  position.  The  computation  of  Fload  is  the  same  as  described  in  the 
wheeled  vehicle  version  (1),  with  the  exception  of  including  Fbrake,  which  is  provided  by  the 
driver  model.  Therefore,  Fload  is  computed  as  follows: 

^Load  Faero  Froj|  Fgrafje  +  F|Jrajce 


4.3  Wheel  Model  (Version  1) 

Figure  48  shows  the  physical  model  of  the  vehicle  and  the  driven  wheel,  along  with  all  forces 
acting  on  the  vehicle.  The  following  is  an  explanation  of  all  forces  acting  on  the  vehicle: 

Fioad  [N]:  Sum  of  all  load  forces  opposing  the  motion  of  the  vehicle.  This 

force  comprises  the  Aero  dynamic  load  and  the  gradient  load. 

Ff  [N] :  The  friction  force  acting  on  the  surface  of  the  tires  as  a  result  of  the 

friction  coefficient  between  the  tires  and  the  road, 
m  [Kg] :  The  total  mass  of  the  vehicle  (including  the  wheel  mass) 

x  [m] :  The  linear  displacement  of  the  vehicle 
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X 


Figure  48.  Physical  Model  of  the  Vehicle  and  Wheel 


For  the  system  shown  in  Figure  48,  Newton’s  law  can  be  formulated  as  follows: 


mx  =  Ff  -  Fload  (1) 


Figure  49  shows  the  physical  model  of  the  wheel.  The  following  is  an  explanation  of  the  forces 
and  moments  acting  on  the  wheel: 


Td  [Nm]: 
Ff  [N] : 


Jw[kgm2]: 
q  [rad]: 


The  driving  torque  on  the  wheel  delivered  from  the  axle 

The  friction  force  acting  on  the  surface  of  the  tires  as  a  result  of  the 

friction  between  the  tires  and  the  road. 

The  mass  moment  of  inertia  of  the  wheel 
The  rotation  in  radians  of  the  wheel 
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Assuming  a  constant  tire  radius,  Newton’s  law  can  be  formulated  for  the  wheel  as  follows: 

Jw9  =  Td-Ffrt  (2) 

From  Equation  (1),  Ff  can  be  expressed  as  follows: 

Ff  =  mx  +  Fload  (3) 

Combining  equations  (2)  and  (3),  it  follows: 

JWe  =  Td-mxr,-Floadr, 
x  =  r,0 

Jw0  =  Td-mr?0-Floadrt  (4) 

(Jw  +  mr?)0  =  Td  -  Fload  rt 

0  =  fioad  gt 

Jw  +  mr? 
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4.4  Wheel  Model  (Version  2) 


This  wheel  model  converts  the  vehicle  velocity  into  the  rotational  speed  of  the  wheel,  which  is 
equal  to  the  rotational  speed  of  the  output  shaft  of  the  axle.  The  axle  or  coupler  connects  the 
electric  motor  to  the  wheel. 

03  wheel  =  V  /  rwhee! 

where: 

“wheel  [rad/sec]:  Rotational  speed  of  the  wheel 
rwheei  [ml:  Wheel  radius 

V  [m/s]:  Vehicle  linear  velocity 

The  driving  torque  acting  on  the  wheel  is  converted  to  driving  force  as  follows: 


'  —  1  axle  _  out 

driving 

^ wheel 


where: 

^driving  [N]: 

Total  driving  force  acting  on  vehicle 

'faxle  out  [Nm]. 

Total  output  torque  delivered  from  axle 

Tvhee!  [^]* 

Wheel  radius 

If  the  entire  mass  of  the  vehicle  is  regarded  as  concentrated  at  the  wheel  mass  center,  then  the 
total  inertia  of  the  wheel  can  be  computed  as  follows: 

m^  *  r^  =  (1  /  2  *  mwheel  *  rtthed  2  +  Jaxle_  out 

where: 

meq  [Kg]:  Equivalent  vehicle  mass 

rwheel[m] :  Wheel  radius 

Jaxle  outpt  [kgm2] :  Output  Inertia  of  axle 
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5.0  COUPLINGS  LIBRARY 


5.1  Motor  Wheel  Coupling  (Version  1) 

Figure  50  shows  the  physical  model  of  the  motor  wheel  coupling  as  well  as  the  free  body 
diagram  for  each  gear.  The  following  is  an  explanation  of  all  forces  and  moments  shown  in 
Figure  50: 


T,  [Nm]: 

Torque  at  the  input  side  of  the  coupler  model 

T2  [Nm]: 

Torque  at  the  output  side  of  the  coupler  model 

F  [N]: 

The  contact  force  between  gear  1  and  gear  2 

q,  [rad]: 

The  rotation  of  the  shaft  at  the  input  side  of  the  coupler  model 

<h  M]: 

The  rotation  of  the  shaft  at  the  output  side  of  the  coupler  model 

J,  [kgm2]: 

Inertia  at  the  input  side  of  the  coupler  model 

J2  [kgm2]: 

Inertia  at  the  output  side  of  the  coupler  model 

ri  [m]: 

Radius  of  input  gear  of  the  coupler  model 

r2  [m]: 

Radius  of  output  gear  of  the  coupler  model 

G  []: 

Gear  ratio  of  the  coupler 

nD: 

Coupler  efficiency 

Figure  50.  The  Coupler  Physical  Model 
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Newton’s  law  can  be  formulated  for  both  gears  as  follows: 


9i 


T.-Fn 

Ji 


02 


Fr2  -  T2 
J2 


(h  =  =  II  =  (} 

02  02  ri 


(1) 


.  F  =  (J202  +  T2) 
r2 


0i  =  02G 


From  equation  (1),  the  coupler  model  can  be  formulated  as  follows  (traction): 

£2GJl-Tl  =  (J2&  +  T2) 

Cj 

T  2  =  [Ti-(GJi  +  (^-)J2^2)]G7 
Cj 

The  output  torque  T2  of  the  coupler  in  regen  mode  is  then: 
T2=[T,-(OJ,+(i)J2S2]G(-t) 

Lr  77 
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5.2  Motor  Wheel  Coupling  (Version  2) 

This  model  computes  the  total  inertia  as  seen  by  the  wheel  as  follows: 


(Ji„+Ji„gear)G2+J 


outgear 


=  J 


wheel 


where: 


4  [kgm2]: 

Jjngear  [kgm2]: 

G  []: 


Inertia  of  the  rotor  of  the  motor 

Inertia  of  the  input  side  of  the  gear  coupler  between  the  motor  and  wheel 
Gear  ratio  of  coupler 
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Joutgear  [kgm2]:  Inertia  of  the  output  side  of  the  gear  coupler  between  the  motor  and  wheel 
JWheei  [kgm2]:  Total  inertia  that  the  wheel  has  to  accelerate 

The  gear  ratio  is  defined  as  follows: 


where: 

Ni,  [rpm]:  Rotational  speed  at  the  input  side  of  the  coupler 

NoJrpm] :  Rotational  speed  at  the  output  side  of  the  coupler 

The  torque  output  of  the  coupler  is  computed  as  follows  (traction  mode): 
To*  =  TinG  rj 

The  torque  output  of  the  coupler  is  computed  as  follows  (regen  mode): 

T.  =  T,„G(i) 

V 


where: 

Tout  [Nm]: 

Torque  to  drive  wheel 

Tin  [Nm]: 

Torque  inputted  into  coupler  from  motor 

G  []: 

Gear  ratio 

T1  []: 

Coupler  efficiency 
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6.0  MISCELLANEOUS  LIBRARY 


The  equations  used  in  this  library  are  briefly  discussed  in  the  User’s  manual.  Therefore,  the 
reader  is  referred  to  the  User’s  Manual  for  the  equations  describing  the  models  of  the 
miscellaneous  library. 
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EXECUTIVE  SUMMARY 


The  following  report  summarizes  the  efforts  to  validate  PATHS,  a  modeling  and  simulation 
software  for  electric  and  hybrid  vehicles.  Due  to  time  and  budget  constraints,  it  was  impossible 
to  validate  all  the  components  of  the  PATHS  software.  Therefore,  only  key  components  of 
PATHS  were  validated.  The  APU,  AC  induction  motor,  lead  acid  and  driver-axle-wheel-vehicle 
combined  models  were  validated.  These  are  the  components  that  are  principally  used  to  model 
a  commercial  hybrid  vehicle. 

Validation  was  performed  by  first  developing  a  computer  model  of  a  22-foot  series  hybrid  shuttle 
bus,  followed  by  comparing  its  output  to  data  collected  from  the  same  bus  on  a  dynamometer. 
Section  1  describes  the  series  hybrid  bus  and  the  driving  profiles  used  during  validation. 

Section  1 . 1  describes  the  results  of  the  APU  model  validation.  Fuel  consumption  and  CO 
(Carbon  Monoxide)  and  NOx  (Nitrogen  Oxide)  emissions  were  measured  and  compared  with 
model  predictions.  Model  and  measured  data  agreed  with  each  other  to  within  less  than  5%, 
while  considerably  larger  deviations  up  to  100%  were  observed  between  measured  and  model 
predicted  CO  and  NOx  emissions.  This  can  be  largely  attributed  to  insufficient  emissions  data 
obtained  from  the  manufacturer. 

Section  1 .2  describes  the  validation  of  the  AC  induction  motor  model.  The  consumed  motor 
current  in  motor  mode  and  regenerated  current  in  generator  mode  were  compared  with  data 
collected  from  the  motors  of  the  bus.  The  model  predicts  transient  and  steady-state  operations 
of  the  motor  within  10%,  except  for  measured  motor  current  spikes  in  current  draw  and 
regeneration. 

Section  1 .3  describes  the  validation  of  the  lead  acid  battery  model.  Battery  state  of  charge  and 
voltage  predicted  by  the  model  were  compared  with  measured  data  from  the  hybrid  bus.  The 
model  predicted  battery  state  of  charge  compared  within  5%  to  the  measured  data  of  the  same. 
Larger  deviations  were  observed  with  the  comparison  of  the  model  predicted  and  measured 
battery  voltage. 

Section  1 .4  describes  the  validation  of  the  combined  driver-axle-wheel-vehicle  model. 

Finally,  Section  2.0  describes  the  hybrid  bus  model  of  the  22-foot  series  hybrid  shuttle  bus  that 
was  used  to  validate  PATHS. 

Section  3.0  provides  concluding  comments  on  this  validation  report. 
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1.0  INTRODUCTION 


This  document  comprises  the  validation  and  testing  of  the  main  components  of  PATHS.  The 
reader  is  referred  to  the  User's  Guide  and  Technical  Notes  for  a  more  detailed  description  of 
PATHS.  The  validation  includes  the  following  main  models  of  PATHS.  These  are  the  models 
that  are  used  at  a  minimum  to  model  a  general  series  hybrid  vehicle. 

■  APU  model  from  the  engine  library 

■  AC  induction  motor  model  (including  inverter)  from  the  electric  machine  library 

■  Lead  acid  battery  model  from  the  energy  storage  library 

■  Driver-axle-wheel-vehicle  combined  model 

The  data  for  the  validation  of  the  above  components  was  obtained  from  a  22-foot  series  hybrid 
shuttle  bus  developed  by  the  Electricore  consortium.  This  bus  has  a  35-kW,  CNG-fueled 
Auxiliary  Power  Unit  (APU)  with  two  battery  packs,  each  containing  26  lead-acid  batteries.  The 
vehicle  is  propelled  by  two  70-kW  (peak),  AC-induction  wheel  motors,  one  on  each  side  of  the 
rear  wheels.  A  fixed  ratio  speed  reducer  of  20:1  is  integrated  into  the  motor  before  the  output 
to  the  wheels. 

The  hybrid  bus  was  placed  on  a  heavy-duty  chassis  dynamometer  at  Southwest  Research 
Institute  (SwRI)  and  operated  on  two  driving  cycles:  the  CBD  and  EPA  schedule  D  driving 
cycles.  Figure  51  shows  the  CBD  driving  profile,  and  Figure  52  illustrates  the  EPA  schedule  D 
driving  profile. 
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Data  was  collected  over  six  consecutive  EPA  schedule  D  and  ten  CBD-14  driving  cycles 
respectively.  A  CBD-14  driving  profile  consists  of  14  consecutive  CBD  segments  (see  Figure 
50  for  one  CBD  driving  segment).  For  each  of  the  aforementioned  driving  cycles,  the  following 
data  was  collected  on  a  second-by-second  basis: 

■  Bus  voltage  (volts) 

■  Vehicle  speed  (miles  per  hour) 

■  APU  engine  speed  (rpm) 

■  Accelerator  pedal  position  (0- 1 ) 

■  Brake  pedal  position  (0-1) 

■  Current  into  and  out  of  left  motor  inverter  (DC  amps) 

■  Current  into  and  out  of  right  motor  inverter  (DC  Amps) 

■  Current  into  and  out  of  left  battery  pack  (DC  Amps) 

■  Current  into  and  out  of  right  battery  pack  (DC  Amps) 

■  Generator  current  (DC  Amps) 

■  Battery  pack  state  of  charge  (%) 

The  following  sections  describe  the  validation  results  of  the  main  components  of  PATHS. 

1.1  APU  Model  Validation 
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Figure  53.  APU  Mode!  Validation 
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The  APU  model  shown  in  Figure  53  is  of  the  table  look-up  type  and  models  a  35  kW  CNG 
engine,  coupled  to  a  permanent  magnet  generator.  Table  27  describes  the  inputs  and  outputs  of 
the  APU  model. 


Table  27.  APU  MODEL 

INPUTS 

OUTPUTS 

Engine  on/off  state 

Mass  of  fuel  consumed  (kg) 

Desired  electric  power  (kW) 

Generator  current  (DC  Amps) 

Battery  voltage  (volts) 

Engine  speed  (rpm) 

APU  inertia  (kgm2) 

The  APU  model  was  tested  by  providing  collected  test  data  on  bus  voltage  and  generator  current 
(the  electric  power  is  equal  to  the  product  of  the  bus  voltage  and  generator  current)  as  inputs  to 
the  model  and  comparing  the  model  outputs  of  fuel  consumed  and  engine  speed  with  measured 
data  of  the  same.  Figure  54  shows  a  comparison  of  fuel  consumed  as  predicted  by  the  APU 
model  and  measured  fuel  consumption  data.  Figure  55a  and  55b  illustrate  a  comparison  of  the 
APU  emissions  (HC  and  Nox)  predicted  by  the  model  and  the  collected  data. 


Figure  54.  APU  Validation  —  6  EPA  Schedule  D 
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1.2  AC  Induction  Motor  Model 
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Figure  56.  AC  Induction  Motor 


The  AC  induction  motor  model  must  have  several  parameters  defined  prior  to  integration  into 
a  hybrid  vehicle  simulation  model.  The  essential  motor  model  parameters  are  described  in  Table 
28. 


Table  28.  AC  INDUCTION  MOTOR  MODEL  PARAMETERS 

P  (hp): 

Rated  power  of  motor 

V  (volts) 

Rated  operating  voltage 

f  (Hz) 

Rated  operating  frequency 

n(rpm) 

Rated  speed 

CL(hp) 

Core  losses  at  rated  voltage  and  speed 

SLL(hp) 

Stray  load  losses  at  rated  voltage  and  speed 

Rs(O) 

Per  phase  stator  resistance 

Rr(Q) 

Per  phase  rotor  resistance 

Lls(H) 

Per  phase  stator  leakage  inductance 

Llr(H) 

Per  phase  rotor  leakage  resistance 

Lm(H) 

Per  phase  magnetizing  inductance 

n(-) 

Inverter  efficiency 

l(kgm2) 

motor  inertia 

sample_N 

Number  of  samples  saved  to  workspace 

sample_T 

Time  between  each  saved  sample 
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The  above  parameters  must  be  matched  for  an  individual  motor  by  means  of  the  full-throttle 
torque  speed  curve  of  the  motor.  The  full  throttle  and  part  throttle  torque  speed  curves  for  the 
70  kW  AC  induction  motor  are  illustrated  in  Figure  57.  An  iterative  procedure  was  used  to 
match  the  motor  parameters  to  produce  the  torque  speed  characteristics  shown  in  Figure  57. 
Figure  58  shows  the  motor  model  predicted  torque  speed  characteristics  superimposed  on  that 
of  Figure  57.  The  final  motor  parameters  that  were  used  to  simulate  the  hybrid  bus  are  defined 
in  the  model  mask,  illustrated  in  Figure  58. 


Figure  57.  System  Voltage  =312 V  -  Measured  Motor  Data 


Figure  58.  System  Voltage  =  312  V  -  Measured  Motor  Data  vs.  Model 
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Block  name:  motl1__module 

Block  type:  SS  Induction  Motor  Sys  (t 


SS  Induction  Motor/Inverter  System 


Rated  [P[hp),V(V),f(Hz)]: 


[114,312,135] 


Rated  [n(rpm),CL(hp],SLL(hp}]: 


[4000,10,10] 


[Rs,Rr  in  ohms,  Lls,Llr,Lm  in  H]: 


[0.01 983,0.01  235,42. 5E-6,42.5E-6,2.5E~3] 


Inverter  Efficiency: 


.98 


Machine  lnertia(kg*m*m): 


.02821 


[sample_N,sample_T]: 


[20000,0.1] 


Figure  59.  SS  Induction  Motor  Sys 


Table  29.  AC  INDUCTION  MOTOR  MODEL 

INPUTS 

OUTPUTS 

Bus  voltage  (V) 

Motor  current  (DC  Amps) 

Motor  speed  (rpm) 

Motor  torque  (Nm) 

Torque  command  (0-1) 

Motor  inertia  (kgm2) 

The  motor  model  was  validated  by  providing  collected  test  data  on  bus  voltage,  torque 
command  from  driver  accelerator  and  brake  pedal,  and  vehicle  speed  as  inputs  to  the  model  and 
comparing  the  model  outputs  of  current  and  torque  with  measured  data  of  the  same.  Figure  60 
shows  a  comparison  of  the  model  predicted  current  and  measured  data  for  the  first  500  seconds 
of  the  EPA  schedule  D  driving  profile  .  Figure  61  shows  the  same  comparison  for  the  last  500 
seconds  of  the  same  driving  profile. 
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Figure  61.  Motor  Model  Validation  -  EPA  Schedule  D  (last  500  seconds^ 
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1.3  Lead  Acid  Battery  Model 
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Figure  62.  Lead  Acid  Battery  Model 


The  lead  acid  battery  model  parameters  are  presented  in  Table  30. 


Table  30.  LEAD  ACID  BATTERY  MODEL  PARAMETERS 

Number  of  batteries 

(-) 

Open  circuit  voltage  of  a  battery 

(volts) 

Internal  resistance  of  the  battery 

(0) 

Polarization  resistance  of  the  battery  (Q) 

(Q) 

Battery  capacitance 

(F) 

Terminal  battery  resistance 

(Cl) 

Terminal  battery  capacitance 

(fi) 

Discharge  capacity  at  a  5  hour  rate  and  30° 
Celsius 

(A-hr) 

Discharge  current  at  a  5  hour  rate  and  30° 
Celsius 

(A) 

Charge  efficiency 

(0-1) 

Discharge  efficiency 

(0-1) 

Operating  temperature 

(°C) 

Initial  polarization  voltage 

(V) 

Initial  state  of  charge 

(0-100) 

The  hybrid  bus  is  equipped  with  two  battery  packs,  each  containing  26  batteries.  The  above 
parameters  needed  for  the  model  were  obtained  from  the  battery  manufacturer  and  are  entered 
into  the  model  mask,  illustrated  in  Figure  63.  Note  the  mask  represents  a  battery  pack  containing 
26  batteries. 
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Figure  63.  Lead-Acid  Battery  -  ESS 


The  input  and  outputs  of  the  battery  model  are  described  in  Table  31. 


Table  31.  BATTERY  MODEL 


Inputs 


Charge/discharge  current  (DC  Amps) 


Outputs 


Battery  voltage  (volts) 
State  of  charge  (0-1) 


The  battery  model  was  validated  by  supplying  to  the  battery  model  input  the  measured 
charge/discharge  current  and  comparing  the  model  predicted  battery  pack  state  of  charge  and 
voltage  with  measured  data.  Figure  64  shows  a  comparison  of  the  battery  state  of  charge  as 
predicted  by  the  model  and  measured  on  the  bus.  Figure  65  shows  a  comparison  of  the  battery 
voltage  as  predicted  by  the  model  and  measured  on  the  bus  using  the  EPA  Schedule  D  test. 
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Battery  State  of  Charge  (%) 


1.4  Driver-Axle-Wheel-Vehicle  Combined  Model 
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Figure  66.  Driver-Axle- Wheel-Vehicle  Combined  Model 


This  combined  model  consists  of  four  component  models.  The  nature  of  the  available  test  data 
required  us  to  combine  these  four  components  into  a  single  model  for  testing  and  validation 
purposes.  The  Tables  32-35  summarize  the  inputs  and  outputs  of  each  of  the  component  models. 


Table  32.  DRIVER  MODEL 

INPUTS 

OUTPUTS 

Desired  vehicle  speed  (m/s) 

Power  demand  (0-1) 

Actual  vehicle  speed  (m/s) 

Table  33.  AXLE  MODEL 

INPUTS 

OUTPUTS 

Input  torque  (Nm) 

Output  torque  (Nm) 

Output  shaft  speed  (rpm) 

Input  shaft  speed  (rpm) 

Input  inertia  (kgm2) 

Equivalent  output  inertia  (kgm2) 
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|  Table  34.  WHEEL  MODEL 

INPUTS 

OUTPUTS 

Maximum  traction  force  (N) 

Rotational  wheel  speed  (rpm) 

Power  demand  from  driver  (0-1) 

Vehicle  speed  (m/s) 

Torque  at  input  (Nm) 

Vehicle  acceleration  (m/s2) 

Total  resistive  load  (N) 

Wheel  moment  of  inertia  (kgm2) 

Total  vehicle  mass  (kg) 

Service  brake  force  (N) 

Equivalent  inertia  (kgm2)  from  axle 


Table  35.  VEHICLE  MODEL 

INPUTS 

OUTPUTS 

Vehicle  speed  (m/s) 

Total  load  on  vehicle  (N) 

Gradient  (radians) 

Maximum  traction  force  (N) 

Fuel  mass  consumed  (kg) 

Vehicle  speed  (m/s) 

Vehicle  acceleration  (m/s2) 

Total  vehicle  mass  (kg) 

Service  brake  force  (N) 


that  predicted  by 


2.0  Hybrid  Bus  Model 
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Figure  68.  Hybrid  Bus  Model 


The  subcomponents  described  in  Sections  1.1  to  1.4  were  combined  to  model  a  22-foot  series 
hybrid  shuttle  bus  developed  by  the  Electricore  consortium.  As  described  in  Section  1,  this  bus 
is  equipped  with  two  70-kW  wheel,  AC-induction  motors,  a  35-kW  CNG-fueled  APU,  and  two 
battery  packs,  each  containing  26  lead-acid  batteries.  The  weight  of  the  hybrid  bus  is 
approximately  14,600  pounds  (empty  weight  plus  an  estimated  10  passengers).  The  inertia 
weight  for  the  chassis  dynamometer  was  computed  from  this  estimate  of  the  bus  weight.  Based 
on  the  frontal  area  and  drag  coefficient,  the  road  load  was  estimated  at  52  hp  at  a  speed  of  50 
mph.  The  rolling  resistance  was  computed  to  be  17  hp  and  the  aerodynamic  losses  were 
computed  as  35  hp.  The  power  management  controller  turned  the  APU  on  when  the  battery  state 
of  charge  fell  below  75%,  and  turned  the  APU  off  when  the  battery  state  of  charge  increased 
above  85%.  The  electrical  power  requested  from  the  APU  was  a  function  of  state  of  charge,  that 
is  user  modifiable.  Figure  69  shows  the  model  of  the  hybrid  bus  using  the  components  described 
in  Sections  1.1  to  1.4. 
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Figure  69  shows  a  comparison  of  the  battery  state  of  charge,  while  the  hybrid  bus  was  operated 
over  6  EPA-Schedule  D  driving  profiles.  Figure  70  shows  the  relative  error  of  the  model 
predicted  battery  voltage  with  respect  to  the  measured  data  during  operation  over  6  EPA- 
Schedule  D  driving  profiles.  Figure  71  compares  the  vehicle  speed  predicted  by  the  model  with 
data  measured  on  the  dynamometer  test  stand. 


Figure  69.  Hybrid  Bus  Model  Validation  -  EPA  Schedule  D 
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Vehicle  Speed  (m/s)  Relative  Voltage  Error  with  Respect  to  Measured  Data  (%) 
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3.0  CONCLUSIONS 


This  document  describes  the  first  attempt  to  validate  PATHS.  The  validation  exercise  was 
limited  to  key  components  from  the  many  libraries  of  PATHS.  The  validation  was  performed 
by  comparing  the  model  predictions  with  data  collected  from  a  hybrid  bus  operating  on  EPA  and 
CBD  driving  profiles.  The  models  that  were  validated  appear  to  agree  with  the  test  data 
reasonably  well.  The  results  that  are  described  in  this  report  lend  significant  credibility  in 
applying  PATHS  to  designing  and  building  hybrid  vehicles. 

It  is  apparent  that  some  models  compare  very  well  with  measured  data,  while  others  do  not.  For 
example,  the  APU  model  was  able  to  predict  fuel  consumption  to  within  1%  of  the  measured 
data,  while  predictions  of  the  HC  emissions  were  off  by  almost  a  100%.  Although  the  predicted 
battery  state  of  charge  compared  quite  well  to  the  measured  data  (within  less  than  5%),  the 
predicted  battery  voltage  was  sometimes  much  larger  than  what  was  measured.  However,  it  was 
remarkable  to  observe  the  trends  were  duplicated  throughout  the  driving  profile. 

Ideally,  changes  to  the  model  would  have  to  be  made  after  performing  the  described  validation 
exercise  in  order  to  refine  and  develop  a  better  model.  There  were  no  funds  remaining  for  this 
task,  but  we  hope  that  funding  will  be  available  for  this  important  task  in  the  future. 

The  validation  results  have  provided  the  much-needed  confidence  to  use  PATHS  in  designing 
and  developing  a  hybrid  or  pure  electric  vehicle.  We  hope  that  it  will  find  more  widespread 
usage  throughout  the  electric/hybrid  vehicle  industry  and  become  an  accepted  design  tool. 
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